esei_OblateSpheroid.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluEllipticIntegralsOblateSpheroid.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_OBLATE_SPHEROID_H
#define MUMECH_ESHELBY_SOLU_ELLIPTIC_INTEGRALS_OBLATE_SPHEROID_H

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


namespace mumech {

class eshelbySoluEllipticIntegralsOblateSpheroid : public eshelbySoluEllipticIntegrals
{
 public:
  eshelbySoluEllipticIntegralsOblateSpheroid (const InclusionRecord3D *i) : eshelbySoluEllipticIntegrals (i) { }
  virtual ~eshelbySoluEllipticIntegralsOblateSpheroid () {}
  
  //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);
  
 private:
  // Function returning integrals Ii and Iij for oblate spheroid (Mura, p.84)
  // For internal fields: lambda = 0
  // For external fields: lambda != 0
  double eI(double a[3], double k, double Theta, double mult);
  double I1(double a1, double a3, double lambda, bool intpoint);
  double I2(double I1);
  double I3(double a1, double a3, double I1, double lambda, bool intpoint);
  double I11(double I12);
  double I12(double a1, double a3, double I13, double lambda, bool intpoint);
  double I13(double a1, double a3, double I1, double I3);
  double I22(double I12);
  double I23(double I13);
  double I33(double a1, double a3, double I13, double lambda, bool intpoint);
  
  // Helper functions
  double get_b(double a1, double a3, double lambda); // Auxiliary number b (Mura, p.94)
  double get_d(double a1, double a3, double lambda); // Auxiliary number d (Mura, p.94)
  double get_delta_lambda(double a1, double a3, double lambda); // Delta lambda (Mura, p.93)
  double getLambda(const double a[3], double x1, double x2, double x3); // Returns lambda for a given point (x1, x2, x3)
};

} // end of namespace mumech

#endif

/*end of file*/