esei_Sphere.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluEllipticIntegralsSphere.h
// author(s):      Jan Novak, Lukasz Kaczmarczyk
// 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_SPHERE_H
#define MUMECH_ESHELBY_SOLU_ELLIPTIC_INTEGRALS_SPHERE_H

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


namespace mumech {

class eshelbySoluEllipticIntegralsSphere : public eshelbySoluEllipticIntegrals
{
 public:
  eshelbySoluEllipticIntegralsSphere (const InclusionRecord3D *i) : eshelbySoluEllipticIntegrals (i) { }
  virtual ~eshelbySoluEllipticIntegralsSphere () {}
  
  //******************************************************************************************************
  //Function gives a component of the Eshelby tensor
  double eshelbyTensCompUniformField( EshelbyTensComponent flag, double nu, double a );
  //******************************************************************************************************
  //Function initialize an 1D array by Eshelby tensor components
  void eshelbyTensUniformFieldSphere( double eshTens[12], double nu, double a );
  //******************************************************************************************************
  //Function gives the values of all elliptic integrals
  void giveEllipticIntegrals (double J[13], double lambda, bool intpoint); // Overloaded. Should not depend on anyting else than this.

  // void giveEllipticIntegrals (double J[13], const double a[3], double lambda, const double loc_x[3]);
  // //******************************************************************************************************
  // //Function gives the values of all elliptic Ferers-Dysons integrals
  // void give_Ii_and_Iij( double J[13], double mult, double norm_x[5], double lambda, double aa );
  // //******************************************************************************************************
  // //Function gives the solution of the integral I prefix 'e' means "from euclidean space"
  // double eI( double mult, double norm_x[5], const double a[3], double lambda , bool intpoint);
  // //******************************************************************************************************
  // //Function gives the solution of the integral Ii
  // double Ii( double mult, double norm_x[5], double lambda, bool intpoint );
  // //******************************************************************************************************
  // //Function gives the solution of the integral Iij
  // double Iij(double mult, double norm_x[5], double lambda, double aa , bool intpoint);
  // //******************************************************************************************************
  
  
#ifndef DEBUG
 protected:
#endif
  
  //******************************************************************************************************
  //Function gives the values of all elliptic integrals' derivatives
  void giveDerivativesOfEllipticIntegrals(Point * point, bool intpoint);
  
  //Ferers-Dysons integrals:

  //Derivatives of Ferers-Dysons integrals:
  //******************************************************************************************************
  //Function gives the first derivative of the integral Ii
  double Ii_i( double mult, double norm_x[5], double x[3], derivativeDirection direction );
  //******************************************************************************************************
  //Function gives the second derivative of the integral Ii
  double Ii_ij( double mult, double norm_x[5], double x[3], derivativeDirection direction_1,
        derivativeDirection direction_2 );
  //******************************************************************************************************
  //Function gives the first derivative of the integral Iij
  double Iij_i( double mult, double norm_x[5], double x[3], derivativeDirection direction );
  //******************************************************************************************************
  //Function gives the second derivative of the integral Iij
  double Iij_ij( double mult, double norm_x[5], double x[3], derivativeDirection direction_1,
         derivativeDirection direction_2 );
  //******************************************************************************************************
  //Function gives all first derivatives of Ferers-Dysons' integrals Ii
  void give_dIi( double dJi[13], double mult, double norm_x[5], double x[3] );
  //******************************************************************************************************
  //Function gives all first derivatives of Ferers-Dysons' integrals Iij
  void give_dIij( double dJij[27], double mult, double norm_x[5], double x[3] );
  //******************************************************************************************************
  void give_ddIi( double ddJi[27], double mult, double norm_x[5], double x[3] );
  //******************************************************************************************************
  //Function gives all second derivatives of Ferers-Dysons' integrals Iij
  void give_ddIij( double ddJij[81], double mult, double norm_x[5], double x[3] );
  //******************************************************************************************************

  //Complex numerical derivatives of Ferers-Dysons integrals:
  //---------------------------------------------------------

  //******************************************************************************************************
  //Function gives the first derivative of the integral Ii by using the complex differentiation
  std::complex<double> cIi_i( double mult, std::complex<double> x[3], derivativeDirection direction );
  //******************************************************************************************************
  //Function gives the first derivative of the integral Iij by using the complex differentiation
  std::complex<double> cIij_i( double mult, std::complex<double> x[3], derivativeDirection direction );
  //******************************************************************************************************
  //Function gives all first derivatives of Ferers-Dysons' integrals Ii by using the complex
  //differentiation
  void give_cdIi( std::complex<double> dJi[13], double mult, double x[3] );
  //******************************************************************************************************
  //Function gives all first derivatives of Ferers-Dysons' integrals Iij by using the complex
  //differentiation
  void give_cdIij( std::complex<double> dJij[13], double mult, double x[3] );
  //******************************************************************************************************
  //Function updates vectors cx according the derivative direction
  void update_cx( std::complex<double> cx[3], derivativeDirection direction );
  //******************************************************************************************************
  //Function gives all second derivatives of Ferers-Dysons' integrals Ii by using the complex
  //differentiation
  void give_cddIi( double ddJi[27], std::complex<double> dJi[13] );
  //******************************************************************************************************
  //Function gives all second derivatives of Ferers-Dysons' integrals Iij by using the complex
  //differentiation
  void give_cddIij( double ddJij[81], std::complex<double> dJij[27] );
  //******************************************************************************************************

  double getLambda(const double a[3], double x1, double x2, double x3);
};

} // end of namespace mumech

#endif

/*end of file*/