include/solvercontrol/Solver.h

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

//#UBLK-BEG-HEADERH
//#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 <cstdio>
#include <mpi.h>
#include "solvercontrol/PartitionData.h"
#include "solvercontrol/SolverParam_S.h"
#include "solvercontrol/TheoreticalSolution_S.h"
#include "ds/Vector.h"
#include "ds/Matrix.h"
#include "ds/OneIndexTable.h"
#include "ds/TwoIndexTable.h"
#include "enumerations/PostProcessorProgram_E.h"
#include "fegroups/FEGroups.h"
#include "model/Model.h"
#include "model/Nodes.h"
#include "symbolicanalyzer/SymbolicParse.h"
#include "enumerations/ContactSolver_E.h"
#include "enumerations/ContactModel_E.h"


class Solver
{
//#UBLK-BEG-CLASSDEF mm:8ed17010-020e-11e0-9503-bdd33c3b382f
//#UBLK-END-CLASSDEF

 public:

   // METHOD=mm:75c8b726-020e-11e0-9503-bdd33c3b382f
   Solver();

   // METHOD=mm:80bf84ed-3aaa-11e0-87c1-fb85dddc7efb
   ~Solver();

   // METHOD=mm:f12724b7-66f5-11e0-94a8-e19117dd202e
   Solver(Solver& Instance);

   // METHOD=mm:3f3c1168-66f6-11e0-94a8-e19117dd202e
   Solver& operator=(Solver& Instance);

   // METHOD=mm:a0db204d-61ec-11e0-99a5-df56fb5dd99d
   void SetSolverAttributes(SolverParam_S& SolverParameters);

   // METHOD=mm:337cc73d-61ef-11e0-99a5-df56fb5dd99d
   void SetSolverAttributes(SolverType_E SolvType, MassStiffType_E MassType, StiffnessType_E StiffType, ContactModel_E ContactMod, ContactSolver_E ContactSol, LinearSystemSolutionMethod_E LSMethod, unsigned long LSMaxIt, double LSPrec, SolutionAlgorithm_E SoluAlgorithm, unsigned long SolMaxIt, double SolPrec, unsigned long NumLoadSteps, AnalysisType_E Atype, OperatorType_E OpType, double to, double tf, unsigned long NTS, unsigned long FrictionCoef, double PenaltyParameter, unsigned long ContactInteOrder, TimeIntegrationMeth_E TimeIntMeth, unsigned long OutFreqNR, unsigned long OutFreqTime);

   // METHOD=mm:19dfad9f-645b-11e0-8923-dd061b2ac5dc
   void SetTheoreticalSolutionParameters(TheoreticalSolution_S& TheoSol);

   // METHOD=mm:b679fb37-0219-11e0-8f07-c00dff923bc5
   LinearSystemSolutionMethod_E GetLSMethod();

   // METHOD=mm:a620d0eb-8c78-487f-8e4e-66a6ce9a725c
   ContactSolver_E GetContactIntegrationMethod();

   // METHOD=mm:69fd560c-a7c8-42c1-b9c7-b9c51c5cdb63
   ContactModel_E GetContactFrictionCondition();

   // METHOD=mm:89fdff20-1654-476b-9425-e1c2cc1cbd29
   unsigned long GetFrictionCoef();

   // METHOD=mm:e6291d77-0219-11e0-8f07-c00dff923bc5
   SolverType_E GetSolverType();

   // METHOD=mm:f6d986a7-0219-11e0-8f07-c00dff923bc5
   OperatorType_E GetOperatorType();

   // METHOD=mm:048f0399-021a-11e0-8f07-c00dff923bc5
   unsigned long GetIterativeSolutionMaxIterations();

   // METHOD=mm:1a043feb-021a-11e0-8f07-c00dff923bc5
   double GetIterativeSolutionPrecision();

   // METHOD=mm:298388eb-021a-11e0-8f07-c00dff923bc5
   SolutionAlgorithm_E GetSolutionAlgorithm();

   // METHOD=mm:4c15ae1b-021a-11e0-8f07-c00dff923bc5
   unsigned long GetSolutionMaxIterations();

   // METHOD=mm:5df746b7-021a-11e0-8f07-c00dff923bc5
   double GetNewtonRaphsonSolutionPrecision();

   // METHOD=mm:b6c6df37-021a-11e0-8f07-c00dff923bc5
   unsigned long GetNumberLoadSteps();

   // METHOD=mm:75c8b72a-020e-11e0-9503-bdd33c3b382f
   void PostProcessing(PostProcessorProgram_E PostProgram);

   // METHOD=mm:75c8b727-020e-11e0-9503-bdd33c3b382f
   void Read(FILE* FEMFile, FILE* DEFFile, FILE* LogFile);

   // METHOD=mm:71579f31-0615-4534-9938-dbe72bfb50a0
   virtual void Run() = 0;


   // asse=mm:94b7d0c2-d1a6-4b4a-a51f-5ddab5e469b8
   TheoreticalSolution_S * theTheoreticalSolution_S;

   // asse=mm:279080a8-7ca6-4d3f-b58b-7a1273020db2
   SolverParam_S * theSolverParam_S;

   // asse=mm:9559a12c-9f2e-443a-ab8a-2edca2ac8547
   PartitionData * thePartitionData;


 protected:

   // METHOD=mm:75c8b7c2-020e-11e0-9503-bdd33c3b382f
   void ExportGiDMeshFile(char* FileName, FEGroups& Groups, Nodes& Nds);

   // METHOD=mm:7db8263d-8008-4e1f-8a44-f943eedbb89f
   void ExportGiDResFile(char* FileName, FEGroups& Groups, Nodes& Nds);

   // METHOD=mm:75c8b7c3-020e-11e0-9503-bdd33c3b382f
   void ExportGMSHFile(char* FileName, FEGroups& Groups, Nodes& Nds, double* Solution);

   // METHOD=mm:75c8b7c4-020e-11e0-9503-bdd33c3b382f
   void ExportHyperMeshFile(char* FileName, FEGroups& Groups, Nodes& Nds, double* Solution);

   // METHOD=mm:f5e88b17-3abd-11e0-87c1-fb85dddc7efb
   void ComputeErrorNorm();

   // METHOD=mm:8c04ba7a-6d6b-4c48-b6c7-a098d3a37159
   void ComputeErrorNorm_Parallel(Vector& UPart);

   // METHOD=mm:65dc182c-eba6-4483-9347-67c2652a16c1
   void ComputeErrorNorm1DMatrices();

   // METHOD=mm:ba356a64-7365-4a68-b5f3-9de96b1da240
   void ComputeErrorNorm1DMatrices_Parallel(Vector& UPart);

   // METHOD=mm:7b3dba8d-3df6-11e0-891a-d4865a8b0383
   void MatrixSchurComplement(SymmetricMatrix& A, unsigned long NBD, unsigned long NID);

   // METHOD=mm:ecf096c5-9f5f-48e1-8fc1-78f49e3a24ea
   void MatrixSchurComplement(SymmetricMatrix& A, double* Aint, unsigned long NBD, unsigned long NID);

   // METHOD=mm:3e68cbc7-5516-4454-ba92-054edd943941
   void MatrixMatrixSchurComplement(SymmetricMatrix& A, double* B, unsigned long NBD, unsigned long NID);

   // METHOD=mm:fcab3091-4b4e-49b7-aca8-02365f6ae746
   void MatrixVectorSchurComplement(SymmetricMatrix& A, double* B, unsigned long NBD, unsigned long NID);

   // METHOD=mm:bf30b211-3df6-11e0-891a-d4865a8b0383
   void VectorSchurComplement(double* A, double* B, unsigned long NBD, unsigned long NID);

   // METHOD=mm:9cd73cdd-2224-476b-9832-9bfe3af9e958
   void VectorSchurComplement(double* A, double* Aint, double* B, unsigned long NBD, unsigned long NID);

   // METHOD=mm:3887f33d-ca1e-47f7-8fee-33b76369afa4
   void RecoverInternalSolSchur(double* A, double* B, double* Usol, unsigned long NBD, unsigned long NID, double c = 1);

   // METHOD=mm:bde69226-cf9f-4f0e-8018-bf556deea2fe
   void RecoverInternalSolSchur(double* A, double* Aint, double* B, double* Usol, unsigned long NBD, unsigned long NID, double c = 1);

   // METHOD=mm:27d4429f-48bd-4844-aaed-5139d56cd011
   void LinearSystemSolutionContact(SymmetricMatrix& A, double* F, unsigned long LSNum);

   // METHOD=mm:22bbf330-8baf-42a3-8555-80edd424863e
   void LinearSystemSolutionContact(SymmetricMatrix& A, double* F, unsigned long LSNum, Vector& U);

   // METHOD=mm:c4032b09-bf50-45e5-b10e-2a5090078b84
   void SetLoadIncidSchur(unsigned long& NLoadedElem);


   // attr=mm:cdc99e7b-3df3-11e0-891a-d4865a8b0383
   TheoreticalSolution_S TheoSolu;

   // attr=mm:d3f9b745-0216-11e0-8f07-c00dff923bc5
   SolverParam_S SolvParam;

   // attr=mm:e70c8399-0216-11e0-8f07-c00dff923bc5
   Model GlobalModel;

   // attr=mm:b7f6499b-8365-4e33-8445-5d094c1b449c
   PartitionData PartitionInformations;

   // attr=mm:66207b0d-69a6-11e0-94a8-e19117dd202e
   OneIndexTable<double> ElemMassMatrices;

   // attr=mm:86e8c0ed-69a6-11e0-94a8-e19117dd202e
   OneIndexTable<double> ElemStiffMatrices;

   // attr=mm:61286a20-7dd5-4d7f-a75a-476714e0f47f
   OneIndexTable<double> ElemEffMatrices;

   // attr=mm:8fa04add-e777-4066-b96f-f65ba0cac569
   OneIndexTable<double> ElemInternalMatrices;

   // attr=mm:4e95bf01-69a7-11e0-94a8-e19117dd202e
   OneIndexTable<double> ElemLoadVectors;

   // attr=mm:c1eaeb35-69a6-11e0-94a8-e19117dd202e
   Vector PatchMeas;

   // attr=mm:f2f8c721-41f1-4477-81f5-08807dfa13ba
   Vector AuxCoords;

   // attr=mm:5a595d9a-4595-4d93-a57e-5420cb90e248
   SymmetricMatrix LocalMassMatrix;

   // attr=mm:3b0cd21a-fec2-4142-8b4d-6397f3593271
   SymmetricMatrix LocalStiffMatrix;

   // attr=mm:6636e96d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> U;
   // attr=mm:337658d5-2b9f-4582-bb7a-ce71732c70d5
   Vector Displacement;

   // attr=mm:7c44ab7d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> Upt;

   // attr=mm:7cdd3f7d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> U2pt;

   // attr=mm:aad3cc4d-6792-11e0-94a8-e19117dd202e
   TwoIndexTable<double> StrainSol;

   // attr=mm:1f7d4cdd-69a3-11e0-94a8-e19117dd202e
   TwoIndexTable<double> StressSol;

   // attr=mm:553a69c2-cb95-40b3-aa48-ee65533e9be2
   SymmetricMatrix ContactStiffMatrix;

   // attr=mm:f90bdba0-f25a-4fd4-b785-0417e3d2cc82
   SymmetricMatrix Kc;

   // attr=mm:3459caaf-69a3-11e0-94a8-e19117dd202e
   TwoIndexTable<double> ExtraSol;

   // attr=mm:89fa251d-3aa5-11e0-87c1-fb85dddc7efb
   Vector incU;
   // attr=mm:f584b794-148b-42b4-9b83-0017902e4ed6
   Vector Fint;

   // attr=mm:aac29087-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> F;

   // attr=mm:ad1c997a-d8a2-48b0-b4cc-ba3ee6c8a45c
   Vector FeAux;

   // attr=mm:ae55ff6d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> Fgr;

   // attr=mm:b44bf87d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> Press;

   // attr=mm:c62d971d-3aa5-11e0-87c1-fb85dddc7efb
   Vector Load;

   // attr=mm:d5e6da02-a1a7-4e98-b452-03b58d3abedf
   Vector ContactResidual;

   // attr=mm:1c826e4f-f24f-49e8-ab7d-a195e6699c9c
   Vector ContStressIP;

   // attr=mm:4a7e2973-5360-4359-a473-da898131ca63
   Vector ContStressCP;

   // attr=mm:04314480-0a4b-4f73-96d1-3bf82eca9b85
   Vector ContGap;

   // attr=mm:bbf3bbbf-5ff8-40d2-887e-2cc5f9cb5a78
   Vector ContCondition;

   // attr=mm:36b5b7cd-4818-4478-ba54-c42efd769414
   Vector StressPosition;

   // attr=mm:1ae28008-b97c-4df8-b1d9-0a94e5385a10
   BuiltInArray<long> LoadIncidSchur;

   // attr=mm:ca89923b-0051-4e6a-a9d9-fb0ecc0a3f92
   Vector ElemLoad;

   // attr=mm:24825d62-66be-4296-92bf-e1053e4a0760
   Vector InternalForce;

   // attr=mm:91714cbb-4685-4dc2-a7f1-fbaa7cd061d9
   Vector Damage;

   // attr=mm:c9c1074d-3aa5-11e0-87c1-fb85dddc7efb
   OneIndexTable<double> React;

   // attr=mm:a9df311d-3aa8-11e0-87c1-fb85dddc7efb
   char* PreFileName; //= NULL;

   // attr=mm:32b47cd1-3abe-11e0-87c1-fb85dddc7efb
   Matrix ErrorNorm;

   // attr=mm:492c4a49-3df6-11e0-891a-d4865a8b0383
   long SchurFlag; //= 0;

   // attr=mm:276ab52b-66fd-4672-b48f-17d3d4bdd38d
   SymmetricMatrix Schur_AA;

   // attr=mm:66e032a1-7c76-4c4f-b878-244d2d9a06dd
   SymmetricMatrix Schur_BB;

   // attr=mm:21621b55-a3b1-4163-bc43-cec9bbd99a49
   Matrix Schur_AB;

   // attr=mm:c79f2d2d-aa23-4e43-88ea-f77ec4973a32
   Matrix Schur_AB_Aux;

   // attr=mm:0acbb9ea-e67d-422c-a64e-8b76d60ebb38
   Matrix Schur_Mat_Aux;

   // attr=mm:c652c50e-ee8e-4bfb-903d-0998a8488bd4
   Vector Schur_Vec_Aux;

   // attr=mm:878f4370-2143-4b42-9383-610e2b6cc8bc
   Vector Schur_Vec_Aux1;

   // attr=mm:95927ba7-a722-4999-8e48-4d4acd58d563
   Vector Schur_fi;

   // attr=mm:1ac5744e-b07d-4c34-905b-54882eb1f30f
   Vector Schur_fb;

   // attr=mm:012a4ea4-cdea-4dc6-9716-d6ff32ab757e
   Vector LoadAux;

   // attr=mm:4e70700e-ffb1-4461-91ae-e26b98c2ad7b
   Vector InternalMassME;


 private:

//#UBLK-BEG-CLASSMEMB mm:8ed17010-020e-11e0-9503-bdd33c3b382f
//#UBLK-END-CLASSMEMB
};

//#UBLK-BEG-GLOBALH
//*** Removed methods: 14 Mar 2012 17:03
//---
//uuid=mm:75c8b72c-020e-11e0-9503-bdd33c3b382f
//   void Run();

//*** Removed methods: 03 Dec 2012 16:18
//---
//uuid=mm:3688029e-1526-4810-b7a4-fddae58b1cf1
//   void ComputeErrorNorm1DMatrices();

//*** Removed methods: 13 May 2014 16:59
//---
//uuid=mm:f5dfe496-c2b9-4339-a667-efca5c4a49f0
//   virtual void Run() = 0;

//*** Removed methods: 13 Mar 2015 16:16
//---
//uuid=mm:d8baf3f9-800f-4039-8a39-8d263cf50629
//   ContactSolver_E GetContactIntegrationMethod();
//---
//uuid=mm:65a885f1-ee67-4534-b593-fa1b89915198
//   ContactModel_E GetContactFrictionCondition();
//---
//uuid=mm:bd26a52f-8120-4c8d-8d9b-713e214bfa5c
//   unsigned long GetFrictionCoef();

//#UBLK-END-GLOBALH


#endif // _SOLVER_H_