selfBalanceAlgorithm.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           selfBalanceAlgorithm.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_SELF_BALANCE_ALGORITHMS_H
#define MUMECH_SELF_BALANCE_ALGORITHMS_H

//included libraries
#include "inclusion.h"
#include "types.h"
#include "matrix.h"


namespace mumech {

//selfBalanceAlgorithm class declaration
class selfBalanceAlgorithm
{
 public:
  selfBalanceAlgorithm(void) { }
  ~selfBalanceAlgorithm( ) { }
  
  //******************************************************************************************************
  // Function gives the total transformation strain field in inclusion centroids regarding
  // the presence of other adjacent inlusion.
  static void totalEigStrainInInclCentroidsUpdate (Problem *p);
  
  
 protected:

  //******************************************************************************************************
  // Protected variables
  
  //******************************************************************************************************
  // Function gives the quadratic norm of multiple eigenstrain vectors calculated
  // from strain perturbation values
  // double giveQuadNormMultEigstrain( Point * point, double * oldStrain, int n );

  //******************************************************************************************************
  // Function copies strain perturbatios in inclusion origins into an auxiliary vector:
  //void copyPertStrainsToAuxVect(Point *point, double * auxVect, int n );
  
  //******************************************************************************************************
  // Function copies total global eigenstrains from InclRecord to an auxiliary vector:
  // auxVect  '@n x 6' ... this is number of vector entries
  static void copyGlobTotalEigStrainsToAuxVector (InclusionRecord3D **inclRec, double *auxVect, int n);
  
  //******************************************************************************************************
  // Function gives the quadratic norm of multiple eigenstrain vectors calculated
  // from total global eigenstrain values.
  static double giveQuadNormMultEigstrain (InclusionRecord3D **inclRec, double *oldStrain, int n);
  
  //******************************************************************************************************
  // Function copies total global eigenstrain in inclusion origins into an auxiliary vector:
  static void copyGlobTotEigStrainsToAuxVect (InclusionRecord3D **inclRec, double * auxVect, int n );
  
  //******************************************************************************************************
  // Function updates total local eigenstrains from their total global counterparts
  static void localEigStrainUpdateTotal (mumech::InclusionRecord3D **inclRec, int n );
  
  //******************************************************************************************************
  // Function prints resulting total transformation strains in inclusion centroids
  static void printTotEigenStrains (InclusionRecord3D **inclRec, int n, const char * notice );
  
  
}; // end of class declaration

} // end of namespace mumech

#endif

/*end of file*/