esei.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluEllipticIntegrals.h
// author(s):      Jan Novak
// 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_H
#define MUMECH_ESHELBY_SOLU_ELLIPTIC_INTEGRALS_H

#include "types.h"
#include "matrix.h"

namespace mumech {

class InclusionRecord3D;

class eshelbySoluEllipticIntegrals
{
 protected:
  const InclusionRecord3D *I;
  
 public:
  eshelbySoluEllipticIntegrals (const InclusionRecord3D *i) { I = i; }
  virtual ~eshelbySoluEllipticIntegrals () {}
  
  virtual void giveEllipticIntegrals (double integrArray[13], double lambda, bool intpoint) = 0;
  
  virtual void giveDerivativesOfEllipticIntegrals (Point *point, bool intpoint) = 0;
  
  virtual double getLambda (const double a[3], double x1, double x2, double x3) = 0;

  void getPerturbatedLambdas (double *lambdas, const double loc_x[3]);


};

} // end of namespace mumech

#endif

/*end of file*/