esei_Cylinder.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluEllipticIntegralsCylinder.h
// author(s):      Jan Novak, Tomas Janda
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 This program is distributed in the hope that it will be useful,
//                 but WITHOUT ANY WARRANTY; without even the implied warranty of
//                 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//                 GNU Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************
#ifndef MUMECH_ESHELBY_SOLU_ELLIPTIC_INTEGRALS_CYLINDER_H
#define MUMECH_ESHELBY_SOLU_ELLIPTIC_INTEGRALS_CYLINDER_H

#include "types.h"
#include "esei.h"
#include <complex>


namespace mumech {

class eshelbySoluEllipticIntegralsCylinder : public eshelbySoluEllipticIntegrals
{
 private:
  bool log;
  
 public:
  eshelbySoluEllipticIntegralsCylinder (const InclusionRecord3D *i) : eshelbySoluEllipticIntegrals (i) { log = true; }
  virtual ~eshelbySoluEllipticIntegralsCylinder () {}
  
  //Function gives the the values of all elliptic integrals
  void giveEllipticIntegrals(double J[13], double lambda, bool intpoint);
  //Function gives the derivatives of Ferers-Dysons integrals
  void giveDerivativesOfEllipticIntegrals(Point * point, bool intpoint);
  // Returns lambda for a given point (x1, x2, x3)
  double getLambda(const double a[3], double x1, double x2, double x3);
  
 private:
  // Function returning integrals Ii and Iij for cylinder (Mura, p.80)
  // For internal fields: lambda = 0
  // For external fields: lambda != 0
  
  // I2 = I3 due to the symmetry of cylinder
  double I2(double a2, double lambda);
  // I22 = I23 = I33 due to the symmetry of cylinder
  double I22(double a1, double lambda);
  
};

} // end of namespace mumech

#endif

/*end of file*/