esei.cpp

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluEllipticalIntegrals.cpp
// 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.
//********************************************************************************************************

// #include <stdio.h>
// #include <stdlib.h>
// #include <math.h>
// #include "macros.h"
#include "esei.h"
// #include "esei_Ellipsoid.h"
// #include "esei_Sphere.h"
// #include "esei_OblateSpheroid.h"
// #include "esei_ProlateSpheroid.h"
// #include "esei_Penny.h"
// #include "esei_FlatEllipsoid.h"
// #include "esei_EllipticCylinder.h"
// #include "esei_Cylinder.h"
#include "inclusion.h"


namespace mumech {

//********************************************************************************************************
//                                            PUBLIC FUNCTIONS
//********************************************************************************************************


// Computes lambda of each perturbated point.
void eshelbySoluEllipticIntegrals :: getPerturbatedLambdas (double *lambdas, const double loc_x[3])
{
  // Define coordinates of current point
  double x = loc_x[0];
  double y = loc_x[1];
  double z = loc_x[2];
  double xp = x + I->ndiff_1;
  double xm = x - I->ndiff_1;
  double yp = y + I->ndiff_1;
  double ym = y - I->ndiff_1;
  double zp = z + I->ndiff_1;
  double zm = z - I->ndiff_1;
  
  // Fill the list
  lambdas[_CENTER_] = getLambda(I->a, x, y, z);
  // X
  lambdas[_Xp_] = getLambda(I->a, xp, y, z);
  lambdas[_Xm_] = getLambda(I->a, xm, y, z);
  // Y
  lambdas[_Yp_] = getLambda(I->a, x, yp, z);
  lambdas[_Ym_] = getLambda(I->a, x, ym, z);
  // Z
  lambdas[_Zp_] = getLambda(I->a, x, y, zp);
  lambdas[_Zm_] = getLambda(I->a, x, y, zm);
  // X, Y
  lambdas[_XpYp_] = getLambda(I->a, xp, yp, z);
  lambdas[_XpYm_] = getLambda(I->a, xp, ym, z);
  lambdas[_XmYp_] = getLambda(I->a, xm, yp, z);
  lambdas[_XmYm_] = getLambda(I->a, xm, ym, z);
  // X, Z
  lambdas[_XpZp_] = getLambda(I->a, xp, y, zp);
  lambdas[_XpZm_] = getLambda(I->a, xp, y, zm);
  lambdas[_XmZp_] = getLambda(I->a, xm, y, zp);
  lambdas[_XmZm_] = getLambda(I->a, xm, y, zm);
  // Y, Z
  lambdas[_YpZp_] = getLambda(I->a, x, yp, zp);
  lambdas[_YpZm_] = getLambda(I->a, x, yp, zm);
  lambdas[_YmZp_] = getLambda(I->a, x, ym, zp);
  lambdas[_YmZm_] = getLambda(I->a, x, ym, zm);
}


} // end of namespace mumech

/*end of file*/