include/finiteelement/FiniteElement.h

#ifndef _FINITEELEMENT_H_
#define _FINITEELEMENT_H_
// ---------------------------------------------
// company   : 
// user      : 
// date      : 22 Mar 2015 17:22
// file      : FiniteElement.h
// markers   : Full
// language  : C++
// generator : Metamill 6.0
// ---------------------------------------------

//#UBLK-BEG-HEADERH
//==============================================================
// Description : 
//--------------------------------------------------------------
// Author      : Gilberto Luis
// Created     : 17 Aug 2011
//--------------------------------------------------------------
// Change history : 
//   17 Aug 2011 (Gilberto Luis) Initial version generated
//
//==============================================================
//#UBLK-END-HEADERH

/*=================================================================================================
=                                                                                                 =
=       Copyright 2010-2015  Marco Lucio Bittencourt / Fabiano Fernandes Bargos                       =
=                        Gilberto Luis Valente / Jorge Luis Suzuki / Allan Patrick Cordeiro Dias  =
=                                                                                                 =
=       Licensed under the Apache License, Version 2.0 (the "License");                               =
=       you may not use this file except in compliance with the License.                              =
=       You may obtain a copy of the License at                                                       =
=                                                                                                 =
=               http://www.apache.org/licenses/LICENSE-2.0                                                =
=                                                                                                 =
=       Unless required by applicable law or agreed to in writing, software                           =
=       distributed under the License is distributed on an "AS IS" BASIS,                             =
=       WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.                      =
=       See the License for the specific language governing permissions and                           =
=       limitations under the License.                                                                =
=                                                                                                 =
=================================================================================================*/


#include <sys/time.h>
#include <sys/resource.h>
#include <iostream>
#include <cmath>
#include <string>
#include <cstdio>
#include "ds/DS_defs.h"
#include "acdp/acdp.h"
#include "ds/Vector.h"
#include "ds/Matrix.h"
#include "ds/OneIndexTable.h"
#include "ds/TwoIndexTable.h"
#include "enumerations/ElementShape_E.h"
#include "enumerations/Entity_E.h"
#include "enumerations/ProblemType_E.h"
#include "enumerations/Kinematics_E.h"
#include "enumerations/StrainMeasure_E.h"
#include "enumerations/CalculationProcedure_E.h"
#include "enumerations/MeshType_E.h"
#include "enumerations/MappingType_E.h"
#include "enumerations/MassStiffType_E.h"
#include "enumerations/StiffnessType_E.h"
#include "enumerations/FiniteFormulation_E.h"
#include "enumerations/pMesh_E.h"
#include "enumerations/GeometricProperty_E.h"
#include "interpolation/shapefunctions/ShapeFunctions.h"
#include "finiteelement/GradientMeasure.h"
#include "finiteelement/FiniteElementAttributes_S.h"
#include "fegroups/GeometricProperties.h"
#include "fegroups/Mapping.h"
#include "loadsets/LoadIntensity.h"
#include "material/Material.h"
#include "symbolicanalyzer/SymbolicParse.h"


using namespace std;


class FiniteElement
{
//#UBLK-BEG-CLASSDEF mm:079de445-f650-11df-817b-b59fee0c7729
//#UBLK-END-CLASSDEF

 public:

   // METHOD=mm:551b4e59-f652-11df-817b-b59fee0c7729
   FiniteElement();

   // METHOD=mm:e9d066db-f64f-11df-817b-b59fee0c7729
   FiniteElement(FiniteElement& Instance, ElementShape_E ElementShape, ProblemType_E ProblemType, unsigned long MaxSolutionOrder, unsigned long Dim, unsigned long MappingOrder, CalculationProcedure_E CalculationProcedure = STANDARD, MeshType_E MeshType = NON_STRUCTURED, Kinematics_E Kinematics = INFINITE, StrainMeasure_E StrainMeasure = INFINITES, FiniteFormulation_E FiniteFormulation = INFINITESIMAL, pMesh_E pMeshDistribution = pUNIFORM, MappingType_E MappingType = ISOPARAMETRIC, MassStiffType_E MassType = CONSIST, StiffnessType_E StiffType = CONSISTENT, MaterialType_E MatType = HOOKE);

   // METHOD=mm:b853c175-0239-11e0-9699-fe1341938431
   FiniteElement(FiniteElement& Instance);

   // METHOD=mm:ae31cfbb-0239-11e0-9699-fe1341938431
   ~FiniteElement();

   // METHOD=mm:acd5a08b-557f-11e0-8a13-cd08645cc9e7
   FiniteElement& operator=(FiniteElement& Instance);

   // METHOD=mm:f2469a89-5580-11e0-8a13-cd08645cc9e7
   operator GradientMeasure&();

   // METHOD=mm:a693f451-5580-11e0-8a13-cd08645cc9e7
   operator ShapeFunctions&();

   // METHOD=mm:5f864c0f-02d6-11e0-89cb-b0660132333d
   void Print(FILE* File);

   // METHOD=mm:2681172b-856e-11e0-9373-ac70ba4510db
   FiniteElementAttributes_S& GetFiniteElementAttributes();

   // METHOD=mm:831bffd0-1fad-4380-8620-35e3841c76ba
   void SetFiniteElementAttributes(FiniteElementAttributes_S& FiniteElementAttributes);

   // METHOD=mm:936c0d78-84da-42eb-9b8c-256820d4e73a
   void SetFiniteElementAttributes(ElementShape_E ElementShape, ProblemType_E ProblemType, unsigned long MaxSolutionOrder, unsigned long Dim, CalculationProcedure_E CalculationProcedure, MeshType_E MeshType, Kinematics_E Kinematics, StrainMeasure_E StrainMeasure, FiniteFormulation_E FiniteFormulation, pMesh_E pMeshDistribution, MappingType_E MappingType, unsigned long MappingOrder, MassStiffType_E MassType, StiffnessType_E StiffType, MaterialType_E MatType);

   // METHOD=mm:ebe921ba-6d4f-4bf8-864c-5f2dedf93162
   unsigned long GetMaxNumberElemDOFs();

   // METHOD=mm:9fa48dd2-ac37-42c1-a51a-e917a0d2d35a
   unsigned long GetMaxNumberElemPDOFs();

   // METHOD=mm:e9d066dc-f64f-11df-817b-b59fee0c7729
   void GetMassMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinates, Mapping& ElementMapping, GeometricProperties& GeoProp, SymmetricMatrix& MassMatrix, double& ElementMeasure, double Density = 1);

   // METHOD=mm:87952c8e-77a3-4165-89b6-7ee3e103efa3
   void GetLocalMassMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinates, GeometricProperties& GeoProp, SymmetricMatrix& MassMatrix, double Density = 1);

   // METHOD=mm:e9d066df-f64f-11df-817b-b59fee0c7729
   void GetLinearStiffnessMatrix(unsigned long ElementNumber, unsigned long ElementOrder, Mapping& ElementMapping, GeometricProperties& GeoProp, Material& Mater, unsigned long NumElemCoordinates, double* ElementCoordinates, SymmetricMatrix& StiffnessMatrix, double* Fe = NULL);

   // METHOD=mm:88e096aa-83da-4a6a-9a05-59e5a3f6abb1
   void GetTangentStiffnessMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce);

   // METHOD=mm:db653021-34a9-47b9-bf7d-caf5ed011cfa
   void GetDTangentStiffnessMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce, Vector& Damage);

   // METHOD=mm:3645245f-f26d-4425-8c13-71165655399f
   void GetQITangentStiffnessMatrix(unsigned long ElementNumber, unsigned long DisplacementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce);

   // METHOD=mm:f4680acc-b012-433a-a7ee-30c7d08b0924
   void GetDQITangentStiffnessMatrix(unsigned long ElementNumber, unsigned long DisplacementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce, Vector& Damage);

   // METHOD=mm:4a5abc90-1f76-472c-af5d-219d1b19414a
   void GetPbar(unsigned long NumIntPoints, double kappa, double* Jacobi, double* detJi, double* Jacobc, Vector& pbar);

   // METHOD=mm:416ac79d-9cd1-4758-b63f-75a21926b1d0
   void GetPproj(unsigned long ElemNum, double kappa, unsigned long DisplacementOrder, unsigned long PressureOrder, unsigned long NumElemCoords, double* ElemCoords, double* ElemCoords1, double* ShapeFunctions, unsigned long NumberSFU, unsigned long NumberSFP, double* Weights, unsigned long NumIntPoints, double* Jacobi, double* detJi, double* Jacobc);

   // METHOD=mm:e9d066de-f64f-11df-817b-b59fee0c7729
   void GetInternalForceVector(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, GeometricProperties& GeoProp, Material& Mater, double* Ue, double* Fe);

   // METHOD=mm:e9d066e3-f64f-11df-817b-b59fee0c7729
   void GetEdgeLoad(unsigned long ElementNumber, unsigned long ElementOrder, LoadIntensity& Load, Mapping& ElementMapping, unsigned long NumElemCoords, double* ElementCoordinates, unsigned long LocalEdge, double* Fe);

   // METHOD=mm:4c8346c3-f66b-11df-8531-ea83aef40aa3
   void GetFaceLoad(unsigned long ElementNumber, unsigned long ElementOrder, LoadIntensity& Load, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, unsigned long LocalFace, double* Fe);

   // METHOD=mm:e9d066e2-f64f-11df-817b-b59fee0c7729
   void GetBodyLoadVector(unsigned long ElementNumber, unsigned long ElementOrder, LoadIntensity& Load, unsigned long NumElemCoordinates, double* ElementCoordinates, Mapping& ElementMapping, double* Fe, double t = 0);

   // METHOD=mm:768818e7-f66a-11df-8531-ea83aef40aa3
   void GetEdgePressure(unsigned long ElementNumber, unsigned long ElementOrder, LoadIntensity& Load, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, double* Fe);

   // METHOD=mm:cf0be7b7-f66b-11df-8531-ea83aef40aa3
   void GetFacePressure(unsigned long ElementNumber, unsigned long ElementOrder, LoadIntensity& Load, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, double* Fe);

   // METHOD=mm:05953d1c-f66b-11df-8531-ea83aef40aa3
   void GetEdgePressureStiffnessMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumberEdgeNodes, LoadIntensity& Load, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, SymmetricMatrix& Ke, double* Fe);

   // METHOD=mm:04184f95-f66c-11df-8531-ea83aef40aa3
   void GetFacePressureStiffnessMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumberFaceNodes, LoadIntensity& Load, unsigned long NumElemCoords, double* ElementCoordinates, Mapping& ElementMapping, SymmetricMatrix& Ke, double* Fe);

   // METHOD=mm:6e8509ca-cce2-11e0-8c35-faf3955cbabc
   void Run1DLocalMatrices(BuiltInArray<unsigned long>& ElementOrder);

   // METHOD=mm:e9d0677a-f64f-11df-817b-b59fee0c7729
   double* Get1DLocalMassMatrix(unsigned long ElementOrder, unsigned long IntegOrder, unsigned long& TotalNumberSF, unsigned long& MatrixOrder, unsigned long& NumberIP);

   // METHOD=mm:e9d0677b-f64f-11df-817b-b59fee0c7729
   double* Get1DLocalStiffnessMatrix(unsigned long ElementOrder, unsigned long IntegOrder, unsigned long& TotalNumberSF, unsigned long& MatrixOrder, unsigned long& NumberIP);

   // METHOD=mm:e9d0677c-f64f-11df-817b-b59fee0c7729
   double* Get1DLocalMixedMatrix(unsigned long ElementOrder, unsigned long IntegOrder, unsigned long& TotalNumberSF, unsigned long& MatrixOrder, unsigned long& NumberIP);

   // METHOD=mm:1fe72d7d-8882-4e38-a8f3-ac1ee896dee4
   double* Interpolation(Vector& UeSol, unsigned long CoefOrder, unsigned long PosOrder, unsigned long& UePosSize);


 protected:

   // METHOD=mm:101ea912-60e4-4451-9db0-fab596c9ac36
   void Run1DLocalMatrices_copy();


   // attr=mm:727c798d-8567-11e0-9373-ac70ba4510db
   FiniteElementAttributes_S FEAtr;

   // attr=mm:4aac6a99-f651-11df-817b-b59fee0c7729
   TwoIndexTable<double> Local1DMassMatrix;

   // attr=mm:5b5cee9f-f651-11df-817b-b59fee0c7729
   TwoIndexTable<double> Local1DStiffMatrix;

   // attr=mm:85aeafcb-f651-11df-817b-b59fee0c7729
   TwoIndexTable<double> Local1DMixedMatrix;

   // attr=mm:2309d179-f652-11df-817b-b59fee0c7729
   GradientMeasure GM;

   // attr=mm:5a2f452f-3b9f-11e0-8f62-9cadaf298b92
   ShapeFunctions SF;

   // attr=mm:9926e099-43f2-40d1-8504-0146b3e38115
   Matrix F;

   // attr=mm:b0ee1b95-5d42-4bfb-8a80-cf123d2403e1
   Matrix gradU;

   // attr=mm:47b6f2a9-9df0-41ab-98c8-96d53a371c16
   SymmetricMatrix S;

   // attr=mm:fce116df-c117-482a-ab06-fb492deea131
   Vector Sv;

   // attr=mm:8c1e525c-c4b4-42c7-b3ce-eb2088520f33
   SymmetricMatrix E;

   // attr=mm:750f2110-2512-47cd-8aaf-6892e795941c
   Vector Ev;

   // attr=mm:3768f20d-4fe0-4f02-be6b-820144c73492
   SymmetricMatrix DMatrix;

   // attr=mm:8eacdd73-8e1d-11e0-9278-b12ff911d20a
   Matrix BMatrix;

   // attr=mm:9ded21d3-119e-4a9f-893c-378260c65b0e
   Matrix NLBMatrix;

   // attr=mm:d2d553b8-6aa6-4584-b397-c0b75f470d3c
   Matrix ISBMatrix;

   // attr=mm:3976ccd9-d789-405e-9cbe-5c98071b33c6
   Vector dJdE;

   // attr=mm:a0bf27a4-46f1-47f6-9039-a93fdc1b7450
   Vector GVector;

   // attr=mm:eb86eaae-abd1-442f-b22e-47ceb3c12ecf
   Vector AuxSF;

   // attr=mm:95de9b8e-d111-406c-8b94-09f630b1b078
   Vector DispP;

   // attr=mm:014441ae-e8f7-45f6-abb1-27e7da050564
   Vector ProjP;

   // attr=mm:0dc8dfa5-b322-4f6c-ac1a-a2c1e1b8a598
   Vector Pcoefs;

   // attr=mm:a027aecf-50c7-4890-99f1-9a66e34dbb86
   Vector Auxv;

   // attr=mm:8f456f11-8e1d-11e0-9278-b12ff911d20a
   SymmetricMatrix Ke;

   // attr=mm:caaa3372-709a-444f-af09-a7954fdbfee0
   SymmetricMatrix Kteaux;

   // attr=mm:92d16e69-bca3-11e0-8f79-b28c8a2c6811
   SymmetricMatrix Me;

   // attr=mm:1437a1c7-a677-49dc-92de-bca516d86f5b
   SymmetricMatrix Meaux;

   // attr=mm:37e95907-9369-11e0-9502-bbd3be26df5e
   SymmetricMatrix Kte;

   // attr=mm:0e1a1fd3-5758-4524-9c06-148f4017bd48
   Matrix Kp;

   // attr=mm:7d5aa117-4095-4003-abcb-342d819f5984
   Matrix Kpaux;

   // attr=mm:291434b7-ddd2-466d-a572-73a4ead2f461
   Matrix Kpaux1;

   // attr=mm:ce227a25-6cf5-4f0d-9e60-de3a918a2e6c
   Vector Re;

   // attr=mm:cd4d0ab8-433b-4cfd-91ec-5b2330deacfb
   Vector Rhs;

   // attr=mm:4adbf07f-dc31-4788-9041-0a6edb848531
   Vector UePos;

   // attr=mm:c391490d-cdb4-11e0-83dd-d9f41665dfa4
   OneIndexTable<double> Mpq;

   // attr=mm:8a59eebd-cdb6-11e0-83dd-d9f41665dfa4
   OneIndexTable<double> Kpq;

   // attr=mm:8a59f6b5-cdb6-11e0-83dd-d9f41665dfa4
   OneIndexTable<double> Dpq;

   // attr=mm:972a4c94-ed01-41ff-a6fa-8295f110e802
   BuiltInArray<long> ElementOrderSets;

   // attr=mm:d7cf80fc-1eda-455d-8a15-f176379810f1
   BuiltInArray<unsigned long> SFIncid;


 private:

//#UBLK-BEG-CLASSMEMB mm:079de445-f650-11df-817b-b59fee0c7729
//#UBLK-END-CLASSMEMB
};

//#UBLK-BEG-GLOBALH
//*** Removed methods: 18 Aug 2011 11:44
//---
//uuid=mm:b6b0ff63-8c7c-11e0-83b2-eb81abfc0587
//   operator NumericalIntegration&();

//*** Removed methods: 23 Aug 2011 15:19
//---
//uuid=mm:e9d0677e-f64f-11df-817b-b59fee0c7729
//   void Read(FILE* DEFFile, ElementShape_E ElementShape, ProblemType_E ProblemType, unsigned short NumberDOFsNode, unsigned short MaxSolutionOrder, CalculationProcedure_E CalculationProcedure = STANDARD, MeshType_E MeshType = pUNIFORM, Kinematics_E Kinematics = INFINITESIMAL, StrainMeasure_E StrainMeasure = INFINITESIMAL, FiniteFormulation_E FiniteFormulation = INFINITESIMAL, MassStiffType_E MassType = CONSISTENT, StiffnessType_E StiffnessType = CONSISTENT);

//*** Removed methods: 30 Nov 2012 15:36
//---
//uuid=mm:2681218f-856e-11e0-9373-ac70ba4510db
//   void SetFiniteElementAttributes(FiniteElementAttributes_S& FiniteElementAttributes);
//---
//uuid=mm:f422f52f-856c-11e0-9373-ac70ba4510db
//   void SetFiniteElementAttributes(ElementShape_E ElementShape, ProblemType_E ProblemType, unsigned long MaxSolutionOrder, unsigned long Dim, CalculationProcedure_E CalculationProcedure, MeshType_E MeshType, Kinematics_E Kinematics, StrainMeasure_E StrainMeasure, FiniteFormulation_E FiniteFormulation, pMesh_E pMeshDistribution, MassStiffType_E MassType, StiffnessType_E StiffType);
//---
//uuid=mm:dbb9156b-3df5-4002-85eb-c7ebc322686a
//   void GetLocalMassMatrix(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinates, GeometricProperties& GeoProp, SymmetricMatrix& MassMatrix, double Density = 1);

//*** Removed methods: 13 May 2014 16:57
//---
//uuid=mm:f06af471-f667-11df-8531-ea83aef40aa3
//   void GetTangentStiffnessMatrix_copy(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce);
//---
//uuid=mm:b389dcfe-973e-41c6-aa43-564e50956b72
//   void GetQITangentStiffnessMatrix_copy(unsigned long ElementNumber, unsigned long ElementOrder, unsigned long NumElemCoordinates, double* ElementCoordinatesi, double* ElementCoordinatesc, Mapping& ElementMap, GeometricProperties& GeoProp, Material& Mater, SymmetricMatrix& StiffnessMatrix, Vector& InternalForce);

//#UBLK-END-GLOBALH


#endif // _FINITEELEMENT_H_