include/solvercontrol/GlobalSolver.h

#ifndef _GLOBALSOLVER_H_
#define _GLOBALSOLVER_H_
// ---------------------------------------------
// company   : 
// user      : 
// date      : 22 Mar 2015 17:25
// file      : GlobalSolver.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 "ds/SymmetricSparse.h"
#include "solvercontrol/Solver.h"


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

 public:

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

   // METHOD=mm:38d5c64f-69aa-11e0-94a8-e19117dd202e
   GlobalSolver(GlobalSolver& Instance);

   // METHOD=mm:6884bbea-69aa-11e0-94a8-e19117dd202e
   GlobalSolver& operator=(GlobalSolver& Instance);

   // METHOD=mm:17183bb1-0238-11e0-8f07-c00dff923bc5
   ~GlobalSolver();

   // METHOD=mm:75c8b863-020e-11e0-9503-bdd33c3b382f
   void Run();


 protected:

   // METHOD=mm:30b399be-022c-11e0-8f07-c00dff923bc5
   void NewtonRaphsonSolver();

   // METHOD=mm:3063b667-5591-47ad-a270-145584e5273c
   void NewtonRaphsonContactSolver();

   // METHOD=mm:593b83e3-022c-11e0-8f07-c00dff923bc5
   void DynamicLinearSolver();

   // METHOD=mm:a7e930ca-022c-11e0-8f07-c00dff923bc5
   void ProjectionSolver();

   // METHOD=mm:d830b816-022c-11e0-8f07-c00dff923bc5
   void StaticLinearSolver();

   // METHOD=mm:c8fd38c5-022b-11e0-8f07-c00dff923bc5
   void CentralDifferenceGlobalMethod();

   // METHOD=mm:1083c600-022c-11e0-8f07-c00dff923bc5
   void NewmarkGlobalMethod();

   // METHOD=mm:1af288c8-022d-11e0-8f07-c00dff923bc5
   void SetIterativeSparseTopology(SymmetricSparse& SPAIt, unsigned long NoBCFlag = 0);

   // METHOD=mm:d4e90856-6310-42f3-8782-3567539a519f
   void SetIterativeSparseTopology(SymmetricSparse& SPAItM, SymmetricSparse& SPAItK);

   // METHOD=mm:75c8b8fa-020e-11e0-9503-bdd33c3b382f
   void AssemblingGlbLoadVector(LoadSet& LoadSets, unsigned long LoadSetNumber, double  t = 0);

   // METHOD=mm:47f37d95-022e-11e0-8f07-c00dff923bc5
   void AssemblingGlbNonLinearPressureVector(unsigned long LoadSets, unsigned long LoadSetNumber, unsigned long Delta);

   // METHOD=mm:69ed0ebb-022e-11e0-8f07-c00dff923bc5
   void AssemblingGlbPressureVector(unsigned long LoadSets, unsigned long LoadSetNumber);

   // METHOD=mm:e8b2e6c6-aae7-4b6e-8831-56184ee1546d
   void AssemblingGlbContactData(LoadSet& LoadSets, unsigned long LoadSetNumber);

   // METHOD=mm:75c8b8f9-020e-11e0-9503-bdd33c3b382f
   void AssemblingGlobalStiffnessMatrix();

   // METHOD=mm:81c48e6f-022e-11e0-8f07-c00dff923bc5
   void AssemblingGlobalNonLinearStiffnessMatrix();

   // METHOD=mm:a22569fe-9a4c-4ee5-819d-d71ef2999d34
   void AssemblingGlobalContactStiffnessMatrix(LoadSet& LoadSets, unsigned long LoadSetNumber);

   // METHOD=mm:9f91d281-022e-11e0-8f07-c00dff923bc5
   void AssemblingGlobalStiffnessMatrixLoadVector();

   // METHOD=mm:75c8b8fb-020e-11e0-9503-bdd33c3b382f
   void AssemblingGlobalMassMatrix();

   // METHOD=mm:a082543b-3946-49c0-a301-c2905eada212
   void AssemblingGlobalEffectiveMatrix();

   // METHOD=mm:d36cfd79-022e-11e0-8f07-c00dff923bc5
   void ApplyNHDirichletBCLinear(double DeltaInc, unsigned short Iter = 0);

   // METHOD=mm:ea13631f-022e-11e0-8f07-c00dff923bc5
   void ApplyNHDirichletBCNonLinear();

   // METHOD=mm:75c8b868-020e-11e0-9503-bdd33c3b382f
   void LinearSystemSolution(SymmetricSparse& A, double* F, unsigned long LSNum);

   // METHOD=mm:da7fca40-98f5-4dbb-bf4d-46c195ab4537
   void LinearSystemSolution(SymmetricSparse& A, double* b, unsigned long LSNum, double* x);

   // METHOD=mm:d6152512-1336-4485-8eb7-263af49627d6
   void LinearSystemSolution(SymmetricSparse& A, double* b, unsigned long LSNum, Vector& x, unsigned long& NIT, double& NormResidual);


   // attr=mm:75c8b7c8-020e-11e0-9503-bdd33c3b382f
   SymmetricSparse StiffMatrix;

   // attr=mm:75c8b7c9-020e-11e0-9503-bdd33c3b382f
   SymmetricSparse MassMatrix;

   // attr=mm:322db6d5-8f73-4370-b0eb-0d3b26041043
   SymmetricSparse EffMatrix;


 private:

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

//#UBLK-BEG-GLOBALH
//*** Removed methods: 10 Nov 2014 17:07
//---
//uuid=mm:5b236e7c-1222-4c90-bd5c-be184b6ec897
//   void SetIterativeSparseTopology(SymmetricSparse& SPAIt, unsigned long Size);

//*** Removed methods: 13 Mar 2015 16:16
//---
//uuid=mm:956d1e74-a1d2-4151-8550-4678b60465f0
//   void NewtonRaphsonContactSolver();
//---
//uuid=mm:caa4a027-a301-4f32-8afe-5eb7b67943ea
//   void AssemblingGlbContactData(LoadSet& LoadSets, unsigned long LoadSetNumber);
//---
//uuid=mm:8da57291-031a-49d5-886c-19cf67689ae1
//   void AssemblingGlobalContactStiffnessMatrix(LoadSet& LoadSets, unsigned long LoadSetNumber);

//#UBLK-END-GLOBALH


#endif // _GLOBALSOLVER_H_