-// Copyright (C) 2007-2015 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2022 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
#include <math.h>
#include <algorithm>
#include <sstream>
+#include <cmath>
using namespace INTERP_KERNEL;
+const double GaussInfo::SEG2A_REF[2]={-1.0, 1.0};
+
+const double GaussInfo::SEG2B_REF[2]={0., 1.0};
+
+const double GaussInfo::SEG3_REF[3]={-1.0, 1.0, 0.0};
+
+const double GaussInfo::TRIA3A_REF[6]={-1.0, 1.0, -1.0, -1.0, 1.0, -1.0};
+
+const double GaussInfo::TRIA3B_REF[6]={0.0, 0.0, 1.0, 0.0, 0.0, 1.0};
+
+const double GaussInfo::TRIA6A_REF[12]={-1.0, 1.0, -1.0, -1.0, 1.0, -1.0, -1.0, 1.0, 0.0, -1.0, 0.0, 0.0};
+
+const double GaussInfo::TRIA6B_REF[12]={0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.5, 0.0, 0.5, 0.5, 0.0, 0.5};
+
+const double GaussInfo::TRIA7A_REF[14]={0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.5, 0.0, 0.5, 0.5, 0.0, 0.5, 0.3333333333333333, 0.3333333333333333};
+
+const double GaussInfo::QUAD4A_REF[8]={-1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0};
+
+const double GaussInfo::QUAD4B_REF[8]={-1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0};
+
+const double GaussInfo::QUAD8A_REF[16]={-1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 1.0};
+
+const double GaussInfo::QUAD8B_REF[16]={-1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0, 0.0, -1.0, 1.0, 0.0, 0.0, 1.0, -1.0, 0.0};
+
+const double GaussInfo::QUAD9A_REF[18]={-1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0, 0.0, -1.0, 1.0, 0.0, 0.0, 1.0, -1.0, 0.0, 0.0, 0.0};
+
+const double GaussInfo::TETRA4A_REF[12]={0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0};
+
+const double GaussInfo::TETRA4B_REF[12]={0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0};
+
+const double GaussInfo::TETRA10A_REF[30]={0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.5, 0.5, 0.0, 0.0, 0.5, 0.0, 0.5, 0.0, 0.5, 0.5, 0.0, 0.5, 0.0, 0.5, 0.5, 0.0, 0.0};
+
+const double GaussInfo::TETRA10B_REF[30]={0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0, 0.0, 0.5, 0.0, 0.0, 0.0, 0.5, 0.0, 0.5, 0.5, 0.5, 0.5, 0.0, 0.5, 0.0, 0.0, 0.5, 0.0, 0.5};
+
+const double GaussInfo::PYRA5A_REF[15]={1.0, 0.0, 0.0, 0.0, 1.0, 0.0, -1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0};
+
+const double GaussInfo::PYRA5B_REF[15]={1.0, 0.0, 0.0, 0.0, -1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0};
+
+const double GaussInfo::PYRA13A_REF[39]={1.0, 0.0, 0.0, 0.0, 1.0, 0.0, -1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.5, 0.5, 0.0, -0.5, 0.5, 0.0, -0.5, -0.5, 0.0, 0.5, -0.5, 0.0, 0.5, 0.0, 0.5, 0.0, 0.5, 0.5, -0.5, 0.0, 0.5, 0.0, -0.5, 0.5};
+
+const double GaussInfo::PYRA13B_REF[39]={1.0, 0.0, 0.0, 0.0, -1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.5, -0.5, 0.0, -0.5, -0.5, 0.0, -0.5, 0.5, 0.0, 0.5, 0.5, 0.0, 0.5, 0.0, 0.5, 0.0, -0.5, 0.5, -0.5, 0.0, 0.5, 0.0, 0.5, 0.5};
+
+const double GaussInfo::PENTA6A_REF[18]={-1.0, 1.0, 0.0, -1.0, -0.0, 1.0, -1.0, 0.0, 0.0, 1.0, 1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 0.0, 0.0};
+
+const double GaussInfo::PENTA6B_REF[18]={-1.0, 1.0, 0.0, -1.0, 0.0, 0.0, -1.0, -0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 1.0};
+
+const double GaussInfo::PENTA15A_REF[45]={-1.0, 1.0, 0.0, -1.0, -0.0, 1.0, -1.0, 0.0, 0.0, 1.0, 1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 0.0, 0.0, -1.0, 0.5, 0.5, -1.0, 0.0, 0.5, -1.0, 0.5, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.5, 0.5, 1.0, 0.0, 0.5, 1.0, 0.5, 0.0};
+
+const double GaussInfo::PENTA15B_REF[45]={-1.0, 1.0, 0.0, -1.0, 0.0, 0.0, -1.0, -0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 1.0, -1.0, 0.5, 0.0, -1.0, 0.0, 0.5, -1.0, 0.5, 0.5, 1.0, 0.5, 0.0, 1.0, 0.0, 0.5, 1.0, 0.5, 0.5, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0};
+
+const double GaussInfo::PENTA18A_REF[54]={-1.0, 1.0, 0.0, -1.0, -0.0, 1.0, -1.0, 0.0, 0.0, 1.0, 1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 0.0, 0.0, -1.0, 0.5, 0.5, -1.0, 0.0, 0.5, -1.0, 0.5, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.5, 0.5, 1.0, 0.0, 0.5, 1.0, 0.5, 0.0, 0.0, 0.5, 0.5, 0.0, 0.0, 0.5, 0.0, 0.5, 0.0};
+
+const double GaussInfo::PENTA18B_REF[54]={-1.0, 1.0, 0.0, -1.0, 0.0, 0.0, -1.0, -0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, 1.0, -1.0, 0.5, 0.0, -1.0, 0.0, 0.5, -1.0, 0.5, 0.5, 1.0, 0.5, 0.0, 1.0, 0.0, 0.5, 1.0, 0.5, 0.5, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.5, 0.0, 0.0, 0.0, 0.5, 0.0, 0.5, 0.5};
+
+const double GaussInfo::HEXA8A_REF[24]={-1.0, -1.0, -1.0, 1.0, -1.0, -1.0, 1.0, 1.0, -1.0, -1.0, 1.0, -1.0, -1.0, -1.0, 1.0, 1.0, -1.0, 1.0, 1.0, 1.0, 1.0, -1.0, 1.0, 1.0};
+
+const double GaussInfo::HEXA8B_REF[24]={-1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0, -1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, -1.0, 1.0};
+
+const double GaussInfo::HEXA20A_REF[60]={-1.0, -1.0, -1.0, 1.0, -1.0, -1.0, 1.0, 1.0, -1.0, -1.0, 1.0, -1.0, -1.0, -1.0, 1.0, 1.0, -1.0, 1.0, 1.0, 1.0, 1.0, -1.0, 1.0, 1.0, 0.0, -1.0, -1.0, 1.0, 0.0, -1.0, 0.0, 1.0, -1.0, -1.0, 0.0, -1.0, -1.0, -1.0, 0.0, 1.0, -1.0, 0.0, 1.0, 1.0, 0.0, -1.0, 1.0, 0.0, 0.0, -1.0, 1.0, 1.0, 0.0, 1.0, 0.0, 1.0, 1.0, -1.0, 0.0, 1.0};
+
+const double GaussInfo::HEXA20B_REF[60]={-1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0, -1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 0.0, -1.0, 0.0, 1.0, -1.0, 1.0, 0.0, -1.0, 0.0, -1.0, -1.0, -1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 0.0, -1.0, 1.0, -1.0, -1.0, 0.0, -1.0, 1.0, 0.0, 1.0, 1.0, 0.0, 1.0, -1.0, 0.0};
+
+const double GaussInfo::HEXA27A_REF[81]={-1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, -1.0, 1.0, -1.0, -1.0, -1.0, -1.0, 1.0, -1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 0.0, -1.0, 0.0, 1.0, -1.0, 1.0, 0.0, -1.0, 0.0, -1.0, -1.0, -1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 1.0, 0.0, 1.0, 0.0, -1.0, 1.0, -1.0, -1.0, 0.0, -1.0, 1.0, 0.0, 1.0, 1.0, 0.0, 1.0, -1.0, 0.0, 0.0, 0.0, -1.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0};
+
//Define common part of the code in the MACRO
//---------------------------------------------------------------
#define LOCAL_COORD_MACRO_BEGIN \
}
//---------------------------------------------------------------
-#define SHAPE_FUN_MACRO_BEGIN \
+#define SHAPE_FUN_MACRO_BEGIN \
for( int gaussId = 0 ; gaussId < _my_nb_gauss ; gaussId++ ) \
- { \
- double* funValue = &_my_function_value[ gaussId * _my_nb_ref ]; \
+ { \
+ double* funValue = &_my_function_value[ gaussId * _my_nb_ref ]; \
const double* gc = &_my_gauss_coord[ gaussId * getGaussCoordDim() ];
//---------------------------------------------------------------
-#define SHAPE_FUN_MACRO_END \
+#define SHAPE_FUN_MACRO_END \
}
+#define DEV_SHAPE_FUN_MACRO_BEGIN \
+ for( int gaussId = 0 ; gaussId < _my_nb_gauss ; gaussId++ ) \
+ { \
+ double *devFunValue = _my_derivative_func_value.data() + gaussId * getReferenceCoordDim() * _my_nb_ref; \
+ const double *gc = _my_gauss_coord.data() + gaussId * getGaussCoordDim();
+
+#define DEV_SHAPE_FUN_MACRO_END \
+ }
+
#define CHECK_MACRO \
if( ! aSatify ) \
{ \
_my_nb_ref(theNbRef),
_my_reference_coord(theReferenceCoord)
{
-
//Allocate shape function values
_my_function_value.resize( _my_nb_gauss * _my_nb_ref );
+ _my_derivative_func_value.resize( _my_nb_gauss * _my_nb_ref * getReferenceCoordDim() );
}
/*!
{
if( _my_nb_gauss )
{
- return _my_gauss_coord.size()/_my_nb_gauss;
+ return (int)_my_gauss_coord.size()/_my_nb_gauss;
}
else
{
{
if( _my_nb_ref )
{
- return _my_reference_coord.size()/_my_nb_ref;
+ return (int)(_my_reference_coord.size()/_my_nb_ref);
}
else
{
return _my_nb_ref;
}
+GaussInfo GaussInfo::convertToLinear() const
+{
+ switch(_my_geometry)
+ {
+ case NORM_SEG3:
+ {
+ std::vector<double> a(SEG3_REF,SEG3_REF+3);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(SEG2A_REF,SEG2A_REF+2);
+ return GaussInfo(NORM_SEG2,_my_gauss_coord,getNbGauss(),c,2);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for SEG3 !");
+ }
+ case NORM_TRI6:
+ {
+ std::vector<double> a(TRIA6A_REF,TRIA6A_REF+12),b(TRIA6B_REF,TRIA6B_REF+12);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(TRIA3A_REF,TRIA3A_REF+6);
+ return GaussInfo(NORM_TRI3,_my_gauss_coord,getNbGauss(),c,3);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(TRIA3B_REF,TRIA3B_REF+6);
+ return GaussInfo(NORM_TRI3,_my_gauss_coord,getNbGauss(),c,3);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for TRI6 !");
+ }
+ case NORM_TRI7:
+ {
+ std::vector<double> a(TRIA7A_REF,TRIA7A_REF+14);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(TRIA3B_REF,TRIA3B_REF+6);
+ return GaussInfo(NORM_TRI3,_my_gauss_coord,getNbGauss(),c,3);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for TRI7 !");
+ }
+ case NORM_QUAD8:
+ {
+ std::vector<double> a(QUAD8A_REF,QUAD8A_REF+16),b(QUAD8B_REF,QUAD8B_REF+16);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(QUAD4A_REF,QUAD4A_REF+8);
+ return GaussInfo(NORM_QUAD4,_my_gauss_coord,getNbGauss(),c,4);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(QUAD4B_REF,QUAD4B_REF+8);
+ return GaussInfo(NORM_QUAD4,_my_gauss_coord,getNbGauss(),c,4);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for QUAD8 !");
+ }
+ case NORM_QUAD9:
+ {
+ std::vector<double> a(QUAD9A_REF,QUAD9A_REF+18);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(QUAD4B_REF,QUAD4B_REF+8);
+ return GaussInfo(NORM_QUAD4,_my_gauss_coord,getNbGauss(),c,4);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for QUAD9 !");
+ }
+ case NORM_TETRA10:
+ {
+ std::vector<double> a(TETRA10A_REF,TETRA10A_REF+30),b(TETRA10B_REF,TETRA10B_REF+30);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(TETRA4A_REF,TETRA4A_REF+12);
+ return GaussInfo(NORM_TETRA4,_my_gauss_coord,getNbGauss(),c,4);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(TETRA4B_REF,TETRA4B_REF+12);
+ return GaussInfo(NORM_TETRA4,_my_gauss_coord,getNbGauss(),c,4);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for TETRA10 !");
+ }
+ case NORM_PYRA13:
+ {
+ std::vector<double> a(PYRA13A_REF,PYRA13A_REF+39),b(PYRA13B_REF,PYRA13B_REF+39);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(PYRA5A_REF,PYRA5A_REF+15);
+ return GaussInfo(NORM_PYRA5,_my_gauss_coord,getNbGauss(),c,5);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(PYRA5B_REF,PYRA5B_REF+15);
+ return GaussInfo(NORM_PYRA5,_my_gauss_coord,getNbGauss(),c,5);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for PYRA13 !");
+ }
+ case NORM_PENTA15:
+ {
+ std::vector<double> a(PENTA15A_REF,PENTA15A_REF+45),b(PENTA15B_REF,PENTA15B_REF+45);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(PENTA6A_REF,PENTA6A_REF+18);
+ return GaussInfo(NORM_PENTA6,_my_gauss_coord,getNbGauss(),c,6);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(PENTA6B_REF,PENTA6B_REF+18);
+ return GaussInfo(NORM_PENTA6,_my_gauss_coord,getNbGauss(),c,6);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for PENTA15 !");
+ }
+ case NORM_PENTA18:
+ {
+ std::vector<double> a(PENTA18A_REF,PENTA18A_REF+54),b(PENTA18B_REF,PENTA18B_REF+54);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(PENTA6A_REF,PENTA6A_REF+18);
+ return GaussInfo(NORM_PENTA6,_my_gauss_coord,getNbGauss(),c,6);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(PENTA6B_REF,PENTA6B_REF+18);
+ return GaussInfo(NORM_PENTA6,_my_gauss_coord,getNbGauss(),c,6);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for PENTA18 !");
+ }
+ case NORM_HEXA20:
+ {
+ std::vector<double> a(HEXA20A_REF,HEXA20A_REF+60),b(HEXA20B_REF,HEXA20B_REF+60);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(HEXA8A_REF,HEXA8A_REF+24);
+ return GaussInfo(NORM_HEXA8,_my_gauss_coord,getNbGauss(),c,8);
+ }
+ if(IsSatisfy(b,_my_reference_coord))
+ {
+ std::vector<double> c(HEXA8B_REF,HEXA8B_REF+24);
+ return GaussInfo(NORM_HEXA8,_my_gauss_coord,getNbGauss(),c,8);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for HEXA20 !");
+ }
+ case NORM_HEXA27:
+ {
+ std::vector<double> a(HEXA27A_REF,HEXA27A_REF+81);
+ if(IsSatisfy(a,_my_reference_coord))
+ {
+ std::vector<double> c(HEXA8B_REF,HEXA8B_REF+24);
+ return GaussInfo(NORM_HEXA8,_my_gauss_coord,getNbGauss(),c,8);
+ }
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not recognized pattern for HEXA27 !");
+ }
+ default:
+ throw INTERP_KERNEL::Exception("GaussInfo::convertToLinear : not implemented yet for other types than TETRA10, HEXA20, HEXA27, TRI3, QUAD8, QUAD8, PYRA13, PENTA15 !");
+ }
+}
+
+
+bool GaussInfo::IsSatisfy(const std::vector<double>& ref1, const std::vector<double>& ref2)
+{
+ std::size_t sz(ref1.size());
+ if(sz!=ref2.size())
+ return false;
+ for(std::size_t i=0;i<sz;i++)
+ if(!IsEqual(ref1[i],ref2[i]))
+ return false;
+ return true;
+}
+
/*!
* Check coordinates
*/
return ret;
}
+std::vector<double> GaussInfo::GetDefaultReferenceCoordinatesOf(NormalizedCellType ct)
+{
+ switch(ct)
+ {
+ case INTERP_KERNEL::NORM_SEG2:
+ return std::vector<double>(SEG2A_REF,SEG2A_REF+sizeof(SEG2A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_SEG3:
+ return std::vector<double>(SEG3_REF,SEG3_REF+sizeof(SEG3_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_TRI3:
+ return std::vector<double>(TRIA3A_REF,TRIA3A_REF+sizeof(TRIA3A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_TRI6:
+ return std::vector<double>(TRIA6A_REF,TRIA6A_REF+sizeof(TRIA6A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_TRI7:
+ return std::vector<double>(TRIA7A_REF,TRIA7A_REF+sizeof(TRIA7A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_QUAD4:
+ return std::vector<double>(QUAD4A_REF,QUAD4A_REF+sizeof(QUAD4A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_QUAD8:
+ return std::vector<double>(QUAD8A_REF,QUAD8A_REF+sizeof(QUAD8A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_QUAD9:
+ return std::vector<double>(QUAD9A_REF,QUAD9A_REF+sizeof(QUAD9A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_TETRA4:
+ return std::vector<double>(TETRA4A_REF,TETRA4A_REF+sizeof(TETRA4A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_TETRA10:
+ return std::vector<double>(TETRA10A_REF,TETRA10A_REF+sizeof(TETRA10A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_PYRA5:
+ return std::vector<double>(PYRA5A_REF,PYRA5A_REF+sizeof(PYRA5A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_PYRA13:
+ return std::vector<double>(PYRA13A_REF,PYRA13A_REF+sizeof(PYRA13A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_PENTA6:
+ return std::vector<double>(PENTA6A_REF,PENTA6A_REF+sizeof(PENTA6A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_PENTA15:
+ return std::vector<double>(PENTA15A_REF,PENTA15A_REF+sizeof(PENTA15A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_PENTA18:
+ return std::vector<double>(PENTA18A_REF,PENTA18A_REF+sizeof(PENTA18A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_HEXA8:
+ return std::vector<double>(HEXA8A_REF,HEXA8A_REF+sizeof(HEXA8A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_HEXA20:
+ return std::vector<double>(HEXA20A_REF,HEXA20A_REF+sizeof(HEXA20A_REF)/sizeof(double));
+ case INTERP_KERNEL::NORM_HEXA27:
+ return std::vector<double>(HEXA27A_REF,HEXA27A_REF+sizeof(HEXA27A_REF)/sizeof(double));
+ default:
+ THROW_IK_EXCEPTION("Input type " << ct << "is not managed by GetDefaultReferenceCoordinatesOf")
+ }
+}
+
+/*!
+ * Returns true if \a ptInRefCoo is in reference cell of type \a ct (regarding GetDefaultReferenceCoordinatesOf)
+ * \sa GetDefaultReferenceCoordinatesOf
+ */
+bool GaussInfo::IsInOrOutForReference(NormalizedCellType ct, const double *ptInRefCoo, double eps)
+{
+ switch(ct)
+ {
+ case INTERP_KERNEL::NORM_SEG2:
+ case INTERP_KERNEL::NORM_SEG3:
+ {
+ return std::fabs(ptInRefCoo[0]) < 1.0+eps;
+ }
+ case INTERP_KERNEL::NORM_QUAD4:
+ case INTERP_KERNEL::NORM_QUAD8:
+ case INTERP_KERNEL::NORM_QUAD9:
+ {
+ return std::find_if(ptInRefCoo,ptInRefCoo+2,[eps](double v){ return std::fabs(v) > 1.0+eps; }) == ptInRefCoo+2;
+ }
+ case INTERP_KERNEL::NORM_HEXA8:
+ case INTERP_KERNEL::NORM_HEXA20:
+ case INTERP_KERNEL::NORM_HEXA27:
+ {
+ return std::find_if(ptInRefCoo,ptInRefCoo+3,[eps](double v){ return std::fabs(v) > 1.0+eps; }) == ptInRefCoo+3;
+ }
+ default:
+ THROW_IK_EXCEPTION("IsInOrOutForReference : not implemented for this geo type !");
+ }
+}
+
typedef void (*MapToShapeFunction)(GaussInfo& obj);
/*!
case NORM_SEG2:
_my_local_ref_dim = 1;
_my_local_nb_ref = 2;
- seg2Init();
+ seg2aInit();
aSatify = isSatisfy();
- CHECK_MACRO;
+ if(!aSatify)
+ {
+ seg2bInit();
+ aSatify = isSatisfy();
+ CHECK_MACRO;
+ }
break;
case NORM_SEG3:
break;
}
-
- _my_local_ref_dim = 3;
+/* _my_local_ref_dim = 3;
_my_local_nb_ref = 6;
penta6aInit();
aSatify = isSatisfy();
CHECK_MACRO;
}
break;
-
+*/
case NORM_PENTA15:
{
_my_local_ref_dim = 3;
break;
}
+ case NORM_PENTA18:
+ {
+ _my_local_ref_dim = 3;
+ _my_local_nb_ref = 18;
+ MapToShapeFunction PENTA18PTR[]={Penta18aInit,Penta18bInit};
+ std::size_t NB_OF_PENTA18PTR(sizeof(PENTA18PTR)/sizeof(MapToShapeFunction));
+ for(std::size_t i=0;i<NB_OF_PENTA18PTR && !aSatify;i++)
+ {
+ (PENTA18PTR[i])(*this);
+ aSatify = isSatisfy();
+ }
+ CHECK_MACRO;
+ break;
+ }
+
case NORM_HEXA8:
{
_my_local_ref_dim = 3;
/**
* Return shape function value by node id
*/
-const double* GaussInfo::getFunctionValues( const int theGaussId ) const
+const double *GaussInfo::getFunctionValues( const int theGaussId ) const
+{
+ return _my_function_value.data() + _my_nb_ref*theGaussId ;
+}
+
+/**
+ * Return the derivative of shape function value by node id
+ */
+const double *GaussInfo::getDerivativeOfShapeFunctionAt( const int theGaussId ) const
{
- return &_my_function_value[ _my_nb_ref*theGaussId ];
+ return _my_derivative_func_value.data() + _my_nb_ref*getReferenceCoordDim()*theGaussId;
}
void GaussInfo::point1Init()
/*!
* Init Segment 2 Reference coordinates ans Shape function.
*/
-void GaussInfo::seg2Init()
+void GaussInfo::seg2aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
+ case 0:
+ coords[0] = SEG2A_REF[0];
break;
- case 1:
- coords[0] = 1.0;
+ case 1:
+ coords[0] = SEG2A_REF[1];
break;
- LOCAL_COORD_MACRO_END;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*(1.0 - gc[0]);
+ funValue[1] = 0.5*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+}
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*(1.0 - gc[0]);
- funValue[1] = 0.5*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
+void GaussInfo::seg2bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = SEG2B_REF[0];
+ break;
+ case 1:
+ coords[0] = SEG2B_REF[1];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 1.0 - gc[0];
+ funValue[1] = gc[0];
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::seg3Init()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
+ case 0:
+ coords[0] = SEG3_REF[0];
break;
- case 1:
- coords[0] = 1.0;
+ case 1:
+ coords[0] = SEG3_REF[1];
break;
- case 2:
- coords[0] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = -0.5*(1.0 - gc[0])*gc[0];
- funValue[1] = 0.5*(1.0 + gc[0])*gc[0];
- funValue[2] = (1.0 + gc[0])*(1.0 - gc[0]);
- SHAPE_FUN_MACRO_END;
+ case 2:
+ coords[0] = SEG3_REF[2];
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = -0.5*(1.0 - gc[0])*gc[0];
+ funValue[1] = 0.5*(1.0 + gc[0])*gc[0];
+ funValue[2] = (1.0 + gc[0])*(1.0 - gc[0]);
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tria3aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = TRIA3A_REF[0];
+ coords[1] = TRIA3A_REF[1];
+ break;
+ case 1:
+ coords[0] = TRIA3A_REF[2];
+ coords[1] = TRIA3A_REF[3];
+ break;
+ case 2:
+ coords[0] = TRIA3A_REF[4];
+ coords[1] = TRIA3A_REF[5];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*(1.0 + gc[1]);
+ funValue[1] = -0.5*(gc[0] + gc[1]);
+ funValue[2] = 0.5*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*(1.0 + gc[1]);
- funValue[1] = -0.5*(gc[0] + gc[1]);
- funValue[2] = 0.5*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[0] = 0.0 ;
+ devFunValue[1] = 0.5 ;
+
+ devFunValue[2] = -0.5;
+ devFunValue[3] = -0.5;
+
+ devFunValue[4] = 0.5;
+ devFunValue[5] = 0.0;
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tria3bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 0.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = TRIA3B_REF[0];
+ coords[1] = TRIA3B_REF[1];
+ break;
+ case 1:
+ coords[0] = TRIA3B_REF[2];
+ coords[1] = TRIA3B_REF[3];
+ break;
+ case 2:
+ coords[0] = TRIA3B_REF[4];
+ coords[1] = TRIA3B_REF[5];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 1.0 - gc[0] - gc[1];
+ funValue[1] = gc[0];
+ funValue[2] = gc[1];
+ SHAPE_FUN_MACRO_END;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 1.0 - gc[0] - gc[1];
- funValue[1] = gc[0];
- funValue[2] = gc[1];
- SHAPE_FUN_MACRO_END;
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = -1.0 ;
+ devFunValue[1] = -1.0 ;
+
+ devFunValue[2] = 1.0 ;
+ devFunValue[3] = 0.0 ;
+
+ devFunValue[4] = 0.0 ;
+ devFunValue[5] = 1.0 ;
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tria6aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = -1.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = TRIA6A_REF[0];
+ coords[1] = TRIA6A_REF[1];
+ break;
+ case 1:
+ coords[0] = TRIA6A_REF[2];
+ coords[1] = TRIA6A_REF[3];
+ break;
+ case 2:
+ coords[0] = TRIA6A_REF[4];
+ coords[1] = TRIA6A_REF[5];
+ break;
+ case 3:
+ coords[0] = TRIA6A_REF[6];
+ coords[1] = TRIA6A_REF[7];
+ break;
+ case 4:
+ coords[0] = TRIA6A_REF[8];
+ coords[1] = TRIA6A_REF[9];
+ break;
+ case 5:
+ coords[0] = TRIA6A_REF[10];
+ coords[1] = TRIA6A_REF[11];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*(1.0 + gc[1])*gc[1];
+ funValue[1] = 0.5*(gc[0] + gc[1])*(gc[0] + gc[1] + 1);
+ funValue[2] = 0.5*(1.0 + gc[0])*gc[0];
+ funValue[3] = -1.0*(1.0 + gc[1])*(gc[0] + gc[1]);
+ funValue[4] = -1.0*(1.0 + gc[0])*(gc[0] + gc[1]);
+ funValue[5] = (1.0 + gc[1])*(1.0 + gc[1]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*(1.0 + gc[1])*gc[1];
- funValue[1] = 0.5*(gc[0] + gc[1])*(gc[0] + gc[1] + 1);
- funValue[2] = 0.5*(1.0 + gc[0])*gc[0];
- funValue[3] = -1.0*(1.0 + gc[1])*(gc[0] + gc[1]);
- funValue[4] = -1.0*(1.0 + gc[0])*(gc[0] + gc[1]);
- funValue[5] = (1.0 + gc[1])*(1.0 + gc[1]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[0] = 0.0;
+ devFunValue[1] = 0.5*( 2*gc[1] + 1.0 );
+
+ devFunValue[2] = 0.5*( 2*gc[0] + 2.0*gc[1] + 1.0);
+ devFunValue[3] = 0.5*( 2*gc[1] + 2.0*gc[0] + 1.0);
+
+ devFunValue[4] = gc[0] + 0.5;
+ devFunValue[5] = 0.0;
+
+ devFunValue[6] = -1.0*(1.0 + gc[1]);
+ devFunValue[7] = -1.0*(2*gc[1]+gc[0]+1.0);
+
+ devFunValue[8] = -1.0*(2*gc[0]+gc[1]+1.0);
+ devFunValue[9] = -1.0*(1.0 + gc[0]);
+
+ devFunValue[10] = 0.0;
+ devFunValue[11] = (2*gc[1]+2.0);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tria6bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 0.0;
- break;
+ case 0:
+ coords[0] = TRIA6B_REF[0];
+ coords[1] = TRIA6B_REF[1];
+ break;
+ case 1:
+ coords[0] = TRIA6B_REF[2];
+ coords[1] = TRIA6B_REF[3];
+ break;
+ case 2:
+ coords[0] = TRIA6B_REF[4];
+ coords[1] = TRIA6B_REF[5];
+ break;
+ case 3:
+ coords[0] = TRIA6B_REF[6];
+ coords[1] = TRIA6B_REF[7];
+ break;
+ case 4:
+ coords[0] = TRIA6B_REF[8];
+ coords[1] = TRIA6B_REF[9];
+ break;
+ case 5:
+ coords[0] = TRIA6B_REF[10];
+ coords[1] = TRIA6B_REF[11];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = (1.0 - gc[0] - gc[1])*(1.0 - 2.0*gc[0] - 2.0*gc[1]);
+ funValue[1] = gc[0]*(2.0*gc[0] - 1.0);
+ funValue[2] = gc[1]*(2.0*gc[1] - 1.0);
+ funValue[3] = 4.0*gc[0]*(1.0 - gc[0] - gc[1]);
+ funValue[4] = 4.0*gc[0]*gc[1];
+ funValue[5] = 4.0*gc[1]*(1.0 - gc[0] - gc[1]);
+ SHAPE_FUN_MACRO_END;
- case 1:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- case 2:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
+ devFunValue[0] = 4*gc[0] + 4*gc[1] - 3.0 ;
+ devFunValue[1] = 4*gc[1] + 4*gc[0] - 3.0 ;
- case 3:
- coords[0] = 0.5;
- coords[1] = 0.0;
- break;
+ devFunValue[2] = 4*gc[0] - 1.0 ;
+ devFunValue[3] = 0.0 ;
- case 4:
- coords[0] = 0.5;
- coords[1] = 0.5;
- break;
+ devFunValue[4] = 0.0 ;
+ devFunValue[5] = 4*gc[1] - 1.0 ;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.5;
- break;
+ devFunValue[6] = -8.0*gc[0] - 4.0 * gc[1] + 4.0 ;
+ devFunValue[7] = -4.0*gc[0] ;
- LOCAL_COORD_MACRO_END;
+ devFunValue[8] = 4.0*gc[1];
+ devFunValue[9] = 4.0*gc[0];
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = (1.0 - gc[0] - gc[1])*(1.0 - 2.0*gc[0] - 2.0*gc[1]);
- funValue[1] = gc[0]*(2.0*gc[0] - 1.0);
- funValue[2] = gc[1]*(2.0*gc[1] - 1.0);
- funValue[3] = 4.0*gc[0]*(1.0 - gc[0] - gc[1]);
- funValue[4] = 4.0*gc[0]*gc[1];
- funValue[5] = 4.0*gc[1]*(1.0 - gc[0] - gc[1]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[10] = -4.0*gc[1] ;
+ devFunValue[11] = -8.0*gc[1] - 4.0*gc[0] + 4.0 ;
+
+ DEV_SHAPE_FUN_MACRO_END;
}
void GaussInfo::tria7aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 0.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
- case 3:
- coords[0] = 0.5;
- coords[1] = 0.0;
- break;
- case 4:
- coords[0] = 0.5;
- coords[1] = 0.5;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.5;
- break;
- case 6:
- coords[0] = 0.3333333333333333;
- coords[1] = 0.3333333333333333;
- break;
-
+ case 0:
+ coords[0] = TRIA7A_REF[0];
+ coords[1] = TRIA7A_REF[1];
+ break;
+ case 1:
+ coords[0] = TRIA7A_REF[2];
+ coords[1] = TRIA7A_REF[3];
+ break;
+ case 2:
+ coords[0] = TRIA7A_REF[4];
+ coords[1] = TRIA7A_REF[5];
+ break;
+ case 3:
+ coords[0] = TRIA7A_REF[6];
+ coords[1] = TRIA7A_REF[7];
+ break;
+ case 4:
+ coords[0] = TRIA7A_REF[8];
+ coords[1] = TRIA7A_REF[9];
+ break;
+ case 5:
+ coords[0] = TRIA7A_REF[10];
+ coords[1] = TRIA7A_REF[11];
+ break;
+ case 6:
+ coords[0] = TRIA7A_REF[12];
+ coords[1] = TRIA7A_REF[13];
+ break;
LOCAL_COORD_MACRO_END;
SHAPE_FUN_MACRO_BEGIN;
void GaussInfo::quad4aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- break;
+ case 0:
+ coords[0] = QUAD4A_REF[0];
+ coords[1] = QUAD4A_REF[1];
+ break;
+ case 1:
+ coords[0] = QUAD4A_REF[2];
+ coords[1] = QUAD4A_REF[3];
+ break;
+ case 2:
+ coords[0] = QUAD4A_REF[4];
+ coords[1] = QUAD4A_REF[5];
+ break;
+ case 3:
+ coords[0] = QUAD4A_REF[6];
+ coords[1] = QUAD4A_REF[7];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0]);
+ funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0]);
+ funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0]);
+ funValue[3] = 0.25*(1.0 + gc[1])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- LOCAL_COORD_MACRO_END;
+ devFunValue[0] = -0.25*(1.0 + gc[1]);
+ devFunValue[1] = 0.25*(1.0 - gc[0]);
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0]);
- funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0]);
- funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0]);
- funValue[3] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[2] = -0.25*(1.0 - gc[1]);
+ devFunValue[3] = -0.25*(1.0 - gc[0]);
+
+ devFunValue[4] = 0.25*(1.0 - gc[1]);
+ devFunValue[5] = -0.25*(1.0 + gc[0]);
+
+ devFunValue[6] = 0.25*(1.0 + gc[1]);
+ devFunValue[7] = 0.25*(1.0 + gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::quad4bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
+ case 0:
+ coords[0] = QUAD4B_REF[0];
+ coords[1] = QUAD4B_REF[1];
+ break;
+ case 1:
+ coords[0] = QUAD4B_REF[2];
+ coords[1] = QUAD4B_REF[3];
+ break;
+ case 2:
+ coords[0] = QUAD4B_REF[4];
+ coords[1] = QUAD4B_REF[5];
+ break;
+ case 3:
+ coords[0] = QUAD4B_REF[6];
+ coords[1] = QUAD4B_REF[7];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
+ funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
+ funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- LOCAL_COORD_MACRO_END;
+ devFunValue[0] = -0.25*(1.0 - gc[1]);
+ devFunValue[1] = -0.25*(1.0 - gc[0]);
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
- funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
- funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[2] = 0.25*(1.0 - gc[1]);
+ devFunValue[3] = -0.25*(1.0 + gc[0]);
+
+ devFunValue[4] = 0.25*(1.0 + gc[1]);
+ devFunValue[5] = 0.25*(1.0 + gc[0]);
+
+ devFunValue[6] = -0.25*(1.0 + gc[1]);
+ devFunValue[7] = 0.25*(1.0 - gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
void GaussInfo::quad4cInit()
funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
funValue[3] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = -0.25*(1.0 - gc[1]);
+ devFunValue[1] = -0.25*(1.0 - gc[0]);
+
+ devFunValue[2] = -0.25*(1.0 + gc[1]);
+ devFunValue[3] = 0.25*(1.0 - gc[0]);
+
+ devFunValue[4] = 0.25*(1.0 + gc[0]);
+ devFunValue[5] = 0.25*(1.0 + gc[1]);
+
+ devFunValue[6] = 0.25*(1.0 - gc[1]);
+ devFunValue[7] = -0.25*(1.0 + gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
- * This shapefunc map is same as degenerated seg2Init
+ * This shapefunc map is same as degenerated seg2aInit
*/
void GaussInfo::quad4DegSeg2Init()
{
void GaussInfo::quad8aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = 0.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = -1.0;
- break;
- case 6:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
- case 7:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = QUAD8A_REF[0];
+ coords[1] = QUAD8A_REF[1];
+ break;
+ case 1:
+ coords[0] = QUAD8A_REF[2];
+ coords[1] = QUAD8A_REF[3];
+ break;
+ case 2:
+ coords[0] = QUAD8A_REF[4];
+ coords[1] = QUAD8A_REF[5];
+ break;
+ case 3:
+ coords[0] = QUAD8A_REF[6];
+ coords[1] = QUAD8A_REF[7];
+ break;
+ case 4:
+ coords[0] = QUAD8A_REF[8];
+ coords[1] = QUAD8A_REF[9];
+ break;
+ case 5:
+ coords[0] = QUAD8A_REF[10];
+ coords[1] = QUAD8A_REF[11];
+ break;
+ case 6:
+ coords[0] = QUAD8A_REF[12];
+ coords[1] = QUAD8A_REF[13];
+ break;
+ case 7:
+ coords[0] = QUAD8A_REF[14];
+ coords[1] = QUAD8A_REF[15];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0])*(gc[1] - gc[0] - 1.0);
+ funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0])*(-gc[1] - gc[0] - 1.0);
+ funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0])*(-gc[1] + gc[0] - 1.0);
+ funValue[3] = 0.25*(1.0 + gc[1])*(1.0 + gc[0])*(gc[1] + gc[0] - 1.0);
+ funValue[4] = 0.5*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[1]);
+ funValue[5] = 0.5*(1.0 - gc[1])*(1.0 - gc[0])*(1.0 + gc[0]);
+ funValue[6] = 0.5*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[1]);
+ funValue[7] = 0.5*(1.0 + gc[1])*(1.0 - gc[0])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0])*(gc[1] - gc[0] - 1.0);
- funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0])*(-gc[1] - gc[0] - 1.0);
- funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0])*(-gc[1] + gc[0] - 1.0);
- funValue[3] = 0.25*(1.0 + gc[1])*(1.0 + gc[0])*(gc[1] + gc[0] - 1.0);
- funValue[4] = 0.5*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[1]);
- funValue[5] = 0.5*(1.0 - gc[1])*(1.0 - gc[0])*(1.0 + gc[0]);
- funValue[6] = 0.5*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[1]);
- funValue[7] = 0.5*(1.0 + gc[1])*(1.0 - gc[0])*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.25*(1.0 + gc[1])*(2*gc[0]-gc[1]);
+ devFunValue[1] = 0.25*(1.0 - gc[0])*(2*gc[1]-gc[0]);
+
+ devFunValue[2] = 0.25*(1.0 - gc[1])*(2*gc[0]+gc[1]);
+ devFunValue[3] = 0.25*(1.0 - gc[0])*(2*gc[1]+gc[0]);
+
+ devFunValue[4] = 0.25*(1.0 - gc[1])*(2*gc[0]-gc[1]);
+ devFunValue[5] = 0.25*(1.0 + gc[0])*(2*gc[1]-gc[0]);
+
+ devFunValue[6] = 0.25*(1.0 + gc[1])*(2*gc[0]+gc[1]);
+ devFunValue[7] = 0.25*(1.0 + gc[0])*(2*gc[1]+gc[0]);
+
+ devFunValue[8] = -0.5*(1.0 - gc[1])*(1.0 + gc[1]);
+ devFunValue[9] = 0.5*(1.0 - gc[0])*(-2*gc[1]);
+
+ devFunValue[10] = 0.5*(1.0 - gc[1])*(-2*gc[0]);
+ devFunValue[11] = -0.5*(1.0 - gc[0])*(1.0 + gc[0]);
+
+ devFunValue[12] = 0.5*(1.0 - gc[1])*(1.0 + gc[1]);
+ devFunValue[13] = 0.5*(1.0 + gc[0])*(-2*gc[1]);
+
+ devFunValue[14] = 0.5*(1.0 + gc[1])*(-2*gc[0]);
+ devFunValue[15] = 0.5*(1.0 - gc[0])*(1.0 + gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::quad8bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = -1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1])*(-1.0 - gc[0] - gc[1]);
- funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1])*(-1.0 + gc[0] - gc[1]);
- funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1])*(-1.0 + gc[0] + gc[1]);
- funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1])*(-1.0 - gc[0] + gc[1]);
- funValue[4] = 0.5*(1.0 - gc[0]*gc[0])*(1.0 - gc[1]);
- funValue[5] = 0.5*(1.0 - gc[1]*gc[1])*(1.0 + gc[0]);
- funValue[6] = 0.5*(1.0 - gc[0]*gc[0])*(1.0 + gc[1]);
- funValue[7] = 0.5*(1.0 - gc[1]*gc[1])*(1.0 - gc[0]);
- SHAPE_FUN_MACRO_END;
-}
-
-void GaussInfo::quad9aInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = -1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 1.0;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 0.0;
- break;
- case 8:
- coords[0] = 0.0;
- coords[1] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = QUAD8B_REF[0];
+ coords[1] = QUAD8B_REF[1];
+ break;
+ case 1:
+ coords[0] = QUAD8B_REF[2];
+ coords[1] = QUAD8B_REF[3];
+ break;
+ case 2:
+ coords[0] = QUAD8B_REF[4];
+ coords[1] = QUAD8B_REF[5];
+ break;
+ case 3:
+ coords[0] = QUAD8B_REF[6];
+ coords[1] = QUAD8B_REF[7];
+ break;
+ case 4:
+ coords[0] = QUAD8B_REF[8];
+ coords[1] = QUAD8B_REF[9];
+ break;
+ case 5:
+ coords[0] = QUAD8B_REF[10];
+ coords[1] = QUAD8B_REF[11];
+ break;
+ case 6:
+ coords[0] = QUAD8B_REF[12];
+ coords[1] = QUAD8B_REF[13];
+ break;
+ case 7:
+ coords[0] = QUAD8B_REF[14];
+ coords[1] = QUAD8B_REF[15];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1])*(-1.0 - gc[0] - gc[1]);
+ funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1])*(-1.0 + gc[0] - gc[1]);
+ funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1])*(-1.0 + gc[0] + gc[1]);
+ funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1])*(-1.0 - gc[0] + gc[1]);
+ funValue[4] = 0.5*(1.0 - gc[0]*gc[0])*(1.0 - gc[1]);
+ funValue[5] = 0.5*(1.0 - gc[1]*gc[1])*(1.0 + gc[0]);
+ funValue[6] = 0.5*(1.0 - gc[0]*gc[0])*(1.0 + gc[1]);
+ funValue[7] = 0.5*(1.0 - gc[1]*gc[1])*(1.0 - gc[0]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*gc[0]*gc[1]*(gc[0]-1.)*(gc[1]-1.);
- funValue[1] = 0.25*gc[0]*gc[1]*(gc[0]+1.)*(gc[1]-1.);
- funValue[2] = 0.25*gc[0]*gc[1]*(gc[0]+1.)*(gc[1]+1.);
- funValue[3] = 0.25*gc[0]*gc[1]*(gc[0]-1.)*(gc[1]+1.);
- funValue[4] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.);
- funValue[5] = 0.5*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1]);
- funValue[6] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.);
- funValue[7] = 0.5*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1]);
- funValue[8] = (1.-gc[0]*gc[0])*(1.-gc[1]*gc[1]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[0] = 0.25*(1.0 - gc[1])*(2*gc[0] + gc[1] );
+ devFunValue[1] = 0.25*(1.0 - gc[0])*(2*gc[1] + gc[0] );
+
+ devFunValue[2] = 0.25*(1.0 - gc[1])*(2*gc[0] - gc[1] );
+ devFunValue[3] = 0.25*(1.0 + gc[0])*(2*gc[1] - gc[0] );
+
+ devFunValue[4] = 0.25*(1.0 + gc[1])*(2*gc[0] + gc[1]);
+ devFunValue[5] = 0.25*(1.0 + gc[0])*(2*gc[1] + gc[0]);
+
+ devFunValue[6] = 0.25*( 1.0 + gc[1])*(2*gc[0] - gc[1] );
+ devFunValue[7] = 0.25*(1.0 - gc[0])*(2*gc[1] - gc[0]);
+
+ devFunValue[8] = -(1.0 - gc[1])*gc[0];
+ devFunValue[9] = -0.5*(1.0 - gc[0]*gc[0]);
+
+ devFunValue[10] = 0.5*(1.0 - gc[1]*gc[1]);
+ devFunValue[11] = -(1.0 + gc[0])*gc[1];
+
+ devFunValue[12] = -(1.0 + gc[1])*gc[0];
+ devFunValue[13] = 0.5*(1.0 - gc[0]*gc[0]);
+
+ devFunValue[14] = -0.5*(1.0 - gc[1]*gc[1]);
+ devFunValue[15] = -(1.0 - gc[0])*gc[1];
+
+ DEV_SHAPE_FUN_MACRO_END;
+}
+
+void GaussInfo::quad9aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = QUAD9A_REF[0];
+ coords[1] = QUAD9A_REF[1];
+ break;
+ case 1:
+ coords[0] = QUAD9A_REF[2];
+ coords[1] = QUAD9A_REF[3];
+ break;
+ case 2:
+ coords[0] = QUAD9A_REF[4];
+ coords[1] = QUAD9A_REF[5];
+ break;
+ case 3:
+ coords[0] = QUAD9A_REF[6];
+ coords[1] = QUAD9A_REF[7];
+ break;
+ case 4:
+ coords[0] = QUAD9A_REF[8];
+ coords[1] = QUAD9A_REF[9];
+ break;
+ case 5:
+ coords[0] = QUAD9A_REF[10];
+ coords[1] = QUAD9A_REF[11];
+ break;
+ case 6:
+ coords[0] = QUAD9A_REF[12];
+ coords[1] = QUAD9A_REF[13];
+ break;
+ case 7:
+ coords[0] = QUAD9A_REF[14];
+ coords[1] = QUAD9A_REF[15];
+ break;
+ case 8:
+ coords[0] = QUAD9A_REF[16];
+ coords[1] = QUAD9A_REF[17];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*gc[0]*gc[1]*(gc[0]-1.)*(gc[1]-1.);
+ funValue[1] = 0.25*gc[0]*gc[1]*(gc[0]+1.)*(gc[1]-1.);
+ funValue[2] = 0.25*gc[0]*gc[1]*(gc[0]+1.)*(gc[1]+1.);
+ funValue[3] = 0.25*gc[0]*gc[1]*(gc[0]-1.)*(gc[1]+1.);
+ funValue[4] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.);
+ funValue[5] = 0.5*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1]);
+ funValue[6] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.);
+ funValue[7] = 0.5*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1]);
+ funValue[8] = (1.-gc[0]*gc[0])*(1.-gc[1]*gc[1]);
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tetra4aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = gc[1];
- funValue[1] = gc[2];
- funValue[2] = 1.0 - gc[0] - gc[1] - gc[2];
- funValue[3] = gc[0];
- SHAPE_FUN_MACRO_END;
+ case 0:
+ coords[0] = TETRA4A_REF[0];
+ coords[1] = TETRA4A_REF[1];
+ coords[2] = TETRA4A_REF[2];
+ break;
+ case 1:
+ coords[0] = TETRA4A_REF[3];
+ coords[1] = TETRA4A_REF[4];
+ coords[2] = TETRA4A_REF[5];
+ break;
+ case 2:
+ coords[0] = TETRA4A_REF[6];
+ coords[1] = TETRA4A_REF[7];
+ coords[2] = TETRA4A_REF[8];
+ break;
+ case 3:
+ coords[0] = TETRA4A_REF[9];
+ coords[1] = TETRA4A_REF[10];
+ coords[2] = TETRA4A_REF[11];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = gc[1];
+ funValue[1] = gc[2];
+ funValue[2] = 1.0 - gc[0] - gc[1] - gc[2];
+ funValue[3] = gc[0];
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tetra4bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = gc[1];
- funValue[2] = gc[2];
- funValue[1] = 1.0 - gc[0] - gc[1] - gc[2];
- funValue[3] = gc[0];
- SHAPE_FUN_MACRO_END;
-
+ case 0:
+ coords[0] = TETRA4B_REF[0];
+ coords[1] = TETRA4B_REF[1];
+ coords[2] = TETRA4B_REF[2];
+ break;
+ case 1:
+ coords[0] = TETRA4B_REF[3];
+ coords[1] = TETRA4B_REF[4];
+ coords[2] = TETRA4B_REF[5];
+ break;
+ case 2:
+ coords[0] = TETRA4B_REF[6];
+ coords[1] = TETRA4B_REF[7];
+ coords[2] = TETRA4B_REF[8];
+ break;
+ case 3:
+ coords[0] = TETRA4B_REF[9];
+ coords[1] = TETRA4B_REF[10];
+ coords[2] = TETRA4B_REF[11];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = gc[1];
+ funValue[2] = gc[2];
+ funValue[1] = 1.0 - gc[0] - gc[1] - gc[2];
+ funValue[3] = gc[0];
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tetra10aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = 0.5;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 8:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 9:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = TETRA10A_REF[0];
+ coords[1] = TETRA10A_REF[1];
+ coords[2] = TETRA10A_REF[2];
+ break;
+ case 1:
+ coords[0] = TETRA10A_REF[3];
+ coords[1] = TETRA10A_REF[4];
+ coords[2] = TETRA10A_REF[5];
+ break;
+ case 2:
+ coords[0] = TETRA10A_REF[6];
+ coords[1] = TETRA10A_REF[7];
+ coords[2] = TETRA10A_REF[8];
+ break;
+ case 3:
+ coords[0] = TETRA10A_REF[9];
+ coords[1] = TETRA10A_REF[10];
+ coords[2] = TETRA10A_REF[11];
+ break;
+ case 4:
+ coords[0] = TETRA10A_REF[12];
+ coords[1] = TETRA10A_REF[13];
+ coords[2] = TETRA10A_REF[14];
+ break;
+ case 5:
+ coords[0] = TETRA10A_REF[15];
+ coords[1] = TETRA10A_REF[16];
+ coords[2] = TETRA10A_REF[17];
+ break;
+ case 6:
+ coords[0] = TETRA10A_REF[18];
+ coords[1] = TETRA10A_REF[19];
+ coords[2] = TETRA10A_REF[20];
+ break;
+ case 7:
+ coords[0] = TETRA10A_REF[21];
+ coords[1] = TETRA10A_REF[22];
+ coords[2] = TETRA10A_REF[23];
+ break;
+ case 8:
+ coords[0] = TETRA10A_REF[24];
+ coords[1] = TETRA10A_REF[25];
+ coords[2] = TETRA10A_REF[26];
+ break;
+ case 9:
+ coords[0] = TETRA10A_REF[27];
+ coords[1] = TETRA10A_REF[28];
+ coords[2] = TETRA10A_REF[29];
+ break;
+ LOCAL_COORD_MACRO_END;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = gc[1]*(2.0*gc[1] - 1.0);
- funValue[1] = gc[2]*(2.0*gc[2] - 1.0);
- funValue[2] = (1.0 - gc[0] - gc[1] - gc[2])*(1.0 - 2.0*gc[0] - 2.0*gc[1] - 2.0*gc[2]);
- funValue[3] = gc[0]*(2.0*gc[0] - 1.0);
- funValue[4] = 4.0*gc[1]*gc[2];
- funValue[5] = 4.0*gc[2]*(1.0 - gc[0] - gc[1] - gc[2]);
- funValue[6] = 4.0*gc[1]*(1.0 - gc[0] - gc[1] - gc[2]);
- funValue[7] = 4.0*gc[0]*gc[1];
- funValue[8] = 4.0*gc[0]*gc[2];
- funValue[9] = 4.0*gc[0]*(1.0 - gc[0] - gc[1] - gc[2]);
- SHAPE_FUN_MACRO_END;
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = gc[1]*(2.0*gc[1] - 1.0);
+ funValue[1] = gc[2]*(2.0*gc[2] - 1.0);
+ funValue[2] = (1.0 - gc[0] - gc[1] - gc[2])*(1.0 - 2.0*gc[0] - 2.0*gc[1] - 2.0*gc[2]);
+ funValue[3] = gc[0]*(2.0*gc[0] - 1.0);
+ funValue[4] = 4.0*gc[1]*gc[2];
+ funValue[5] = 4.0*gc[2]*(1.0 - gc[0] - gc[1] - gc[2]);
+ funValue[6] = 4.0*gc[1]*(1.0 - gc[0] - gc[1] - gc[2]);
+ funValue[7] = 4.0*gc[0]*gc[1];
+ funValue[8] = 4.0*gc[0]*gc[2];
+ funValue[9] = 4.0*gc[0]*(1.0 - gc[0] - gc[1] - gc[2]);
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::tetra10bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = 0.5;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 9:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 8:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = gc[1]*(2.0*gc[1] - 1.0);
- funValue[2] = gc[2]*(2.0*gc[2] - 1.0);
- funValue[1] = (1.0 - gc[0] - gc[1] - gc[2])*(1.0 - 2.0*gc[0] - 2.0*gc[1] - 2.0*gc[2]);
- funValue[3] = gc[0]*(2.0*gc[0] - 1.0);
- funValue[6] = 4.0*gc[1]*gc[2];
- funValue[5] = 4.0*gc[2]*(1.0 - gc[0] - gc[1] - gc[2]);
- funValue[4] = 4.0*gc[1]*(1.0 - gc[0] - gc[1] - gc[2]);
- funValue[7] = 4.0*gc[0]*gc[1];
- funValue[9] = 4.0*gc[0]*gc[2];
- funValue[8] = 4.0*gc[0]*(1.0 - gc[0] - gc[1] - gc[2]);
- SHAPE_FUN_MACRO_END;
+ case 0:
+ coords[0] = TETRA10B_REF[0];
+ coords[1] = TETRA10B_REF[1];
+ coords[2] = TETRA10B_REF[2];
+ break;
+ case 1:
+ coords[0] = TETRA10B_REF[3];
+ coords[1] = TETRA10B_REF[4];
+ coords[2] = TETRA10B_REF[5];
+ break;
+ case 2:
+ coords[0] = TETRA10B_REF[6];
+ coords[1] = TETRA10B_REF[7];
+ coords[2] = TETRA10B_REF[8];
+ break;
+ case 3:
+ coords[0] = TETRA10B_REF[9];
+ coords[1] = TETRA10B_REF[10];
+ coords[2] = TETRA10B_REF[11];
+ break;
+ case 4:
+ coords[0] = TETRA10B_REF[12];
+ coords[1] = TETRA10B_REF[13];
+ coords[2] = TETRA10B_REF[14];
+ break;
+ case 5:
+ coords[0] = TETRA10B_REF[15];
+ coords[1] = TETRA10B_REF[16];
+ coords[2] = TETRA10B_REF[17];
+ break;
+ case 6:
+ coords[0] = TETRA10B_REF[18];
+ coords[1] = TETRA10B_REF[19];
+ coords[2] = TETRA10B_REF[20];
+ break;
+ case 7:
+ coords[0] = TETRA10B_REF[21];
+ coords[1] = TETRA10B_REF[22];
+ coords[2] = TETRA10B_REF[23];
+ break;
+ case 8:
+ coords[0] = TETRA10B_REF[24];
+ coords[1] = TETRA10B_REF[25];
+ coords[2] = TETRA10B_REF[26];
+ break;
+ case 9:
+ coords[0] = TETRA10B_REF[27];
+ coords[1] = TETRA10B_REF[28];
+ coords[2] = TETRA10B_REF[29];
+ break;
+ LOCAL_COORD_MACRO_END;
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = gc[1]*(2.0*gc[1] - 1.0);
+ funValue[2] = gc[2]*(2.0*gc[2] - 1.0);
+ funValue[1] = (1.0 - gc[0] - gc[1] - gc[2])*(1.0 - 2.0*gc[0] - 2.0*gc[1] - 2.0*gc[2]);
+ funValue[3] = gc[0]*(2.0*gc[0] - 1.0);
+ funValue[6] = 4.0*gc[1]*gc[2];
+ funValue[5] = 4.0*gc[2]*(1.0 - gc[0] - gc[1] - gc[2]);
+ funValue[4] = 4.0*gc[1]*(1.0 - gc[0] - gc[1] - gc[2]);
+ funValue[7] = 4.0*gc[0]*gc[1];
+ funValue[9] = 4.0*gc[0]*gc[2];
+ funValue[8] = 4.0*gc[0]*(1.0 - gc[0] - gc[1] - gc[2]);
+ SHAPE_FUN_MACRO_END;
}
/*!
void GaussInfo::pyra5aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = PYRA5A_REF[0];
+ coords[1] = PYRA5A_REF[1];
+ coords[2] = PYRA5A_REF[2];
+ break;
+ case 1:
+ coords[0] = PYRA5A_REF[3];
+ coords[1] = PYRA5A_REF[4];
+ coords[2] = PYRA5A_REF[5];
+ break;
+ case 2:
+ coords[0] = PYRA5A_REF[6];
+ coords[1] = PYRA5A_REF[7];
+ coords[2] = PYRA5A_REF[8];
+ break;
+ case 3:
+ coords[0] = PYRA5A_REF[9];
+ coords[1] = PYRA5A_REF[10];
+ coords[2] = PYRA5A_REF[11];
+ break;
+ case 4:
+ coords[0] = PYRA5A_REF[12];
+ coords[1] = PYRA5A_REF[13];
+ coords[2] = PYRA5A_REF[14];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(-gc[0] + gc[1] - 1.0)*(-gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[1] = 0.25*(-gc[0] - gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[2] = 0.25*(+gc[0] + gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[3] = 0.25*(+gc[0] + gc[1] - 1.0)*(-gc[0] + gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[4] = gc[2];
+ SHAPE_FUN_MACRO_END;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(-gc[0] + gc[1] - 1.0)*(-gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[1] = 0.25*(-gc[0] - gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[2] = 0.25*(+gc[0] + gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[3] = 0.25*(+gc[0] + gc[1] - 1.0)*(-gc[0] + gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[4] = gc[2];
- SHAPE_FUN_MACRO_END;
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.5*(gc[0]+1)*(1.0 - gc[2]);
+ devFunValue[1] = -0.5*gc[1]*(1.0 - gc[2]);
+ devFunValue[2] = -0.25*(-gc[0] + gc[1] - 1.0)*(-gc[0] - gc[1] - 1.0);
+
+ devFunValue[3] = -0.5*gc[0]*(1.0 - gc[2]);
+ devFunValue[4] = 0.5*(gc[1]+1)*(1.0 - gc[2]);
+ devFunValue[5] = -0.25*(-gc[0] - gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0);
+
+ devFunValue[6] = 0.5*(gc[0]-1)*(1.0 - gc[2]);
+ devFunValue[7] = -0.5*gc[1]*(1.0 - gc[2]);
+ devFunValue[8] = -0.25*(+gc[0] + gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0);
+
+ devFunValue[9] = -0.5*gc[0]*(1.0 - gc[2]);
+ devFunValue[10] = 0.5*(1*gc[1]-1)*(1.0 - gc[2]);
+ devFunValue[11] = -0.25*(+gc[0] + gc[1] - 1.0)*(-gc[0] + gc[1] - 1.0);
+
+ devFunValue[12] = 0.0;
+ devFunValue[13] = 0.0;
+ devFunValue[14] = 1.0;
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
* Init Pyramid Reference coordinates ans Shape function.
void GaussInfo::pyra5bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(-gc[0] + gc[1] - 1.0)*(-gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[3] = 0.25*(-gc[0] - gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[2] = 0.25*(+gc[0] + gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[1] = 0.25*(+gc[0] + gc[1] - 1.0)*(-gc[0] + gc[1] - 1.0)*(1.0 - gc[2]);
- funValue[4] = gc[2];
- SHAPE_FUN_MACRO_END;
+ case 0:
+ coords[0] = PYRA5B_REF[0];
+ coords[1] = PYRA5B_REF[1];
+ coords[2] = PYRA5B_REF[2];
+ break;
+ case 1:
+ coords[0] = PYRA5B_REF[3];
+ coords[1] = PYRA5B_REF[4];
+ coords[2] = PYRA5B_REF[5];
+ break;
+ case 2:
+ coords[0] = PYRA5B_REF[6];
+ coords[1] = PYRA5B_REF[7];
+ coords[2] = PYRA5B_REF[8];
+ break;
+ case 3:
+ coords[0] = PYRA5B_REF[9];
+ coords[1] = PYRA5B_REF[10];
+ coords[2] = PYRA5B_REF[11];
+ break;
+ case 4:
+ coords[0] = PYRA5B_REF[12];
+ coords[1] = PYRA5B_REF[13];
+ coords[2] = PYRA5B_REF[14];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(-gc[0] + gc[1] - 1.0)*(-gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[3] = 0.25*(-gc[0] - gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[2] = 0.25*(+gc[0] + gc[1] - 1.0)*(+gc[0] - gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[1] = 0.25*(+gc[0] + gc[1] - 1.0)*(-gc[0] + gc[1] - 1.0)*(1.0 - gc[2]);
+ funValue[4] = gc[2];
+ SHAPE_FUN_MACRO_END;
}
/*!
* Case A.
*/
void GaussInfo::pyra13aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PYRA13A_REF[0];
+ coords[1] = PYRA13A_REF[1];
+ coords[2] = PYRA13A_REF[2];
+ break;
+ case 1:
+ coords[0] = PYRA13A_REF[3];
+ coords[1] = PYRA13A_REF[4];
+ coords[2] = PYRA13A_REF[5];
+ break;
+ case 2:
+ coords[0] = PYRA13A_REF[6];
+ coords[1] = PYRA13A_REF[7];
+ coords[2] = PYRA13A_REF[8];
+ break;
+ case 3:
+ coords[0] = PYRA13A_REF[9];
+ coords[1] = PYRA13A_REF[10];
+ coords[2] = PYRA13A_REF[11];
+ break;
+ case 4:
+ coords[0] = PYRA13A_REF[12];
+ coords[1] = PYRA13A_REF[13];
+ coords[2] = PYRA13A_REF[14];
+ break;
+ case 5:
+ coords[0] = PYRA13A_REF[15];
+ coords[1] = PYRA13A_REF[16];
+ coords[2] = PYRA13A_REF[17];
+ break;
+ case 6:
+ coords[0] = PYRA13A_REF[18];
+ coords[1] = PYRA13A_REF[19];
+ coords[2] = PYRA13A_REF[20];
+ break;
+ case 7:
+ coords[0] = PYRA13A_REF[21];
+ coords[1] = PYRA13A_REF[22];
+ coords[2] = PYRA13A_REF[23];
+ break;
+ case 8:
+ coords[0] = PYRA13A_REF[24];
+ coords[1] = PYRA13A_REF[25];
+ coords[2] = PYRA13A_REF[26];
+ break;
+ case 9:
+ coords[0] = PYRA13A_REF[27];
+ coords[1] = PYRA13A_REF[28];
+ coords[2] = PYRA13A_REF[29];
+ break;
+ case 10:
+ coords[0] = PYRA13A_REF[30];
+ coords[1] = PYRA13A_REF[31];
+ coords[2] = PYRA13A_REF[32];
+ break;
+ case 11:
+ coords[0] = PYRA13A_REF[33];
+ coords[1] = PYRA13A_REF[34];
+ coords[2] = PYRA13A_REF[35];
+ break;
+ case 12:
+ coords[0] = PYRA13A_REF[36];
+ coords[1] = PYRA13A_REF[37];
+ coords[2] = PYRA13A_REF[38];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)*
+ (gc[0] - 0.5)/(1.0 - gc[2]);
+ funValue[1] = 0.5*(-gc[0] - gc[1] + gc[2] - 1.0)*(+gc[0] - gc[1] + gc[2] - 1.0)*
+ (gc[1] - 0.5)/(1.0 - gc[2]);
+ funValue[2] = 0.5*(+gc[0] - gc[1] + gc[2] - 1.0)*(+gc[0] + gc[1] + gc[2] - 1.0)*
+ (-gc[0] - 0.5)/(1.0 - gc[2]);
+ funValue[3] = 0.5*(+gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)*
+ (-gc[1] - 0.5)/(1.0 - gc[2]);
+
+ funValue[4] = 2.0*gc[2]*(gc[2] - 0.5);
+
+ funValue[5] = 0.5*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)*
+ (gc[0] - gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
+ funValue[6] = 0.5*(-gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] - gc[1] + gc[2] - 1.0)*
+ (gc[0] + gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
+ funValue[7] = 0.5*(gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] + gc[1] + gc[2] - 1.0)*
+ (-gc[0] + gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
+ funValue[8] = 0.5*(gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)*
+ (-gc[0] - gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
+
+ funValue[9] = 0.5*gc[2]*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)/
+ (1.0 - gc[2]);
+ funValue[10] = 0.5*gc[2]*(-gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] - gc[1] + gc[2] - 1.0)/
+ (1.0 - gc[2]);
+ funValue[11] = 0.5*gc[2]*(gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] + gc[1] + gc[2] - 1.0)/
+ (1.0 - gc[2]);
+ funValue[12] = 0.5*gc[2]*(gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)/
+ (1.0 - gc[2]);
+ SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * Init Quadratic Pyramid Reference coordinates ans Shape function.
+ * Case B.
+ */
+void GaussInfo::pyra13bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PYRA13B_REF[0];
+ coords[1] = PYRA13B_REF[1];
+ coords[2] = PYRA13B_REF[2];
+ break;
+ case 1:
+ coords[0] = PYRA13B_REF[3];
+ coords[1] = PYRA13B_REF[4];
+ coords[2] = PYRA13B_REF[5];
+ break;
+ case 2:
+ coords[0] = PYRA13B_REF[6];
+ coords[1] = PYRA13B_REF[7];
+ coords[2] = PYRA13B_REF[8];
+ break;
+ case 3:
+ coords[0] = PYRA13B_REF[9];
+ coords[1] = PYRA13B_REF[10];
+ coords[2] = PYRA13B_REF[11];
+ break;
+ case 4:
+ coords[0] = PYRA13B_REF[12];
+ coords[1] = PYRA13B_REF[13];
+ coords[2] = PYRA13B_REF[14];
+ break;
+ case 5:
+ coords[0] = PYRA13B_REF[15];
+ coords[1] = PYRA13B_REF[16];
+ coords[2] = PYRA13B_REF[17];
+ break;
+ case 6:
+ coords[0] = PYRA13B_REF[18];
+ coords[1] = PYRA13B_REF[19];
+ coords[2] = PYRA13B_REF[20];
+ break;
+ case 7:
+ coords[0] = PYRA13B_REF[21];
+ coords[1] = PYRA13B_REF[22];
+ coords[2] = PYRA13B_REF[23];
+ break;
+ case 8:
+ coords[0] = PYRA13B_REF[24];
+ coords[1] = PYRA13B_REF[25];
+ coords[2] = PYRA13B_REF[26];
+ break;
+ case 9:
+ coords[0] = PYRA13B_REF[27];
+ coords[1] = PYRA13B_REF[28];
+ coords[2] = PYRA13B_REF[29];
+ break;
+ case 10:
+ coords[0] = PYRA13B_REF[30];
+ coords[1] = PYRA13B_REF[31];
+ coords[2] = PYRA13B_REF[32];
+ break;
+ case 11:
+ coords[0] = PYRA13B_REF[33];
+ coords[1] = PYRA13B_REF[34];
+ coords[2] = PYRA13B_REF[35];
+ break;
+ case 12:
+ coords[0] = PYRA13B_REF[36];
+ coords[1] = PYRA13B_REF[37];
+ coords[2] = PYRA13B_REF[38];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] =0.5*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-0.5)/(1.0-gc[2]);
+ funValue[1] =0.5*(+gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[1]-0.5)/(1.0-gc[2]);
+ funValue[2] =0.5*(+gc[0]-gc[1]+gc[2]-1.0)*(+gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-0.5)/(1.0-gc[2]);
+ funValue[3] =0.5*(-gc[0]-gc[1]+gc[2]-1.0)*(+gc[0]-gc[1]+gc[2]-1.0)*(gc[1]-0.5)/(1.0-gc[2]);
+
+ funValue[4] =2.*gc[2]*(gc[2]-0.5);
+
+ funValue[5] =-0.5*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[6] =-0.5*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[7] =-0.5*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[8] =-0.5*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+
+ funValue[9] =gc[2]*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[10]=gc[2]*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[11]=gc[2]*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+ funValue[12]=gc[2]*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
+
+ SHAPE_FUN_MACRO_END;
+}
+
+
+/*!
+ * Init Pentahedron Reference coordinates and Shape function.
+ * Case A.
+ */
+void GaussInfo::penta6aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PENTA6A_REF[0];
+ coords[1] = PENTA6A_REF[1];
+ coords[2] = PENTA6A_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA6A_REF[3];
+ coords[1] = PENTA6A_REF[4];
+ coords[2] = PENTA6A_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA6A_REF[6];
+ coords[1] = PENTA6A_REF[7];
+ coords[2] = PENTA6A_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA6A_REF[9];
+ coords[1] = PENTA6A_REF[10];
+ coords[2] = PENTA6A_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA6A_REF[12];
+ coords[1] = PENTA6A_REF[13];
+ coords[2] = PENTA6A_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA6A_REF[15];
+ coords[1] = PENTA6A_REF[16];
+ coords[2] = PENTA6A_REF[17];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0]);
+ funValue[1] = 0.5*gc[2]*(1.0 - gc[0]);
+ funValue[2] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+
+ funValue[3] = 0.5*gc[1]*(gc[0] + 1.0);
+ funValue[4] = 0.5*gc[2]*(gc[0] + 1.0);
+ funValue[5] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.5*gc[1]*(-1.0);
+ devFunValue[1] = 0.5*(1.0 - gc[0]);
+ devFunValue[2] = 0.0;
+
+ devFunValue[3] = 0.5*gc[2]*(-1.0);
+ devFunValue[4] = 0.0;
+ devFunValue[5] = 0.5*(1.0 - gc[0]);
+
+ devFunValue[6] = 0.5*(1.0 - gc[1] - gc[2])*(-1.0);
+ devFunValue[7] = 0.5*(-1.0)*(1.0 - gc[0]);
+ devFunValue[8] = 0.5*(-1.0)*(1.0 - gc[0]);
+
+ devFunValue[9] = 0.5*gc[1];
+ devFunValue[10] = 0.5*(gc[0] + 1.0);
+ devFunValue[11] = 0.0;
+
+ devFunValue[12] = 0.5*gc[2];
+ devFunValue[13] = 0.0;
+ devFunValue[14] = 0.5*(gc[0] + 1.0);
+
+ devFunValue[15] = 0.5*(1.0 - gc[1] - gc[2]);
+ devFunValue[16] = 0.5*(-1.0)*(1.0 + gc[0]);
+ devFunValue[17] = 0.5*(-1.0)*(1.0 + gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * Init Pentahedron Reference coordinates and Shape function.
+ * Case B.
+ */
+void GaussInfo::penta6bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PENTA6B_REF[0];
+ coords[1] = PENTA6B_REF[1];
+ coords[2] = PENTA6B_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA6B_REF[3];
+ coords[1] = PENTA6B_REF[4];
+ coords[2] = PENTA6B_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA6B_REF[6];
+ coords[1] = PENTA6B_REF[7];
+ coords[2] = PENTA6B_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA6B_REF[9];
+ coords[1] = PENTA6B_REF[10];
+ coords[2] = PENTA6B_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA6B_REF[12];
+ coords[1] = PENTA6B_REF[13];
+ coords[2] = PENTA6B_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA6B_REF[15];
+ coords[1] = PENTA6B_REF[16];
+ coords[2] = PENTA6B_REF[17];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0]);
+ funValue[2] = 0.5*gc[2]*(1.0 - gc[0]);
+ funValue[1] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ funValue[3] = 0.5*gc[1]*(gc[0] + 1.0);
+ funValue[5] = 0.5*gc[2]*(gc[0] + 1.0);
+ funValue[4] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * This shapefunc map is same as degenerated tria3aInit
+ */
+void GaussInfo::penta6DegTria3aInit()
{
LOCAL_COORD_MACRO_BEGIN;
case 0:
- coords[0] = 1.0;
- coords[1] = 0.0;
+ coords[0] = -1.0;
+ coords[1] = 1.0;
coords[2] = 0.0;
break;
case 1:
- coords[0] = 0.0;
- coords[1] = 1.0;
+ coords[0] = -1.0;
+ coords[1] = -1.0;
coords[2] = 0.0;
break;
case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
+ coords[0] = 1.0;
+ coords[1] = -1.0;
coords[2] = 0.0;
break;
case 3:
coords[0] = 0.0;
- coords[1] = -1.0;
+ coords[1] = 0.0;
coords[2] = 0.0;
break;
case 4:
coords[0] = 0.0;
coords[1] = 0.0;
- coords[2] = 1.0;
- break;
-
- case 5:
- coords[0] = 0.5;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 6:
- coords[0] = -0.5;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = -0.5;
- coords[1] = -0.5;
- coords[2] = 0.0;
- break;
- case 8:
- coords[0] = 0.5;
- coords[1] = -0.5;
coords[2] = 0.0;
break;
- case 9:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 10:
+ case 5:
coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 11:
- coords[0] = -0.5;
coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 12:
- coords[0] = 0.0;
- coords[1] = -0.5;
- coords[2] = 0.5;
+ coords[2] = 0.0;
break;
LOCAL_COORD_MACRO_END;
SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)*
- (gc[0] - 0.5)/(1.0 - gc[2]);
- funValue[1] = 0.5*(-gc[0] - gc[1] + gc[2] - 1.0)*(+gc[0] - gc[1] + gc[2] - 1.0)*
- (gc[1] - 0.5)/(1.0 - gc[2]);
- funValue[2] = 0.5*(+gc[0] - gc[1] + gc[2] - 1.0)*(+gc[0] + gc[1] + gc[2] - 1.0)*
- (-gc[0] - 0.5)/(1.0 - gc[2]);
- funValue[3] = 0.5*(+gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)*
- (-gc[1] - 0.5)/(1.0 - gc[2]);
-
- funValue[4] = 2.0*gc[2]*(gc[2] - 0.5);
-
- funValue[5] = 0.5*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)*
- (gc[0] - gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
- funValue[6] = 0.5*(-gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] - gc[1] + gc[2] - 1.0)*
- (gc[0] + gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
- funValue[7] = 0.5*(gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] + gc[1] + gc[2] - 1.0)*
- (-gc[0] + gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
- funValue[8] = 0.5*(gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)*
- (-gc[0] - gc[1] + gc[2] - 1.0)/(1.0 - gc[2]);
-
- funValue[9] = 0.5*gc[2]*(-gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] - gc[1] + gc[2] - 1.0)/
- (1.0 - gc[2]);
- funValue[10] = 0.5*gc[2]*(-gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] - gc[1] + gc[2] - 1.0)/
- (1.0 - gc[2]);
- funValue[11] = 0.5*gc[2]*(gc[0] - gc[1] + gc[2] - 1.0)*(gc[0] + gc[1] + gc[2] - 1.0)/
- (1.0 - gc[2]);
- funValue[12] = 0.5*gc[2]*(gc[0] + gc[1] + gc[2] - 1.0)*(-gc[0] + gc[1] + gc[2] - 1.0)/
- (1.0 - gc[2]);
+ funValue[0] = 0.5*(1.0 + gc[1]);
+ funValue[1] = -0.5*(gc[0] + gc[1]);
+ funValue[2] = 0.5*(1.0 + gc[0]);
+ funValue[3] = 0.;
+ funValue[4] = 0.;
+ funValue[5] = 0.;
SHAPE_FUN_MACRO_END;
}
/*!
- * Init Quadratic Pyramid Reference coordinates ans Shape function.
- * Case B.
+ * This shapefunc map is same as degenerated tria3bInit
*/
-void GaussInfo::pyra13bInit()
+void GaussInfo::penta6DegTria3bInit()
{
LOCAL_COORD_MACRO_BEGIN;
case 0:
- coords[0] = 1.0;
+ coords[0] = 0.0;
coords[1] = 0.0;
coords[2] = 0.0;
break;
case 1:
- coords[0] = 0.0;
- coords[1] = -1.0;
+ coords[0] = 1.0;
+ coords[1] = 0.0;
coords[2] = 0.0;
break;
case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
+ coords[0] = 0.0;
+ coords[1] = 1.0;
coords[2] = 0.0;
break;
case 3:
coords[0] = 0.0;
- coords[1] = 1.0;
+ coords[1] = 0.0;
coords[2] = 0.0;
break;
case 4:
coords[0] = 0.0;
coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = 0.5;
- coords[1] = -0.5;
- coords[2] = 0.0;
- break;
- case 6:
- coords[0] = -0.5;
- coords[1] = -0.5;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = -0.5;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 8:
- coords[0] = 0.5;
- coords[1] = 0.5;
coords[2] = 0.0;
break;
- case 9:
- coords[0] = 0.5;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 10:
+ case 5:
coords[0] = 0.0;
- coords[1] = -0.5;
- coords[2] = 0.5;
- break;
- case 11:
- coords[0] = -0.5;
coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 12:
- coords[0] = 0.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
+ coords[2] = 0.0;
break;
LOCAL_COORD_MACRO_END;
SHAPE_FUN_MACRO_BEGIN;
- funValue[0] =0.5*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-0.5)/(1.0-gc[2]);
- funValue[1] =0.5*(+gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[1]-0.5)/(1.0-gc[2]);
- funValue[2] =0.5*(+gc[0]-gc[1]+gc[2]-1.0)*(+gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-0.5)/(1.0-gc[2]);
- funValue[3] =0.5*(-gc[0]-gc[1]+gc[2]-1.0)*(+gc[0]-gc[1]+gc[2]-1.0)*(gc[1]-0.5)/(1.0-gc[2]);
-
- funValue[4] =2.*gc[2]*(gc[2]-0.5);
-
- funValue[5] =-0.5*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[6] =-0.5*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[7] =-0.5*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[8] =-0.5*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
-
- funValue[9] =gc[2]*(-gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[10]=gc[2]*(gc[0]+gc[1]+gc[2]-1.0)*(-gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[11]=gc[2]*(gc[0]-gc[1]+gc[2]-1.0)*(gc[0]+gc[1]+gc[2]-1.0)/(1.0-gc[2]);
- funValue[12]=gc[2]*(-gc[0]-gc[1]+gc[2]-1.0)*(gc[0]-gc[1]+gc[2]-1.0)/(1.0-gc[2]);
-
+ funValue[0] = 1.0 - gc[0] - gc[1];
+ funValue[1] = gc[0];
+ funValue[2] = gc[1];
+ funValue[3] = 0.;
+ funValue[4] = 0.;
+ funValue[5] = 0.;
SHAPE_FUN_MACRO_END;
}
-
/*!
* Init Pentahedron Reference coordinates and Shape function.
* Case A.
*/
-void GaussInfo::penta6aInit()
+void GaussInfo::penta15aInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -0.0;
- coords[2] = 1.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = PENTA15A_REF[0];
+ coords[1] = PENTA15A_REF[1];
+ coords[2] = PENTA15A_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA15A_REF[3];
+ coords[1] = PENTA15A_REF[4];
+ coords[2] = PENTA15A_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA15A_REF[6];
+ coords[1] = PENTA15A_REF[7];
+ coords[2] = PENTA15A_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA15A_REF[9];
+ coords[1] = PENTA15A_REF[10];
+ coords[2] = PENTA15A_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA15A_REF[12];
+ coords[1] = PENTA15A_REF[13];
+ coords[2] = PENTA15A_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA15A_REF[15];
+ coords[1] = PENTA15A_REF[16];
+ coords[2] = PENTA15A_REF[17];
+ break;
+ case 6:
+ coords[0] = PENTA15A_REF[18];
+ coords[1] = PENTA15A_REF[19];
+ coords[2] = PENTA15A_REF[20];
+ break;
+ case 7:
+ coords[0] = PENTA15A_REF[21];
+ coords[1] = PENTA15A_REF[22];
+ coords[2] = PENTA15A_REF[23];
+ break;
+ case 8:
+ coords[0] = PENTA15A_REF[24];
+ coords[1] = PENTA15A_REF[25];
+ coords[2] = PENTA15A_REF[26];
+ break;
+ case 9:
+ coords[0] = PENTA15A_REF[27];
+ coords[1] = PENTA15A_REF[28];
+ coords[2] = PENTA15A_REF[29];
+ break;
+ case 10:
+ coords[0] = PENTA15A_REF[30];
+ coords[1] = PENTA15A_REF[31];
+ coords[2] = PENTA15A_REF[32];
+ break;
+ case 11:
+ coords[0] = PENTA15A_REF[33];
+ coords[1] = PENTA15A_REF[34];
+ coords[2] = PENTA15A_REF[35];
+ break;
+ case 12:
+ coords[0] = PENTA15A_REF[36];
+ coords[1] = PENTA15A_REF[37];
+ coords[2] = PENTA15A_REF[38];
+ break;
+ case 13:
+ coords[0] = PENTA15A_REF[39];
+ coords[1] = PENTA15A_REF[40];
+ coords[2] = PENTA15A_REF[41];
+ break;
+ case 14:
+ coords[0] = PENTA15A_REF[42];
+ coords[1] = PENTA15A_REF[43];
+ coords[2] = PENTA15A_REF[44];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
+ funValue[1] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
+ funValue[2] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
+ funValue[4] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
+ funValue[5] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[6] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
+ funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ funValue[8] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+
+ funValue[9] = gc[1]*(1.0 - gc[0]*gc[0]);
+ funValue[10] = gc[2]*(1.0 - gc[0]*gc[0]);
+ funValue[11] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
+
+ funValue[12] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
+ funValue[13] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ funValue[14] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*gc[1]*(1.0 - gc[0]);
- funValue[1] = 0.5*gc[2]*(1.0 - gc[0]);
- funValue[2] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ DEV_SHAPE_FUN_MACRO_BEGIN;
- funValue[3] = 0.5*gc[1]*(gc[0] + 1.0);
- funValue[4] = 0.5*gc[2]*(gc[0] + 1.0);
- funValue[5] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
+ devFunValue[0] = 0.5*gc[1]*(2 * gc[0] - 2 * gc[1] + 1.0 );
+ devFunValue[1] = (1.0 - gc[0])*(2.0*gc[1] -1 - 0.5*gc[0]);
+ devFunValue[2] = 0.0;
+
+ devFunValue[3] = 0.5*gc[2]*(2 * gc[0] - 2 * gc[2] + 1.0 );
+ devFunValue[4] = 0.0;
+ devFunValue[5] = (1.0 - gc[0])*(2.0*gc[2] -1 - 0.5*gc[0]);
+
+ devFunValue[6] = 0.5*(1.0 - gc[1] - gc[2])*(2*gc[0] -1.0 + 2*gc[1] + 2*gc[2]);
+ devFunValue[7] = 0.5*(gc[0] - 1.0)* ( -4*gc[1] -gc[0] -4*gc[2] + 2.0);
+ devFunValue[8] = 0.5*(gc[0] - 1.0)* ( -4*gc[2] -4*gc[1] -gc[0] + 2.0);
+
+ devFunValue[9] = 0.5*gc[1]*( 2*gc[0] + 2*gc[1] -1.0 );
+ devFunValue[10] = (1.0 + gc[0])*(2.0*gc[1] - 1.0 + 0.5*gc[0]);
+ devFunValue[11] = 0.0;
+
+ devFunValue[12] = 0.5*gc[2]*( 2*gc[0] + 2*gc[2] -1.0 );
+ devFunValue[13] = 0.0;
+ devFunValue[14] = (1.0 + gc[0])*(2.0*gc[2] - 1.0 + 0.5*gc[0]);
+
+ devFunValue[15] = 0.5*(1.0 - gc[1] - gc[2])*(2 * gc[0] - 2.0 * gc[1] - 2.0 * gc[2] + 1.0 );
+ devFunValue[16] = 0.5*(-gc[0] - 1.0) * (-4.0*gc[1] + gc[0] - 4.0 * gc[2] + 2.0) ;
+ devFunValue[17] = 0.5*(-gc[0] - 1.0) * (-4.0*gc[2] - 4*gc[1] + gc[0] + 2.0 );
+
+ devFunValue[18] = - 2.0*gc[1]*gc[2];
+ devFunValue[19] = 2.0*gc[2]*(1.0 - gc[0]);
+ devFunValue[20] = 2.0*gc[1]*(1.0 - gc[0]);
+
+ devFunValue[21] = -2.0*gc[2]*(1.0 - gc[1] - gc[2]);
+ devFunValue[22] = -2.0*gc[2]*(1.0 - gc[0]);
+ devFunValue[23] = (-4*gc[2]-2*gc[1]+2.0)*(1.0 - gc[0]);
+
+ devFunValue[24] = -2.0*gc[1]*(1.0 - gc[1] - gc[2]);
+ devFunValue[25] = (-4*gc[1]-2*gc[2]+2)*(1.0 - gc[0]);
+ devFunValue[26] = -2.0*gc[1]*(1.0 - gc[0]);
+
+ devFunValue[27] = gc[1]*(- 2.0 * gc[0]);
+ devFunValue[28] = (1.0 - gc[0]*gc[0]);
+ devFunValue[29] = 0.0;
+
+ devFunValue[30] = -2.0 * gc[2] *gc[0];
+ devFunValue[31] = 0.0;
+ devFunValue[32] = (1.0 - gc[0]*gc[0]);
+
+ devFunValue[33] = -2.0 * (1.0 - gc[1] - gc[2]) * gc[0];
+ devFunValue[34] = -1.0*(1.0 - gc[0]*gc[0]);
+ devFunValue[35] = -1.0*(1.0 - gc[0]*gc[0]);
+
+ devFunValue[36] = 2.0*gc[1]*gc[2];
+ devFunValue[37] = 2.0*gc[2]*(1.0 + gc[0]);
+ devFunValue[38] = 2.0*gc[1]*(1.0 + gc[0]);
+
+ devFunValue[39] = 2.0*gc[2]*(1.0 - gc[1] - gc[2]);
+ devFunValue[40] = -2.0*gc[2]*(1.0 + gc[0]);
+ devFunValue[41] = (2.0 - 2.0 * gc[1] - 4.0 * gc[2])*(1.0 + gc[0]);
+
+ devFunValue[42] = 2.0*gc[1]*(1.0 - gc[1] - gc[2]);
+ devFunValue[43] = (2.0 - 4*gc[1] - 2*gc[2])*(1.0 + gc[0]);
+ devFunValue[44] = -2.0*gc[1]*(1.0 + gc[0]);
+
+ DEV_SHAPE_FUN_MACRO_END;
}
/*!
- * Init Pentahedron Reference coordinates and Shape function.
+ * Init Qaudratic Pentahedron Reference coordinates and Shape function.
* Case B.
*/
-void GaussInfo::penta6bInit()
+void GaussInfo::penta15bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = -0.0;
- coords[2] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 4:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
+ case 0:
+ coords[0] = PENTA15B_REF[0];
+ coords[1] = PENTA15B_REF[1];
+ coords[2] = PENTA15B_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA15B_REF[3];
+ coords[1] = PENTA15B_REF[4];
+ coords[2] = PENTA15B_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA15B_REF[6];
+ coords[1] = PENTA15B_REF[7];
+ coords[2] = PENTA15B_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA15B_REF[9];
+ coords[1] = PENTA15B_REF[10];
+ coords[2] = PENTA15B_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA15B_REF[12];
+ coords[1] = PENTA15B_REF[13];
+ coords[2] = PENTA15B_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA15B_REF[15];
+ coords[1] = PENTA15B_REF[16];
+ coords[2] = PENTA15B_REF[17];
+ break;
+ case 6:
+ coords[0] = PENTA15B_REF[18];
+ coords[1] = PENTA15B_REF[19];
+ coords[2] = PENTA15B_REF[20];
+ break;
+ case 7:
+ coords[0] = PENTA15B_REF[21];
+ coords[1] = PENTA15B_REF[22];
+ coords[2] = PENTA15B_REF[23];
+ break;
+ case 8:
+ coords[0] = PENTA15B_REF[24];
+ coords[1] = PENTA15B_REF[25];
+ coords[2] = PENTA15B_REF[26];
+ break;
+ case 9:
+ coords[0] = PENTA15B_REF[27];
+ coords[1] = PENTA15B_REF[28];
+ coords[2] = PENTA15B_REF[29];
+ break;
+ case 10:
+ coords[0] = PENTA15B_REF[30];
+ coords[1] = PENTA15B_REF[31];
+ coords[2] = PENTA15B_REF[32];
+ break;
+ case 11:
+ coords[0] = PENTA15B_REF[33];
+ coords[1] = PENTA15B_REF[34];
+ coords[2] = PENTA15B_REF[35];
+ break;
+ case 12:
+ coords[0] = PENTA15B_REF[36];
+ coords[1] = PENTA15B_REF[37];
+ coords[2] = PENTA15B_REF[38];
+ break;
+ case 13:
+ coords[0] = PENTA15B_REF[39];
+ coords[1] = PENTA15B_REF[40];
+ coords[2] = PENTA15B_REF[41];
+ break;
+ case 14:
+ coords[0] = PENTA15B_REF[42];
+ coords[1] = PENTA15B_REF[43];
+ coords[2] = PENTA15B_REF[44];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
+ funValue[2] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
+ funValue[1] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
+ funValue[5] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
+ funValue[4] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[8] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
+ funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ funValue[6] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+
+ funValue[12] = gc[1]*(1.0 - gc[0]*gc[0]);
+ funValue[14] = gc[2]*(1.0 - gc[0]*gc[0]);
+ funValue[13] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
+
+ funValue[11] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
+ funValue[10] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ funValue[9] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ SHAPE_FUN_MACRO_END;
+}
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*gc[1]*(1.0 - gc[0]);
- funValue[2] = 0.5*gc[2]*(1.0 - gc[0]);
- funValue[1] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
- funValue[3] = 0.5*gc[1]*(gc[0] + 1.0);
- funValue[5] = 0.5*gc[2]*(gc[0] + 1.0);
- funValue[4] = 0.5*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
+void GaussInfo::penta18aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PENTA18A_REF[0];
+ coords[1] = PENTA18A_REF[1];
+ coords[2] = PENTA18A_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA18A_REF[3];
+ coords[1] = PENTA18A_REF[4];
+ coords[2] = PENTA18A_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA18A_REF[6];
+ coords[1] = PENTA18A_REF[7];
+ coords[2] = PENTA18A_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA18A_REF[9];
+ coords[1] = PENTA18A_REF[10];
+ coords[2] = PENTA18A_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA18A_REF[12];
+ coords[1] = PENTA18A_REF[13];
+ coords[2] = PENTA18A_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA18A_REF[15];
+ coords[1] = PENTA18A_REF[16];
+ coords[2] = PENTA18A_REF[17];
+ break;
+ case 6:
+ coords[0] = PENTA18A_REF[18];
+ coords[1] = PENTA18A_REF[19];
+ coords[2] = PENTA18A_REF[20];
+ break;
+ case 7:
+ coords[0] = PENTA18A_REF[21];
+ coords[1] = PENTA18A_REF[22];
+ coords[2] = PENTA18A_REF[23];
+ break;
+ case 8:
+ coords[0] = PENTA18A_REF[24];
+ coords[1] = PENTA18A_REF[25];
+ coords[2] = PENTA18A_REF[26];
+ break;
+ case 9:
+ coords[0] = PENTA18A_REF[27];
+ coords[1] = PENTA18A_REF[28];
+ coords[2] = PENTA18A_REF[29];
+ break;
+ case 10:
+ coords[0] = PENTA18A_REF[30];
+ coords[1] = PENTA18A_REF[31];
+ coords[2] = PENTA18A_REF[32];
+ break;
+ case 11:
+ coords[0] = PENTA18A_REF[33];
+ coords[1] = PENTA18A_REF[34];
+ coords[2] = PENTA18A_REF[35];
+ break;
+ case 12:
+ coords[0] = PENTA18A_REF[36];
+ coords[1] = PENTA18A_REF[37];
+ coords[2] = PENTA18A_REF[38];
+ break;
+ case 13:
+ coords[0] = PENTA18A_REF[39];
+ coords[1] = PENTA18A_REF[40];
+ coords[2] = PENTA18A_REF[41];
+ break;
+ case 14:
+ coords[0] = PENTA18A_REF[42];
+ coords[1] = PENTA18A_REF[43];
+ coords[2] = PENTA18A_REF[44];
+ break;
+ case 15:
+ coords[0] = PENTA18A_REF[45];
+ coords[1] = PENTA18A_REF[46];
+ coords[2] = PENTA18A_REF[47];
+ break;
+ case 16:
+ coords[0] = PENTA18A_REF[48];
+ coords[1] = PENTA18A_REF[49];
+ coords[2] = PENTA18A_REF[50];
+ break;
+ case 17:
+ coords[0] = PENTA18A_REF[51];
+ coords[1] = PENTA18A_REF[52];
+ coords[2] = PENTA18A_REF[53];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
+ funValue[1] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
+ funValue[2] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
+ funValue[4] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
+ funValue[5] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[6] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
+ funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ funValue[8] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+
+ funValue[9] = gc[1]*(1.0 - gc[0]*gc[0]);
+ funValue[10] = gc[2]*(1.0 - gc[0]*gc[0]);
+ funValue[11] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
+
+ funValue[12] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
+ funValue[13] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ funValue[14] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+
+ funValue[15] = 4.0*gc[1]*gc[2]*(1.0 - gc[0] * gc[0]);
+ funValue[16] = 4.0*gc[2]*(gc[0] * gc[0] - 1) * ( gc[2] + gc[1] - 1);
+ funValue[17] = 4.0*gc[1]*(gc[0] * gc[0] - 1) * ( gc[2] + gc[1] - 1);
+ SHAPE_FUN_MACRO_END;
+}
+
+void GaussInfo::penta18bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = PENTA18B_REF[0];
+ coords[1] = PENTA18B_REF[1];
+ coords[2] = PENTA18B_REF[2];
+ break;
+ case 1:
+ coords[0] = PENTA18B_REF[3];
+ coords[1] = PENTA18B_REF[4];
+ coords[2] = PENTA18B_REF[5];
+ break;
+ case 2:
+ coords[0] = PENTA18B_REF[6];
+ coords[1] = PENTA18B_REF[7];
+ coords[2] = PENTA18B_REF[8];
+ break;
+ case 3:
+ coords[0] = PENTA18B_REF[9];
+ coords[1] = PENTA18B_REF[10];
+ coords[2] = PENTA18B_REF[11];
+ break;
+ case 4:
+ coords[0] = PENTA18B_REF[12];
+ coords[1] = PENTA18B_REF[13];
+ coords[2] = PENTA18B_REF[14];
+ break;
+ case 5:
+ coords[0] = PENTA18B_REF[15];
+ coords[1] = PENTA18B_REF[16];
+ coords[2] = PENTA18B_REF[17];
+ break;
+ case 6:
+ coords[0] = PENTA18B_REF[18];
+ coords[1] = PENTA18B_REF[19];
+ coords[2] = PENTA18B_REF[20];
+ break;
+ case 7:
+ coords[0] = PENTA18B_REF[21];
+ coords[1] = PENTA18B_REF[22];
+ coords[2] = PENTA18B_REF[23];
+ break;
+ case 8:
+ coords[0] = PENTA18B_REF[24];
+ coords[1] = PENTA18B_REF[25];
+ coords[2] = PENTA18B_REF[26];
+ break;
+ case 9:
+ coords[0] = PENTA18B_REF[27];
+ coords[1] = PENTA18B_REF[28];
+ coords[2] = PENTA18B_REF[29];
+ break;
+ case 10:
+ coords[0] = PENTA18B_REF[30];
+ coords[1] = PENTA18B_REF[31];
+ coords[2] = PENTA18B_REF[32];
+ break;
+ case 11:
+ coords[0] = PENTA18B_REF[33];
+ coords[1] = PENTA18B_REF[34];
+ coords[2] = PENTA18B_REF[35];
+ break;
+ case 12:
+ coords[0] = PENTA18B_REF[36];
+ coords[1] = PENTA18B_REF[37];
+ coords[2] = PENTA18B_REF[38];
+ break;
+ case 13:
+ coords[0] = PENTA18B_REF[39];
+ coords[1] = PENTA18B_REF[40];
+ coords[2] = PENTA18B_REF[41];
+ break;
+ case 14:
+ coords[0] = PENTA18B_REF[42];
+ coords[1] = PENTA18B_REF[43];
+ coords[2] = PENTA18B_REF[44];
+ break;
+ case 15:
+ coords[0] = PENTA18B_REF[45];
+ coords[1] = PENTA18B_REF[46];
+ coords[2] = PENTA18B_REF[47];
+ break;
+ case 16:
+ coords[0] = PENTA18B_REF[48];
+ coords[1] = PENTA18B_REF[49];
+ coords[2] = PENTA18B_REF[50];
+ break;
+ case 17:
+ coords[0] = PENTA18B_REF[51];
+ coords[1] = PENTA18B_REF[52];
+ coords[2] = PENTA18B_REF[53];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
+ funValue[2] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
+ funValue[1] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
+ funValue[5] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
+ funValue[4] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
+
+ funValue[8] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
+ funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+ funValue[6] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
+
+ funValue[12] = gc[1]*(1.0 - gc[0]*gc[0]);
+ funValue[14] = gc[2]*(1.0 - gc[0]*gc[0]);
+ funValue[13] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
+
+ funValue[11] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
+ funValue[10] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+ funValue[9] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
+
+ funValue[17] = 4.0*gc[1]*gc[2]*(1.0 - gc[0] * gc[0]);
+ funValue[16] = 4.0*gc[2]*(gc[0] * gc[0] - 1) * ( gc[2] + gc[1] - 1);
+ funValue[15] = 4.0*gc[1]*(gc[0] * gc[0] - 1) * ( gc[2] + gc[1] - 1);
+ SHAPE_FUN_MACRO_END;
}
/*!
- * This shapefunc map is same as degenerated tria3aInit
+ * Init Hehahedron Reference coordinates and Shape function.
+ * Case A.
*/
-void GaussInfo::penta6DegTria3aInit()
+void GaussInfo::hexa8aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = HEXA8A_REF[0];
+ coords[1] = HEXA8A_REF[1];
+ coords[2] = HEXA8A_REF[2];
+ break;
+ case 1:
+ coords[0] = HEXA8A_REF[3];
+ coords[1] = HEXA8A_REF[4];
+ coords[2] = HEXA8A_REF[5];
+ break;
+ case 2:
+ coords[0] = HEXA8A_REF[6];
+ coords[1] = HEXA8A_REF[7];
+ coords[2] = HEXA8A_REF[8];
+ break;
+ case 3:
+ coords[0] = HEXA8A_REF[9];
+ coords[1] = HEXA8A_REF[10];
+ coords[2] = HEXA8A_REF[11];
+ break;
+ case 4:
+ coords[0] = HEXA8A_REF[12];
+ coords[1] = HEXA8A_REF[13];
+ coords[2] = HEXA8A_REF[14];
+ break;
+ case 5:
+ coords[0] = HEXA8A_REF[15];
+ coords[1] = HEXA8A_REF[16];
+ coords[2] = HEXA8A_REF[17];
+ break;
+ case 6:
+ coords[0] = HEXA8A_REF[18];
+ coords[1] = HEXA8A_REF[19];
+ coords[2] = HEXA8A_REF[20];
+ break;
+ case 7:
+ coords[0] = HEXA8A_REF[21];
+ coords[1] = HEXA8A_REF[22];
+ coords[2] = HEXA8A_REF[23];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[1] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ funValue[3] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+
+ funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ funValue[7] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.125 * (-1.0) * (1.0 - gc[1])*(1.0 - gc[2]);
+ devFunValue[1] = 0.125 * (1.0 - gc[0]) * (-1.0) * (1.0 - gc[2]);
+ devFunValue[2] = 0.125 * (1.0 - gc[0]) * (1.0 - gc[1]) *(-1.0);
+
+ devFunValue[3] = 0.125*(1.0 - gc[1])*(1.0 - gc[2]);
+ devFunValue[4] = 0.125*(1.0 + gc[0])*(-1.0)*(1.0 - gc[2]);
+ devFunValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(-1.0);
+
+ devFunValue[6] = 0.125*(1.0 + gc[1])*(1.0 - gc[2]);
+ devFunValue[7] = 0.125*(1.0 + gc[0])*(1.0 - gc[2]);
+ devFunValue[8] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(-1.0);
+
+ devFunValue[9 ] = 0.125*(-1.0)*(1.0 + gc[1])*(1.0 - gc[2]);
+ devFunValue[10] = 0.125*(1.0 - gc[0])*(1.0 - gc[2]);
+ devFunValue[11] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(-1.0);
+
+ devFunValue[12] = 0.125*(-1.0)*(1.0 - gc[1])*(1.0 + gc[2]);
+ devFunValue[13] = 0.125*(1.0 - gc[0])*(-1.0)*(1.0 + gc[2]);
+ devFunValue[14] = 0.125*(1.0 - gc[0])*(1.0 - gc[1]);
+
+ devFunValue[15] = 0.125*(1.0 - gc[1])*(1.0 + gc[2]);
+ devFunValue[16] = 0.125*(1.0 + gc[0])*(-1.0)*(1.0 + gc[2]);
+ devFunValue[17] = 0.125*(1.0 + gc[0])*(1.0 - gc[1]);
+
+ devFunValue[18] = 0.125*(1.0 + gc[1])*(1.0 + gc[2]);
+ devFunValue[19] = 0.125*(1.0 + gc[0])*(1.0 + gc[2]);
+ devFunValue[20] = 0.125*(1.0 + gc[0])*(1.0 + gc[1]);
+
+ devFunValue[21] = 0.125*(-1.0)*(1.0 + gc[1])*(1.0 + gc[2]);
+ devFunValue[22] = 0.125*(1.0 - gc[0])*(1.0 + gc[2]);
+ devFunValue[23] = 0.125*(1.0 - gc[0])*(1.0 + gc[1]);
+
+ DEV_SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * Init Hehahedron Reference coordinates and Shape function.
+ * Case B.
+ */
+void GaussInfo::hexa8bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = HEXA8B_REF[0];
+ coords[1] = HEXA8B_REF[1];
+ coords[2] = HEXA8B_REF[2];
+ break;
+ case 1:
+ coords[0] = HEXA8B_REF[3];
+ coords[1] = HEXA8B_REF[4];
+ coords[2] = HEXA8B_REF[5];
+ break;
+ case 2:
+ coords[0] = HEXA8B_REF[6];
+ coords[1] = HEXA8B_REF[7];
+ coords[2] = HEXA8B_REF[8];
+ break;
+ case 3:
+ coords[0] = HEXA8B_REF[9];
+ coords[1] = HEXA8B_REF[10];
+ coords[2] = HEXA8B_REF[11];
+ break;
+ case 4:
+ coords[0] = HEXA8B_REF[12];
+ coords[1] = HEXA8B_REF[13];
+ coords[2] = HEXA8B_REF[14];
+ break;
+ case 5:
+ coords[0] = HEXA8B_REF[15];
+ coords[1] = HEXA8B_REF[16];
+ coords[2] = HEXA8B_REF[17];
+ break;
+ case 6:
+ coords[0] = HEXA8B_REF[18];
+ coords[1] = HEXA8B_REF[19];
+ coords[2] = HEXA8B_REF[20];
+ break;
+ case 7:
+ coords[0] = HEXA8B_REF[21];
+ coords[1] = HEXA8B_REF[22];
+ coords[2] = HEXA8B_REF[23];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[3] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ funValue[1] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+
+ funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[7] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ funValue[5] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * This shapefunc map is same as degenerated quad4bInit
+ */
+void GaussInfo::hexa8DegQuad4aInit()
{
LOCAL_COORD_MACRO_BEGIN;
case 0:
coords[2] = 0.0;
break;
case 3:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*(1.0 + gc[1]);
- funValue[1] = -0.5*(gc[0] + gc[1]);
- funValue[2] = 0.5*(1.0 + gc[0]);
- funValue[3] = 0.;
- funValue[4] = 0.;
- funValue[5] = 0.;
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * This shapefunc map is same as degenerated tria3bInit
- */
-void GaussInfo::penta6DegTria3bInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 1.0 - gc[0] - gc[1];
- funValue[1] = gc[0];
- funValue[2] = gc[1];
- funValue[3] = 0.;
- funValue[4] = 0.;
- funValue[5] = 0.;
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * Init Pentahedron Reference coordinates and Shape function.
- * Case A.
- */
-void GaussInfo::penta15aInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -0.0;
- coords[2] = 1.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
-
- case 6:
- coords[0] = -1.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 8:
- coords[0] = -1.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 9:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 10:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 11:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 12:
- coords[0] = 1.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 13:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 14:
- coords[0] = 1.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
- funValue[1] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
- funValue[2] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
-
- funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
- funValue[4] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
- funValue[5] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
-
- funValue[6] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
- funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
- funValue[8] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
-
- funValue[9] = gc[1]*(1.0 - gc[0]*gc[0]);
- funValue[10] = gc[2]*(1.0 - gc[0]*gc[0]);
- funValue[11] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
-
- funValue[12] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
- funValue[13] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- funValue[14] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * Init Qaudratic Pentahedron Reference coordinates and Shape function.
- * Case B.
- */
-void GaussInfo::penta15bInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = -1.0;
- coords[1] = -0.0;
- coords[2] = 1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 4:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
-
- case 8:
- coords[0] = -1.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 6:
- coords[0] = -1.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- case 12:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 14:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 13:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 11:
- coords[0] = 1.0;
- coords[1] = 0.5;
- coords[2] = 0.5;
- break;
- case 10:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 0.5;
- break;
- case 9:
- coords[0] = 1.0;
- coords[1] = 0.5;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.5*gc[1]*(1.0 - gc[0])*(2.0*gc[1] - 2.0 - gc[0]);
- funValue[2] = 0.5*gc[2]*(1.0 - gc[0])*(2.0*gc[2] - 2.0 - gc[0]);
- funValue[1] = 0.5*(gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(gc[0] + 2.0*gc[1] + 2.0*gc[2]);
-
- funValue[3] = 0.5*gc[1]*(1.0 + gc[0])*(2.0*gc[1] - 2.0 + gc[0]);
- funValue[5] = 0.5*gc[2]*(1.0 + gc[0])*(2.0*gc[2] - 2.0 + gc[0]);
- funValue[4] = 0.5*(-gc[0] - 1.0)*(1.0 - gc[1] - gc[2])*(-gc[0] + 2.0*gc[1] + 2.0*gc[2]);
-
- funValue[8] = 2.0*gc[1]*gc[2]*(1.0 - gc[0]);
- funValue[7] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
- funValue[6] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 - gc[0]);
-
- funValue[12] = gc[1]*(1.0 - gc[0]*gc[0]);
- funValue[14] = gc[2]*(1.0 - gc[0]*gc[0]);
- funValue[13] = (1.0 - gc[1] - gc[2])*(1.0 - gc[0]*gc[0]);
-
- funValue[11] = 2.0*gc[1]*gc[2]*(1.0 + gc[0]);
- funValue[10] = 2.0*gc[2]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- funValue[9] = 2.0*gc[1]*(1.0 - gc[1] - gc[2])*(1.0 + gc[0]);
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * Init Hehahedron Reference coordinates and Shape function.
- * Case A.
- */
-void GaussInfo::hexa8aInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 6:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[1] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
- funValue[3] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
-
- funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- funValue[7] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * Init Hehahedron Reference coordinates and Shape function.
- * Case B.
- */
-void GaussInfo::hexa8bInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 7:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 6:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[3] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
- funValue[1] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
-
- funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[7] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- funValue[5] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * This shapefunc map is same as degenerated quad4bInit
- */
-void GaussInfo::hexa8DegQuad4aInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0]);
- funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0]);
- funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0]);
- funValue[3] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[4] = 0.;
- funValue[5] = 0.;
- funValue[6] = 0.;
- funValue[7] = 0.;
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * This shapefunc map is same as degenerated quad4bInit
- */
-void GaussInfo::hexa8DegQuad4bInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 4:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 5:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 6:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- case 7:
- coords[0] = 0.0;
- coords[1] = 0.0;
- coords[2] = 0.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
- funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
- funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
- funValue[4] = 0.;
- funValue[5] = 0.;
- funValue[6] = 0.;
- funValue[7] = 0.;
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * This shapefunc map is same as degenerated quad4cInit
- */
-void GaussInfo::hexa8DegQuad4cInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = -1.0;
+ coords[0] = 1.0;
+ coords[1] = 1.0;
coords[2] = 0.0;
break;
case 4:
coords[1] = 0.0;
coords[2] = 0.0;
break;
-
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
- funValue[1] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
- funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[3] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
- funValue[4] = 0. ;
- funValue[5] = 0. ;
- funValue[6] = 0. ;
- funValue[7] = 0. ;
- SHAPE_FUN_MACRO_END;
-}
-
-/*!
- * Init Qaudratic Hehahedron Reference coordinates and Shape function.
- * Case A.
- */
-void GaussInfo::hexa20aInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 1:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 3:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 6:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 7:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
-
- case 8:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 9:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 10:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 11:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 12:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 13:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 14:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 15:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 16:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 17:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 18:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 19:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
- funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
- (-2.0 - gc[0] - gc[1] - gc[2]);
- funValue[1] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
- (-2.0 + gc[0] - gc[1] - gc[2]);
- funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
- (-2.0 + gc[0] + gc[1] - gc[2]);
- funValue[3] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
- (-2.0 - gc[0] + gc[1] - gc[2]);
- funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
- (-2.0 - gc[0] - gc[1] + gc[2]);
- funValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
- (-2.0 + gc[0] - gc[1] + gc[2]);
- funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
- (-2.0 + gc[0] + gc[1] + gc[2]);
- funValue[7] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
- (-2.0 - gc[0] + gc[1] + gc[2]);
-
- funValue[8] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[9] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 - gc[2]);
- funValue[10] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
- funValue[11] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 - gc[2]);
- funValue[12] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 - gc[1]);
- funValue[13] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 - gc[1]);
- funValue[14] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[15] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 + gc[1]);
- funValue[16] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[17] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 + gc[2]);
- funValue[18] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- funValue[19] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 + gc[2]);
- SHAPE_FUN_MACRO_END;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 + gc[1])*(1.0 - gc[0]);
+ funValue[1] = 0.25*(1.0 - gc[1])*(1.0 - gc[0]);
+ funValue[2] = 0.25*(1.0 - gc[1])*(1.0 + gc[0]);
+ funValue[3] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[4] = 0.;
+ funValue[5] = 0.;
+ funValue[6] = 0.;
+ funValue[7] = 0.;
+ SHAPE_FUN_MACRO_END;
}
/*!
- * Init Qaudratic Hehahedron Reference coordinates and Shape function.
- * Case B.
- */
-void GaussInfo::hexa20bInit()
-{
- LOCAL_COORD_MACRO_BEGIN;
- case 0:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 3:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 7:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 6:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
-
- case 11:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 10:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 9:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 8:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 16:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 19:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 18:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 17:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 15:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 14:
- coords[0] = 1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- case 13:
- coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 12:
- coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = 1.0;
- break;
- LOCAL_COORD_MACRO_END;
-
- SHAPE_FUN_MACRO_BEGIN;
-
- funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
- (-2.0 - gc[0] - gc[1] - gc[2]);
- funValue[3] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
- (-2.0 + gc[0] - gc[1] - gc[2]);
- funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
- (-2.0 + gc[0] + gc[1] - gc[2]);
- funValue[1] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
- (-2.0 - gc[0] + gc[1] - gc[2]);
- funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
- (-2.0 - gc[0] - gc[1] + gc[2]);
- funValue[7] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
- (-2.0 + gc[0] - gc[1] + gc[2]);
- funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
- (-2.0 + gc[0] + gc[1] + gc[2]);
- funValue[5] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
- (-2.0 - gc[0] + gc[1] + gc[2]);
-
- funValue[11] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
- funValue[10] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 - gc[2]);
- funValue[9] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
- funValue[8] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 - gc[2]);
- funValue[16] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 - gc[1]);
- funValue[19] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 - gc[1]);
- funValue[18] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 + gc[1]);
- funValue[17] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 + gc[1]);
- funValue[15] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
- funValue[14] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 + gc[2]);
- funValue[13] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
- funValue[12] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 + gc[2]);
- SHAPE_FUN_MACRO_END;
-}
-
-void GaussInfo::hexa27aInit()
+ * This shapefunc map is same as degenerated quad4bInit
+ */
+void GaussInfo::hexa8DegQuad4bInit()
{
LOCAL_COORD_MACRO_BEGIN;
- case 0:
+ case 0:
coords[0] = -1.0;
coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 1:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
- break;
- case 2:
- coords[0] = 1.0;
- coords[1] = 1.0;
- coords[2] = -1.0;
+ coords[2] = 0.0;
break;
- case 3:
+ case 1:
coords[0] = 1.0;
coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 4:
- coords[0] = -1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
- break;
- case 5:
- coords[0] = -1.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
+ coords[2] = 0.0;
break;
- case 6:
+ case 2:
coords[0] = 1.0;
coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 7:
- coords[0] = 1.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
+ coords[2] = 0.0;
break;
- case 8:
+ case 3:
coords[0] = -1.0;
- coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 9:
- coords[0] = 0.0;
coords[1] = 1.0;
- coords[2] = -1.0;
+ coords[2] = 0.0;
break;
- case 10:
- coords[0] = 1.0;
+ case 4:
+ coords[0] = 0.0;
coords[1] = 0.0;
- coords[2] = -1.0;
+ coords[2] = 0.0;
break;
- case 11:
+ case 5:
coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = -1.0;
- break;
- case 12:
- coords[0] = -1.0;
coords[1] = 0.0;
- coords[2] = 1.0;
+ coords[2] = 0.0;
break;
- case 13:
+ case 6:
coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 1.0;
- break;
- case 14:
- coords[0] = 1.0;
coords[1] = 0.0;
- coords[2] = 1.0;
+ coords[2] = 0.0;
break;
- case 15:
+ case 7:
coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 1.0;
+ coords[1] = 0.0;
+ coords[2] = 0.0;
break;
- case 16:
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
+ funValue[1] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
+ funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[3] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
+ funValue[4] = 0.;
+ funValue[5] = 0.;
+ funValue[6] = 0.;
+ funValue[7] = 0.;
+ SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * This shapefunc map is same as degenerated quad4cInit
+ */
+void GaussInfo::hexa8DegQuad4cInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
coords[0] = -1.0;
coords[1] = -1.0;
coords[2] = 0.0;
break;
- case 17:
+ case 1:
coords[0] = -1.0;
coords[1] = 1.0;
coords[2] = 0.0;
break;
- case 18:
+ case 2:
coords[0] = 1.0;
coords[1] = 1.0;
coords[2] = 0.0;
break;
- case 19:
+ case 3:
coords[0] = 1.0;
coords[1] = -1.0;
coords[2] = 0.0;
break;
- case 20:
+ case 4:
coords[0] = 0.0;
coords[1] = 0.0;
- coords[2] = -1.0;
- break;
- case 21:
- coords[0] = -1.0;
- coords[1] = 0.0;
coords[2] = 0.0;
break;
- case 22:
+ case 5:
coords[0] = 0.0;
- coords[1] = 1.0;
- coords[2] = 0.0;
- break;
- case 23:
- coords[0] = 1.0;
coords[1] = 0.0;
coords[2] = 0.0;
break;
- case 24:
- coords[0] = 0.0;
- coords[1] = -1.0;
- coords[2] = 0.0;
- break;
- case 25:
+ case 6:
coords[0] = 0.0;
coords[1] = 0.0;
- coords[2] = 1.0;
+ coords[2] = 0.0;
break;
- case 26:
+ case 7:
coords[0] = 0.0;
coords[1] = 0.0;
coords[2] = 0.0;
break;
+
LOCAL_COORD_MACRO_END;
SHAPE_FUN_MACRO_BEGIN;
-
- funValue[0] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
- funValue[1] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
- funValue[2] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
- funValue[3] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
- funValue[4] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
- funValue[5] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
- funValue[6] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
- funValue[7] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
- funValue[8] =0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
- funValue[9] =0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
- funValue[10]=0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
- funValue[11]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
- funValue[12]=0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
- funValue[13]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
- funValue[14]=0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
- funValue[15]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
- funValue[16]=0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
- funValue[17]=0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
- funValue[18]=0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
- funValue[19]=0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
- funValue[20]=0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
- funValue[21]=0.5*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
- funValue[22]=0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
- funValue[23]=0.5*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
- funValue[24]=0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
- funValue[25]=0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
- funValue[26]=(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
-
+ funValue[0] = 0.25*(1.0 - gc[0])*(1.0 - gc[1]);
+ funValue[1] = 0.25*(1.0 - gc[0])*(1.0 + gc[1]);
+ funValue[2] = 0.25*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[3] = 0.25*(1.0 + gc[0])*(1.0 - gc[1]);
+ funValue[4] = 0. ;
+ funValue[5] = 0. ;
+ funValue[6] = 0. ;
+ funValue[7] = 0. ;
SHAPE_FUN_MACRO_END;
}
+/*!
+ * Init Qaudratic Hehahedron Reference coordinates and Shape function.
+ * Case A.
+ */
+void GaussInfo::hexa20aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = HEXA20A_REF[0];
+ coords[1] = HEXA20A_REF[1];
+ coords[2] = HEXA20A_REF[2];
+ break;
+ case 1:
+ coords[0] = HEXA20A_REF[3];
+ coords[1] = HEXA20A_REF[4];
+ coords[2] = HEXA20A_REF[5];
+ break;
+ case 2:
+ coords[0] = HEXA20A_REF[6];
+ coords[1] = HEXA20A_REF[7];
+ coords[2] = HEXA20A_REF[8];
+ break;
+ case 3:
+ coords[0] = HEXA20A_REF[9];
+ coords[1] = HEXA20A_REF[10];
+ coords[2] = HEXA20A_REF[11];
+ break;
+ case 4:
+ coords[0] = HEXA20A_REF[12];
+ coords[1] = HEXA20A_REF[13];
+ coords[2] = HEXA20A_REF[14];
+ break;
+ case 5:
+ coords[0] = HEXA20A_REF[15];
+ coords[1] = HEXA20A_REF[16];
+ coords[2] = HEXA20A_REF[17];
+ break;
+ case 6:
+ coords[0] = HEXA20A_REF[18];
+ coords[1] = HEXA20A_REF[19];
+ coords[2] = HEXA20A_REF[20];
+ break;
+ case 7:
+ coords[0] = HEXA20A_REF[21];
+ coords[1] = HEXA20A_REF[22];
+ coords[2] = HEXA20A_REF[23];
+ break;
+ case 8:
+ coords[0] = HEXA20A_REF[24];
+ coords[1] = HEXA20A_REF[25];
+ coords[2] = HEXA20A_REF[26];
+ break;
+ case 9:
+ coords[0] = HEXA20A_REF[27];
+ coords[1] = HEXA20A_REF[28];
+ coords[2] = HEXA20A_REF[29];
+ break;
+ case 10:
+ coords[0] = HEXA20A_REF[30];
+ coords[1] = HEXA20A_REF[31];
+ coords[2] = HEXA20A_REF[32];
+ break;
+ case 11:
+ coords[0] = HEXA20A_REF[33];
+ coords[1] = HEXA20A_REF[34];
+ coords[2] = HEXA20A_REF[35];
+ break;
+ case 12:
+ coords[0] = HEXA20A_REF[36];
+ coords[1] = HEXA20A_REF[37];
+ coords[2] = HEXA20A_REF[38];
+ break;
+ case 13:
+ coords[0] = HEXA20A_REF[39];
+ coords[1] = HEXA20A_REF[40];
+ coords[2] = HEXA20A_REF[41];
+ break;
+ case 14:
+ coords[0] = HEXA20A_REF[42];
+ coords[1] = HEXA20A_REF[43];
+ coords[2] = HEXA20A_REF[44];
+ break;
+ case 15:
+ coords[0] = HEXA20A_REF[45];
+ coords[1] = HEXA20A_REF[46];
+ coords[2] = HEXA20A_REF[47];
+ break;
+ case 16:
+ coords[0] = HEXA20A_REF[48];
+ coords[1] = HEXA20A_REF[49];
+ coords[2] = HEXA20A_REF[50];
+ break;
+ case 17:
+ coords[0] = HEXA20A_REF[51];
+ coords[1] = HEXA20A_REF[52];
+ coords[2] = HEXA20A_REF[53];
+ break;
+ case 18:
+ coords[0] = HEXA20A_REF[54];
+ coords[1] = HEXA20A_REF[55];
+ coords[2] = HEXA20A_REF[56];
+ break;
+ case 19:
+ coords[0] = HEXA20A_REF[57];
+ coords[1] = HEXA20A_REF[58];
+ coords[2] = HEXA20A_REF[59];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+ funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
+ (-2.0 - gc[0] - gc[1] - gc[2]);
+ funValue[1] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
+ (-2.0 + gc[0] - gc[1] - gc[2]);
+ funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
+ (-2.0 + gc[0] + gc[1] - gc[2]);
+ funValue[3] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
+ (-2.0 - gc[0] + gc[1] - gc[2]);
+ funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
+ (-2.0 - gc[0] - gc[1] + gc[2]);
+ funValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
+ (-2.0 + gc[0] - gc[1] + gc[2]);
+ funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
+ (-2.0 + gc[0] + gc[1] + gc[2]);
+ funValue[7] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
+ (-2.0 - gc[0] + gc[1] + gc[2]);
+
+ funValue[8] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[9] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 - gc[2]);
+ funValue[10] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ funValue[11] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 - gc[2]);
+ funValue[12] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 - gc[1]);
+ funValue[13] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 - gc[1]);
+ funValue[14] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[15] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 + gc[1]);
+ funValue[16] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[17] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 + gc[2]);
+ funValue[18] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ funValue[19] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 + gc[2]);
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.125*(1. + gc[1] + gc[2] + 2* gc[0]) * (1.0 - gc[1])*(1.0 - gc[2]);
+ devFunValue[1] = 0.125*(1.0 - gc[0])*(1. + gc[0] + gc[2] + 2 * gc[1])*(1.0 - gc[2]);
+ devFunValue[2] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[0] + gc[1] + 2 *gc[2]);
+
+ devFunValue[3] = 0.125*(-1.0 - gc[1] - gc[2] + 2*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ devFunValue[4] = 0.125*(1.0 + gc[0])* (1 - gc[0] + gc[2] + 2*gc[1]) *(1.0 - gc[2]);
+ devFunValue[5] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])* (1 - gc[0] + gc[1] + 2*gc[2]);
+
+ devFunValue[6] = 0.125*(-1.0 + gc[1] - gc[2] + 2*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ devFunValue[7] = 0.125*(1.0 + gc[0])* (-1 + gc[0] - gc[2] + 2*gc[1]) *(1.0 - gc[2]);
+ devFunValue[8] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[0] - gc[1] + 2*gc[2]);
+
+ devFunValue[9] = 0.125*(1.0 - gc[1] + gc[2] + 2*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ devFunValue[10] = 0.125*(1.0 - gc[0])*(-1 - gc[0] - gc[2] + 2*gc[1])*(1.0 - gc[2]);
+ devFunValue[11] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[0] - gc[1] + 2*gc[2]);
+
+ devFunValue[12] = 0.125*(1.0 + gc[1] - gc[2] + 2*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ devFunValue[13] = 0.125*(1.0 - gc[0])*(1 + gc[0] - gc[2] + 2*gc[1])*(1.0 + gc[2]);
+ devFunValue[14] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(-1.0 - gc[0] - gc[1] + 2*gc[2]);
+
+ devFunValue[15] = 0.125*(-1.0 - gc[1] + gc[2] + 2*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ devFunValue[16] = 0.125*(1.0 + gc[0])*(1 - gc[0] - gc[2] + 2*gc[1])*(1.0 + gc[2]);
+ devFunValue[17] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(-1.0 + gc[0] - gc[1] + 2*gc[2]);
+
+ devFunValue[18] = 0.125*(-1.0 + gc[1] + gc[2] + 2*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ devFunValue[19] = 0.125*(1.0 + gc[0])*(-1 + gc[0] + gc[2] + 2*gc[1])*(1.0 + gc[2]);
+ devFunValue[20] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(-1.0 + gc[0] + gc[1] + 2*gc[2]);
+
+ devFunValue[21] = 0.125*(1.0 - gc[1] - gc[2] + 2*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ devFunValue[22] = 0.125*(1.0 - gc[0])*(-1 - gc[0] + gc[2] + 2*gc[1])*(1.0 + gc[2]);
+ devFunValue[23] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(-1.0 - gc[0] + gc[1] + 2*gc[2]);
+
+ devFunValue[24] = 0.25*(-2.0*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ devFunValue[25] = 0.25*(1.0 - gc[0]*gc[0])*(-1)*(1.0 - gc[2]);
+ devFunValue[26] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(-1.0);
+
+ devFunValue[27] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[2]);
+ devFunValue[28] = 0.25*(-2.0*gc[1])*(1.0 + gc[0])*(1.0 - gc[2]);
+ devFunValue[29] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(-1.0);
+
+ devFunValue[30] = 0.25*(-2.0*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ devFunValue[31] = 0.25*(1.0 - gc[0]*gc[0])*(1.0)*(1.0 - gc[2]);
+ devFunValue[32] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(-1.0);
+
+ devFunValue[33] = 0.25*(1.0 - gc[1]*gc[1])*(-1.0)*(1.0 - gc[2]);
+ devFunValue[34] = 0.25*(-2.0*gc[1])*(1.0 - gc[0])*(1.0 - gc[2]);
+ devFunValue[35] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(-1.0);
+
+ devFunValue[36] = 0.25*(1.0 - gc[2]*gc[2])*(-1.0)*(1.0 - gc[1]);
+ devFunValue[37] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(-1.0);
+ devFunValue[38] = 0.25*(-2.0*gc[2])*(1.0 - gc[0])*(1.0 - gc[1]);
+
+ devFunValue[39] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[1]);
+ devFunValue[40] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(-1.0);
+ devFunValue[41] = 0.25*(-2.0*gc[2])*(1.0 + gc[0])*(1.0 - gc[1]);
+
+ devFunValue[42] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[1]);
+ devFunValue[43] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0]);
+ devFunValue[44] = 0.25*(-2.0*gc[2])*(1.0 + gc[0])*(1.0 + gc[1]);
+
+ devFunValue[45] = 0.25*(1.0 - gc[2]*gc[2])*(-1.0)*(1.0 + gc[1]);
+ devFunValue[46] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0]);
+ devFunValue[47] = 0.25*(-2.0*gc[2])*(1.0 - gc[0])*(1.0 + gc[1]);
+
+ devFunValue[48] = 0.25*(-2.0*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ devFunValue[49] = 0.25*(1.0 - gc[0]*gc[0])*(-1.0)*(1.0 + gc[2]);
+ devFunValue[50] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1]);
+
+ devFunValue[51] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[2]);
+ devFunValue[52] = 0.25*(-2.0*gc[1])*(1.0 + gc[0])*(1.0 + gc[2]);
+ devFunValue[53] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0]);
+
+ devFunValue[54] = 0.25*(-2.0*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ devFunValue[55] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[2]);
+ devFunValue[56] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1]);
+
+ devFunValue[57] = 0.25*(1.0 - gc[1]*gc[1])*(-1.0)*(1.0 + gc[2]);
+ devFunValue[58] = 0.25*(-2.0*gc[1])*(1.0 - gc[0])*(1.0 + gc[2]);
+ devFunValue[59] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0]);
+
+ SHAPE_FUN_MACRO_END;
+}
+
+/*!
+ * Init Qaudratic Hehahedron Reference coordinates and Shape function.
+ * Case B.
+ */
+void GaussInfo::hexa20bInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = HEXA20B_REF[0];
+ coords[1] = HEXA20B_REF[1];
+ coords[2] = HEXA20B_REF[2];
+ break;
+ case 1:
+ coords[0] = HEXA20B_REF[3];
+ coords[1] = HEXA20B_REF[4];
+ coords[2] = HEXA20B_REF[5];
+ break;
+ case 2:
+ coords[0] = HEXA20B_REF[6];
+ coords[1] = HEXA20B_REF[7];
+ coords[2] = HEXA20B_REF[8];
+ break;
+ case 3:
+ coords[0] = HEXA20B_REF[9];
+ coords[1] = HEXA20B_REF[10];
+ coords[2] = HEXA20B_REF[11];
+ break;
+ case 4:
+ coords[0] = HEXA20B_REF[12];
+ coords[1] = HEXA20B_REF[13];
+ coords[2] = HEXA20B_REF[14];
+ break;
+ case 5:
+ coords[0] = HEXA20B_REF[15];
+ coords[1] = HEXA20B_REF[16];
+ coords[2] = HEXA20B_REF[17];
+ break;
+ case 6:
+ coords[0] = HEXA20B_REF[18];
+ coords[1] = HEXA20B_REF[19];
+ coords[2] = HEXA20B_REF[20];
+ break;
+ case 7:
+ coords[0] = HEXA20B_REF[21];
+ coords[1] = HEXA20B_REF[22];
+ coords[2] = HEXA20B_REF[23];
+ break;
+ case 8:
+ coords[0] = HEXA20B_REF[24];
+ coords[1] = HEXA20B_REF[25];
+ coords[2] = HEXA20B_REF[26];
+ break;
+ case 9:
+ coords[0] = HEXA20B_REF[27];
+ coords[1] = HEXA20B_REF[28];
+ coords[2] = HEXA20B_REF[29];
+ break;
+ case 10:
+ coords[0] = HEXA20B_REF[30];
+ coords[1] = HEXA20B_REF[31];
+ coords[2] = HEXA20B_REF[32];
+ break;
+ case 11:
+ coords[0] = HEXA20B_REF[33];
+ coords[1] = HEXA20B_REF[34];
+ coords[2] = HEXA20B_REF[35];
+ break;
+ case 12:
+ coords[0] = HEXA20B_REF[36];
+ coords[1] = HEXA20B_REF[37];
+ coords[2] = HEXA20B_REF[38];
+ break;
+ case 13:
+ coords[0] = HEXA20B_REF[39];
+ coords[1] = HEXA20B_REF[40];
+ coords[2] = HEXA20B_REF[41];
+ break;
+ case 14:
+ coords[0] = HEXA20B_REF[42];
+ coords[1] = HEXA20B_REF[43];
+ coords[2] = HEXA20B_REF[44];
+ break;
+ case 15:
+ coords[0] = HEXA20B_REF[45];
+ coords[1] = HEXA20B_REF[46];
+ coords[2] = HEXA20B_REF[47];
+ break;
+ case 16:
+ coords[0] = HEXA20B_REF[48];
+ coords[1] = HEXA20B_REF[49];
+ coords[2] = HEXA20B_REF[50];
+ break;
+ case 17:
+ coords[0] = HEXA20B_REF[51];
+ coords[1] = HEXA20B_REF[52];
+ coords[2] = HEXA20B_REF[53];
+ break;
+ case 18:
+ coords[0] = HEXA20B_REF[54];
+ coords[1] = HEXA20B_REF[55];
+ coords[2] = HEXA20B_REF[56];
+ break;
+ case 19:
+ coords[0] = HEXA20B_REF[57];
+ coords[1] = HEXA20B_REF[58];
+ coords[2] = HEXA20B_REF[59];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+
+ funValue[0] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
+ (-2.0 - gc[0] - gc[1] - gc[2]);
+ funValue[3] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 - gc[2])*
+ (-2.0 + gc[0] - gc[1] - gc[2]);
+ funValue[2] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
+ (-2.0 + gc[0] + gc[1] - gc[2]);
+ funValue[1] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 - gc[2])*
+ (-2.0 - gc[0] + gc[1] - gc[2]);
+ funValue[4] = 0.125*(1.0 - gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
+ (-2.0 - gc[0] - gc[1] + gc[2]);
+ funValue[7] = 0.125*(1.0 + gc[0])*(1.0 - gc[1])*(1.0 + gc[2])*
+ (-2.0 + gc[0] - gc[1] + gc[2]);
+ funValue[6] = 0.125*(1.0 + gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
+ (-2.0 + gc[0] + gc[1] + gc[2]);
+ funValue[5] = 0.125*(1.0 - gc[0])*(1.0 + gc[1])*(1.0 + gc[2])*
+ (-2.0 - gc[0] + gc[1] + gc[2]);
+
+ funValue[11] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 - gc[2]);
+ funValue[10] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 - gc[2]);
+ funValue[9] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 - gc[2]);
+ funValue[8] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 - gc[2]);
+ funValue[16] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 - gc[1]);
+ funValue[19] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 - gc[1]);
+ funValue[18] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 + gc[0])*(1.0 + gc[1]);
+ funValue[17] = 0.25*(1.0 - gc[2]*gc[2])*(1.0 - gc[0])*(1.0 + gc[1]);
+ funValue[15] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 - gc[1])*(1.0 + gc[2]);
+ funValue[14] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 + gc[0])*(1.0 + gc[2]);
+ funValue[13] = 0.25*(1.0 - gc[0]*gc[0])*(1.0 + gc[1])*(1.0 + gc[2]);
+ funValue[12] = 0.25*(1.0 - gc[1]*gc[1])*(1.0 - gc[0])*(1.0 + gc[2]);
+ SHAPE_FUN_MACRO_END;
+}
+
+void GaussInfo::hexa27aInit()
+{
+ LOCAL_COORD_MACRO_BEGIN;
+ case 0:
+ coords[0] = HEXA27A_REF[0];
+ coords[1] = HEXA27A_REF[1];
+ coords[2] = HEXA27A_REF[2];
+ break;
+ case 1:
+ coords[0] = HEXA27A_REF[3];
+ coords[1] = HEXA27A_REF[4];
+ coords[2] = HEXA27A_REF[5];
+ break;
+ case 2:
+ coords[0] = HEXA27A_REF[6];
+ coords[1] = HEXA27A_REF[7];
+ coords[2] = HEXA27A_REF[8];
+ break;
+ case 3:
+ coords[0] = HEXA27A_REF[9];
+ coords[1] = HEXA27A_REF[10];
+ coords[2] = HEXA27A_REF[11];
+ break;
+ case 4:
+ coords[0] = HEXA27A_REF[12];
+ coords[1] = HEXA27A_REF[13];
+ coords[2] = HEXA27A_REF[14];
+ break;
+ case 5:
+ coords[0] = HEXA27A_REF[15];
+ coords[1] = HEXA27A_REF[16];
+ coords[2] = HEXA27A_REF[17];
+ break;
+ case 6:
+ coords[0] = HEXA27A_REF[18];
+ coords[1] = HEXA27A_REF[19];
+ coords[2] = HEXA27A_REF[20];
+ break;
+ case 7:
+ coords[0] = HEXA27A_REF[21];
+ coords[1] = HEXA27A_REF[22];
+ coords[2] = HEXA27A_REF[23];
+ break;
+ case 8:
+ coords[0] = HEXA27A_REF[24];
+ coords[1] = HEXA27A_REF[25];
+ coords[2] = HEXA27A_REF[26];
+ break;
+ case 9:
+ coords[0] = HEXA27A_REF[27];
+ coords[1] = HEXA27A_REF[28];
+ coords[2] = HEXA27A_REF[29];
+ break;
+ case 10:
+ coords[0] = HEXA27A_REF[30];
+ coords[1] = HEXA27A_REF[31];
+ coords[2] = HEXA27A_REF[32];
+ break;
+ case 11:
+ coords[0] = HEXA27A_REF[33];
+ coords[1] = HEXA27A_REF[34];
+ coords[2] = HEXA27A_REF[35];
+ break;
+ case 12:
+ coords[0] = HEXA27A_REF[36];
+ coords[1] = HEXA27A_REF[37];
+ coords[2] = HEXA27A_REF[38];
+ break;
+ case 13:
+ coords[0] = HEXA27A_REF[39];
+ coords[1] = HEXA27A_REF[40];
+ coords[2] = HEXA27A_REF[41];
+ break;
+ case 14:
+ coords[0] = HEXA27A_REF[42];
+ coords[1] = HEXA27A_REF[43];
+ coords[2] = HEXA27A_REF[44];
+ break;
+ case 15:
+ coords[0] = HEXA27A_REF[45];
+ coords[1] = HEXA27A_REF[46];
+ coords[2] = HEXA27A_REF[47];
+ break;
+ case 16:
+ coords[0] = HEXA27A_REF[48];
+ coords[1] = HEXA27A_REF[49];
+ coords[2] = HEXA27A_REF[50];
+ break;
+ case 17:
+ coords[0] = HEXA27A_REF[51];
+ coords[1] = HEXA27A_REF[52];
+ coords[2] = HEXA27A_REF[53];
+ break;
+ case 18:
+ coords[0] = HEXA27A_REF[54];
+ coords[1] = HEXA27A_REF[55];
+ coords[2] = HEXA27A_REF[56];
+ break;
+ case 19:
+ coords[0] = HEXA27A_REF[57];
+ coords[1] = HEXA27A_REF[58];
+ coords[2] = HEXA27A_REF[59];
+ break;
+ case 20:
+ coords[0] = HEXA27A_REF[60];
+ coords[1] = HEXA27A_REF[61];
+ coords[2] = HEXA27A_REF[62];
+ break;
+ case 21:
+ coords[0] = HEXA27A_REF[63];
+ coords[1] = HEXA27A_REF[64];
+ coords[2] = HEXA27A_REF[65];
+ break;
+ case 22:
+ coords[0] = HEXA27A_REF[66];
+ coords[1] = HEXA27A_REF[67];
+ coords[2] = HEXA27A_REF[68];
+ break;
+ case 23:
+ coords[0] = HEXA27A_REF[69];
+ coords[1] = HEXA27A_REF[70];
+ coords[2] = HEXA27A_REF[71];
+ break;
+ case 24:
+ coords[0] = HEXA27A_REF[72];
+ coords[1] = HEXA27A_REF[73];
+ coords[2] = HEXA27A_REF[74];
+ break;
+ case 25:
+ coords[0] = HEXA27A_REF[75];
+ coords[1] = HEXA27A_REF[76];
+ coords[2] = HEXA27A_REF[77];
+ break;
+ case 26:
+ coords[0] = HEXA27A_REF[78];
+ coords[1] = HEXA27A_REF[79];
+ coords[2] = HEXA27A_REF[80];
+ break;
+ LOCAL_COORD_MACRO_END;
+
+ SHAPE_FUN_MACRO_BEGIN;
+
+ funValue[0] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ funValue[1] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ funValue[2] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ funValue[3] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ funValue[4] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ funValue[5] =0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ funValue[6] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ funValue[7] =0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ funValue[8] =0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ funValue[9] =0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ funValue[10]=0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ funValue[11]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ funValue[12]=0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ funValue[13]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ funValue[14]=0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ funValue[15]=0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ funValue[16]=0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ funValue[17]=0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ funValue[18]=0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ funValue[19]=0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ funValue[20]=0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ funValue[21]=0.5*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+ funValue[22]=0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ funValue[23]=0.5*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+ funValue[24]=0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ funValue[25]=0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ funValue[26]=(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+
+ SHAPE_FUN_MACRO_END;
+
+ DEV_SHAPE_FUN_MACRO_BEGIN;
+
+ devFunValue[0] = 0.125*(2*gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[1] = 0.125*gc[0]*(gc[0]-1.)*(2*gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[2] = 0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*(2*gc[2]-1.);
+
+ devFunValue[3] = 0.125*(2*gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[4] = 0.125*gc[0]*(gc[0]-1.)*(2*gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[5] = 0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*(2*gc[2]-1.);
+
+ devFunValue[6] = 0.125*(2*gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[7] = 0.125*gc[0]*(gc[0]+1.)*(2*gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[8] = 0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*(2*gc[2]-1.);
+
+ devFunValue[9] = 0.125*(2*gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[10] = 0.125*gc[0]*(gc[0]+1.)*(2*gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[11] = 0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*(2*gc[2]-1.);
+
+ devFunValue[12] = 0.125*(2*gc[0]-1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[13] = 0.125*gc[0]*(gc[0]-1.)*(2*gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[14] = 0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*(2*gc[2]+1.);
+
+ devFunValue[15] = 0.125*(2*gc[0]-1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[16] = 0.125*gc[0]*(gc[0]-1.)*(2*gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[17] = 0.125*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*(2*gc[2]+1.);
+
+ devFunValue[18] = 0.125*(2*gc[0]+1.)*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[19] = 0.125*gc[0]*(gc[0]+1.)*(2*gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[20] = 0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*(2*gc[2]+1.);
+
+ devFunValue[21] = 0.125*(2*gc[0]+1.)*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[22] = 0.125*gc[0]*(gc[0]+1.)*(2*gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[23] = 0.125*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*(2*gc[2]+1.);
+
+ devFunValue[24] = 0.25*(2*gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[25] = 0.25*gc[0]*(gc[0]-1.)*(-2*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[26] = 0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*(2*gc[2]-1.);
+
+ devFunValue[27] = 0.25*(-2*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[28] = 0.25*(1.-gc[0]*gc[0])*(2*gc[1]+1.)*gc[2]*(gc[2]-1.);
+ devFunValue[29] = 0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*(2*gc[2]-1.);
+
+ devFunValue[30] = 0.25*(2*gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[31] = 0.25*gc[0]*(gc[0]+1.)*(-2*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[32] = 0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*(2*gc[2]-1.);
+
+ devFunValue[33] = 0.25*(-2*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[34] = 0.25*(1.-gc[0]*gc[0])*(2*gc[1]-1.)*gc[2]*(gc[2]-1.);
+ devFunValue[35] = 0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*(2*gc[2]-1.);
+
+ devFunValue[36] = 0.25*(2*gc[0]-1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[37] = 0.25*gc[0]*(gc[0]-1.)*(-2*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[38] = 0.25*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*(2*gc[2]+1.);
+
+ devFunValue[39] = 0.25*(-2*gc[0])*gc[1]*(gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[40] = 0.25*(1.-gc[0]*gc[0])*(2*gc[1]+1.)*gc[2]*(gc[2]+1.);
+ devFunValue[41] = 0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*(2*gc[2]+1.);
+
+ devFunValue[42] = 0.25*(2*gc[0]+1.)*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[43] = 0.25*gc[0]*(gc[0]+1.)*(-2*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[44] = 0.25*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*(2*gc[2]+1.);
+
+ devFunValue[45] = 0.25*(-2*gc[0])*gc[1]*(gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[46] = 0.25*(1.-gc[0]*gc[0])*(2*gc[1]-1.)*gc[2]*(gc[2]+1.);
+ devFunValue[47] = 0.25*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*(2*gc[2]+1.);
+
+ devFunValue[48] = 0.25*(2*gc[0]-1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[49] = 0.25*gc[0]*(gc[0]-1.)*(2*gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[50] = 0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]-1.)*(-2*gc[2]);
+
+ devFunValue[51] = 0.25*(2*gc[0]-1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[52] = 0.25*gc[0]*(gc[0]-1.)*(2*gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[53] = 0.25*gc[0]*(gc[0]-1.)*gc[1]*(gc[1]+1.)*(-2*gc[2]);
+
+ devFunValue[54] = 0.25*(2*gc[0]+1.)*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[55] = 0.25*gc[0]*(gc[0]+1.)*(2*gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[56] = 0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]+1.)*(-2*gc[2]);
+
+ devFunValue[57] = 0.25*(2*gc[0]+1.)*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[58] = 0.25*gc[0]*(gc[0]+1.)*(2*gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[59] = 0.25*gc[0]*(gc[0]+1.)*gc[1]*(gc[1]-1.)*(-2*gc[2]);
+
+ devFunValue[60] = 0.5*(-2*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[61] = 0.5*(1.-gc[0]*gc[0])*(-2*gc[1])*gc[2]*(gc[2]-1.);
+ devFunValue[62] = 0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*(2*gc[2]-1.);
+
+ devFunValue[63] = 0.5*(2*gc[0]-1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[64] = 0.5*gc[0]*(gc[0]-1.)*(-2*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[65] = 0.5*gc[0]*(gc[0]-1.)*(1.-gc[1]*gc[1])*(-2*gc[2]);
+
+ devFunValue[66] = 0.5*(-2*gc[0])*gc[1]*(gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[67] = 0.5*(1.-gc[0]*gc[0])*(2*gc[1]+1.)*(1.-gc[2]*gc[2]);
+ devFunValue[68] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]+1.)*(-2*gc[2]);
+
+ devFunValue[69] = 0.5*(2*gc[0]+1.)*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[70] = 0.5*gc[0]*(gc[0]+1.)*(-2*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[71] = 0.5*gc[0]*(gc[0]+1.)*(1.-gc[1]*gc[1])*(-2*gc[2]);
+
+ devFunValue[72] = 0.5*(-2*gc[0])*gc[1]*(gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[73] = 0.5*(1.-gc[0]*gc[0])*(2*gc[1]-1.)*(1.-gc[2]*gc[2]);
+ devFunValue[74] = 0.5*(1.-gc[0]*gc[0])*gc[1]*(gc[1]-1.)*(-2*gc[2]);
+
+ devFunValue[75] = 0.5*(-2*gc[0])*(1.-gc[1]*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[76] = 0.5*(1.-gc[0]*gc[0])*(-2*gc[1])*gc[2]*(gc[2]+1.);
+ devFunValue[77] = 0.5*(1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*(2*gc[2]+1.);
+
+ devFunValue[78] = (-2*gc[0])*(1.-gc[1]*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[79] = (1.-gc[0]*gc[0])*(-2*gc[1])*(1.-gc[2]*gc[2]);
+ devFunValue[80] = (1.-gc[0]*gc[0])*(1.-gc[1]*gc[1])*(-2*gc[2]);
+
+ DEV_SHAPE_FUN_MACRO_END;
+}
+
////////////////////////////////////////////////////////////////////////////////////////////////
// GAUSS COORD CLASS //
////////////////////////////////////////////////////////////////////////////////////////////////
void GaussCoords::addGaussInfo( NormalizedCellType theGeometry,
int coordDim,
const double* theGaussCoord,
- int theNbGauss,
+ mcIdType theNbGauss,
const double* theReferenceCoord,
- int theNbRef)
+ mcIdType theNbRef)
{
GaussInfoVector::iterator it = _my_gauss_info.begin();
for( ; it != _my_gauss_info.end(); it++ )
aReferenceCoord.push_back(theReferenceCoord[i]);
- GaussInfo* info = new GaussInfo( theGeometry, aGaussCoord, theNbGauss, aReferenceCoord, theNbRef);
+ GaussInfo* info = new GaussInfo( theGeometry, aGaussCoord, FromIdType<int>(theNbGauss), aReferenceCoord, FromIdType<int>(theNbRef));
info->initLocalInfo();
//If info with cell type doesn't exist add it
}
else
{
- int index = std::distance(_my_gauss_info.begin(),it);
+ std::size_t index = std::distance(_my_gauss_info.begin(),it);
delete (*it);
_my_gauss_info[index] = info;
}
double* GaussCoords::calculateCoords( NormalizedCellType theGeometry,
const double *theNodeCoords,
const int theSpaceDim,
- const int *theIndex)
+ const mcIdType *theIndex)
{
const GaussInfo *info = getInfoGivenCellType(theGeometry);
int nbCoords = theSpaceDim * info->getNbGauss();
}
-void GaussCoords::calculateCoords( NormalizedCellType theGeometry, const double *theNodeCoords, const int theSpaceDim, const int *theIndex, double *result)
+void GaussCoords::calculateCoords( NormalizedCellType theGeometry, const double *theNodeCoords, const int theSpaceDim, const mcIdType *theIndex, double *result)
{
const GaussInfo *info = getInfoGivenCellType(theGeometry);
calculateCoordsAlg(info,theNodeCoords,theSpaceDim,theIndex,result);
}
-void GaussCoords::calculateCoordsAlg(const GaussInfo *info, const double* theNodeCoords, const int theSpaceDim, const int *theIndex, double *result)
+void GaussCoords::calculateCoordsAlg(const GaussInfo *info, const double* theNodeCoords, const int theSpaceDim, const mcIdType *theIndex, double *result)
{
int aConn = info->getNbRef();
