esei_Penny.cpp

Go to the documentation of this file.

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

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <cmath>
#include <iostream>

#include "macros.h"
#include "esei_Penny.h"
#include "eshelbySoluLambda.h"
#include "legendreIntegrals.h"
#include "elementaryFunctions.h"
#include "inclusion.h"

namespace mumech {
  
// Function gives the values of Ferers-Dysons integrals
// Last edit:   18-Oct-2013
// ATTENTION! The last element in the J array is not I33 but a3^2*I33.
void eshelbySoluEllipticIntegralsPenny :: giveEllipticIntegrals (double J[13], double lambda, bool intpoint)
{
  // Penny shape: a1 = a2 >> a3
  if (I->a[0] != I->a[1])
    {
      _warningg( "Values of semiaxes do not fit the condition for penny shape (a1 = a2). " );
      std::cout << "a1, a2, a3 = " << I->a[0] << ", " << I->a[1] << ", " << I->a[2];
    }
  
  double a = I->a[0];
  
  CleanVector(J, 13);
  // termitovo: toto by se melo predelat na detekci podle intpoint, ale nevim jak to je tady s logikou bodu
  if (lambda > MACHINE_EPS) // External point
    return;
  
  // Internal point
  J[0] = 0.;
  
  J[1] = 0.0; // I1
  J[2] = 0.0; // I2
  J[3] = _4PI_; // I3
  
  J[4] = 0.0; // I11
  J[5] = 0.0; // I12
  J[6] = _4PI_ / SQR(a); // I13
  
  J[7] = 0.0; // I21
  J[8] = 0.0; // I22
  J[9] = _4PI_ / SQR(a); // I23
  
  J[10] = _4PI_ / SQR(a); // I31
  J[11] = _4PI_ / SQR(a); // I32
  J[12] = _4PI_ / 3.0; // a3^2 * I33s
}

// Function gives the derivatives of Ferers-Dysons integrals
// Last edit:   11-Sep-2013
// Flat shapes has zero drtivatives
void eshelbySoluEllipticIntegralsPenny :: giveDerivativesOfEllipticIntegrals (Point *point, bool intpoint)
{
  CleanVector( point->dJi, 9 );
  CleanVector( point->dJij, 27 );
  CleanVector( point->ddJi, 27 );
  CleanVector( point->ddJij, 81 );
}

}