diff --git a/AMDiS/Reinit/src/BoundaryElementDist.h b/AMDiS/Reinit/src/BoundaryElementDist.h index 5088edfd9d3fb15e9bda7ce6a805683f6c9a48f7..f639392dd8bb9f67568e46ed5128321864f8f551 100644 --- a/AMDiS/Reinit/src/BoundaryElementDist.h +++ b/AMDiS/Reinit/src/BoundaryElementDist.h @@ -9,21 +9,17 @@ using namespace AMDiS; class BoundaryElementDist { - public: - MEMORY_MANAGED(BoundaryElementDist); +public: - /** - * Constructor. - */ BoundaryElementDist(ElementLevelSet *elLS_, int dim_) : dim(dim_), - elLS(elLS_) + elLS(elLS_) { FUNCNAME("BoundaryElementDist::BoundaryElementDist()"); TEST_EXIT(dim == 2 || dim == 3) ("function only works for dimension 2 !\n"); - }; + } /** * Virtual destructor. diff --git a/AMDiS/Reinit/src/BoundaryElementEdgeDist.h b/AMDiS/Reinit/src/BoundaryElementEdgeDist.h index 579fa634a6c5a6bd4650812c53635401db443720..3f6ccb2d0d0a4052350db4c21f3becfde22249a6 100644 --- a/AMDiS/Reinit/src/BoundaryElementEdgeDist.h +++ b/AMDiS/Reinit/src/BoundaryElementEdgeDist.h @@ -12,16 +12,12 @@ using namespace AMDiS; class BoundaryElementEdgeDist : public BoundaryElementDist { - public: - MEMORY_MANAGED(BoundaryElementEdgeDist); +public: - /** - * Constructor. - */ BoundaryElementEdgeDist(ElementLevelSet *elLS_, int dim_) : BoundaryElementDist(elLS_, dim_) - {}; + {} /** * Calculates distance from the interface for all vertices of a boundary diff --git a/AMDiS/Reinit/src/BoundaryElementLevelSetDist.h b/AMDiS/Reinit/src/BoundaryElementLevelSetDist.h index f639e207dd34f4bfea8a0015e31e377f4f4692f1..f5f3f976990c1fec31ed190a47ee9cc7ae0547d1 100644 --- a/AMDiS/Reinit/src/BoundaryElementLevelSetDist.h +++ b/AMDiS/Reinit/src/BoundaryElementLevelSetDist.h @@ -3,24 +3,18 @@ #include "ElInfo.h" #include "FixVec.h" - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" using namespace AMDiS; class BoundaryElementLevelSetDist : public BoundaryElementDist { - public: - MEMORY_MANAGED(BoundaryElementLevelSetDist); +public: - /** - * Constructor. - */ BoundaryElementLevelSetDist(ElementLevelSet *elLS_, int dim_) : BoundaryElementDist(elLS_, dim_) - {}; + {} /** * Calculates distance from the interface for all vertices of a boundary diff --git a/AMDiS/Reinit/src/BoundaryElementNormalDist.h b/AMDiS/Reinit/src/BoundaryElementNormalDist.h index a74c2725ac265750c8113369e6a907d55134dd92..5d82fd22a88c35af317c33c23b0e600b451fc011 100644 --- a/AMDiS/Reinit/src/BoundaryElementNormalDist.h +++ b/AMDiS/Reinit/src/BoundaryElementNormalDist.h @@ -3,24 +3,18 @@ #include "ElInfo.h" #include "FixVec.h" - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" using namespace AMDiS; class BoundaryElementNormalDist : public BoundaryElementDist { - public: - MEMORY_MANAGED(BoundaryElementNormalDist); +public: - /** - * Constructor. - */ BoundaryElementNormalDist(ElementLevelSet *elLS_, int dim_) : BoundaryElementDist(elLS_, dim_) - {}; + {} /** * Calculates distance from the interface for all vertices of a boundary diff --git a/AMDiS/Reinit/src/BoundaryElementTopDist.h b/AMDiS/Reinit/src/BoundaryElementTopDist.h index 0adfb23b35f6fb523536260b9f6441285aeab5b7..1724798f83c6bbd1ea8ae077c469f1429a72e6c5 100644 --- a/AMDiS/Reinit/src/BoundaryElementTopDist.h +++ b/AMDiS/Reinit/src/BoundaryElementTopDist.h @@ -3,9 +3,7 @@ #include "ElInfo.h" #include "FixVec.h" - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" #include "VelocityExt.h" @@ -13,24 +11,16 @@ using namespace AMDiS; class BoundaryElementTopDist : public BoundaryElementDist { - public: - MEMORY_MANAGED(BoundaryElementTopDist); - - /** - * Constructor. - */ +public: BoundaryElementTopDist(ElementLevelSet *elLS_, int dim_, VelocityExt *velExt_ = NULL) : BoundaryElementDist(elLS_, dim_), - velExt(velExt_) - {}; + velExt(velExt_) + {} - /** - * Destructor. - */ ~BoundaryElementTopDist() - {}; + {} /** * Calculates distance from the interface for all vertices of a boundary diff --git a/AMDiS/Reinit/src/ElementLevelSet.h b/AMDiS/Reinit/src/ElementLevelSet.h index 5bcf7964fe154057265f19e0750192d7ef4e04d9..a6292e7fd5d08ff0750c8739935f815f45f62259 100644 --- a/AMDiS/Reinit/src/ElementLevelSet.h +++ b/AMDiS/Reinit/src/ElementLevelSet.h @@ -1,17 +1,11 @@ #ifndef AMDIS_ELEMENTLEVELSET_H #define AMDIS_ELEMENTLEVELSET_H +#include "AMDiS_fwd.h" #include "ElementFunction.h" #include "FixVec.h" -// #include "MemoryManager.h" #include "Parameters.h" -namespace AMDiS { - class Element; - class ElInfo; - class Mesh; -} - using namespace AMDiS; using namespace std; @@ -58,26 +52,21 @@ using namespace std; class ElementLevelSet { public: - MEMORY_MANAGED(ElementLevelSet); - - /** - * Constructor. - */ ElementLevelSet(const char *name_, ElementFunction<double> *lSFct_, Mesh *mesh_) : name(name_), - elInfo(NULL), - lastEl(NULL), - level_set_domain(LEVEL_SET_UNDEFINED), - numIntersecPoints(0), - elStatus(LEVEL_SET_UNDEFINED), - numElVertexInterior(0), - numElVertexBoundary(0), - numElVertexExterior(0), - LS_VAL_TOL(1.e-8), - LS_VAL_MIN(1.e-8), - SP_BARY_TOL(1.e-7) + elInfo(NULL), + lastEl(NULL), + level_set_domain(LEVEL_SET_UNDEFINED), + numIntersecPoints(0), + elStatus(LEVEL_SET_UNDEFINED), + numElVertexInterior(0), + numElVertexBoundary(0), + numElVertexExterior(0), + LS_VAL_TOL(1.e-8), + LS_VAL_MIN(1.e-8), + SP_BARY_TOL(1.e-7) { FUNCNAME("ElementLevelSet::ElementLevelSet()"); @@ -108,11 +97,8 @@ class ElementLevelSet TEST_EXIT(LS_VAL_MIN > 0)("illegal LS_VAL_MIN\n"); TEST_EXIT(SP_BARY_TOL > 0)("illegal SP_BARY_TOL\n"); } - }; + } - /** - * Destructor. - */ ~ElementLevelSet() { if (elVertexStatusVec) @@ -121,7 +107,7 @@ class ElementLevelSet delete [] elVertexLevelSetVec; if (elIntersecPoints) DELETE elIntersecPoints; - }; + } /** * Calculates LevelSet-status of element and its intersection points diff --git a/AMDiS/Reinit/src/ElementUpdate.h b/AMDiS/Reinit/src/ElementUpdate.h index c20fb738046e3f42b1ba50a4dcfb5c98dd7b5373..a62839d0ea75ac066e92ae766461e212332c0460 100644 --- a/AMDiS/Reinit/src/ElementUpdate.h +++ b/AMDiS/Reinit/src/ElementUpdate.h @@ -2,28 +2,18 @@ #define ELEMENTUPDATE_H #include "FixVec.h" -// #include "MemoryManager.h" - #include "VelocityExt.h" using namespace AMDiS; class ElementUpdate { - public: - MEMORY_MANAGED(ElementUpdate); - - /** - * Constructor. - */ +public: ElementUpdate(VelocityExt *velExt_) : velExt(velExt_) - {}; + {} - /** - * Virtual destructor. - */ - virtual ~ElementUpdate() {}; + virtual ~ElementUpdate() {} /** * Pure virtual function. diff --git a/AMDiS/Reinit/src/ElementUpdate_2d.h b/AMDiS/Reinit/src/ElementUpdate_2d.h index a5fbbd964c94e0fd2287ed0831f59575c53b28cd..3902bb8add68332445218fc37516b28490d75366 100644 --- a/AMDiS/Reinit/src/ElementUpdate_2d.h +++ b/AMDiS/Reinit/src/ElementUpdate_2d.h @@ -2,8 +2,6 @@ #define ELEMENTUPDATE_2D_H #include "FixVec.h" -// #include "MemoryManager.h" - #include "ElementUpdate.h" #include "VelocityExt.h" @@ -11,15 +9,10 @@ using namespace AMDiS; class ElementUpdate_2d : public ElementUpdate { - public: - MEMORY_MANAGED(ElementUpdate_2d); - - /** - * Constructor. - */ +public: ElementUpdate_2d(VelocityExt *velExt_ = NULL) : ElementUpdate(velExt_) - {}; + {} /** * Realization of ElementUpdate::calcElementUpdate. diff --git a/AMDiS/Reinit/src/ElementUpdate_3d.h b/AMDiS/Reinit/src/ElementUpdate_3d.h index 1801dd731eb33a8d15ab77a7b64dc8dd3c009d26..c58f4dab20cba77de1101150b67bcb0cf176a222 100644 --- a/AMDiS/Reinit/src/ElementUpdate_3d.h +++ b/AMDiS/Reinit/src/ElementUpdate_3d.h @@ -2,8 +2,6 @@ #define ELEMENTUPDATE_3D_H #include "FixVec.h" -// #include "MemoryManager.h" - #include "ElementUpdate.h" #include "ElementUpdate_2d.h" #include "VelocityExt.h" @@ -12,25 +10,17 @@ using namespace AMDiS; class ElementUpdate_3d : public ElementUpdate { - public: - MEMORY_MANAGED(ElementUpdate_3d); - - /** - * Constructor. - */ +public: ElementUpdate_3d(VelocityExt *velExt_ = NULL) : ElementUpdate(velExt_) { - elUpdate2d = NEW ElementUpdate_2d(velExt_); - }; + elUpdate2d = new ElementUpdate_2d(velExt_); + } - /** - * Destructor. - */ ~ElementUpdate_3d() { - DELETE elUpdate2d; - }; + delete elUpdate2d; + } /** * Realization of ElementUpdate::calcElementUpdate. diff --git a/AMDiS/Reinit/src/HL_SignedDist.h b/AMDiS/Reinit/src/HL_SignedDist.h index 0853c6f41a9d92223f003e889ec71f81cec63b5c..ed360609f8f652f3c674ba067a781d15a1a8ad19 100644 --- a/AMDiS/Reinit/src/HL_SignedDist.h +++ b/AMDiS/Reinit/src/HL_SignedDist.h @@ -8,11 +8,8 @@ #include "FileWriter.h" #include "FixVec.h" #include "Flag.h" -// #include "MemoryManager.h" #include "Parameters.h" - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" #include "BoundaryElementLevelSetDist.h" #include "BoundaryElementTopDist.h" @@ -41,12 +38,7 @@ using namespace AMDiS; ////////////////////////////////////////////////////////////////////////////// class HL_SignedDist { - public: - MEMORY_MANAGED(HL_SignedDist); - - /** - * Constructor. - */ +public: HL_SignedDist(const char *name_, int dim_, bool doVelocityExt = false, @@ -77,11 +69,11 @@ class HL_SignedDist // ===== Create functionality for velocity extension. ===== if (doVelocityExt) { if (velExtType.isSet(VEL_EXT)) - velExt = NEW VelocityExt(dim); + velExt = new VelocityExt(dim); else - velExt = NEW VelocityExtFromVelocityField(dim); + velExt = new VelocityExtFromVelocityField(dim); } - }; + } /** * Virtual destructor. @@ -89,19 +81,19 @@ class HL_SignedDist virtual ~HL_SignedDist() { if (elUpdate) - DELETE elUpdate; + delete elUpdate; if (bndElDist) - DELETE bndElDist; + delete bndElDist; if (elLS) - DELETE elLS; + delete elLS; if (phi) - DELETE phi; + delete phi; if (bound_DOF) - DELETE bound_DOF; + delete bound_DOF; - DELETE velExt; + delete velExt; }; /** @@ -251,13 +243,13 @@ class HL_SignedDist */ void printLevelSetFct() { - FileWriter *fileWriter = NEW FileWriter( + FileWriter *fileWriter = new FileWriter( "SignedDist->level set fct output", feSpace->getMesh(), const_cast<DOFVector<double> *>(lS_DOF)); fileWriter->writeFiles(adaptInfo, false); - DELETE fileWriter; + delete fileWriter; }; /** @@ -265,13 +257,13 @@ class HL_SignedDist */ void printSignedDistFct() { - FileWriter *fileWriter = NEW FileWriter( + FileWriter *fileWriter = new FileWriter( "SignedDist->result output", feSpace->getMesh(), sD_DOF); fileWriter->writeFiles(adaptInfo, false); - DELETE fileWriter; + delete fileWriter; }; protected: @@ -306,13 +298,12 @@ class HL_SignedDist */ void printBoundInitFct() { - FileWriter *fileWriter = NEW FileWriter( - "SignedDist->boundary initialization output", - feSpace->getMesh(), - sD_DOF); + FileWriter *fileWriter = new FileWriter("SignedDist->boundary initialization output", + feSpace->getMesh(), + sD_DOF); fileWriter->writeFiles(adaptInfo, false); - DELETE fileWriter; + delete fileWriter; }; public: diff --git a/AMDiS/Reinit/src/HL_SignedDistBornemann.h b/AMDiS/Reinit/src/HL_SignedDistBornemann.h index 648e97a7abd2fb72d9648838a43d4dbcc9c1418e..743095bc4f16dbddf4143bae5a987a750701923c 100644 --- a/AMDiS/Reinit/src/HL_SignedDistBornemann.h +++ b/AMDiS/Reinit/src/HL_SignedDistBornemann.h @@ -1,19 +1,13 @@ #ifndef HL_SIGNEDDISTBORNEMANN #define HL_SIGNEDDISTBORNEMANN +#include <queue> +#include <time.h> #include "ElInfo.h" #include "FixVec.h" -// #include "MemoryManager.h" #include "Traverse.h" -#include <queue> -#include <time.h> - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" -/* #include "BoundaryElementTopDist.h" */ -/* #include "BoundaryElementEdgeDist.h" */ -/* #include "BoundaryElementNormalDist.h" */ #include "ElementUpdate.h" #include "ElementUpdate_2d.h" #include "ElementUpdate_3d.h" @@ -32,25 +26,20 @@ typedef struct class HL_SignedDistBornemann : public HL_SignedDist { - public: - MEMORY_MANAGED(HL_SignedDistBornemann); - - /** - * Constructor. - */ +public: HL_SignedDistBornemann(const char *name_,int dim_) : HL_SignedDist(name_, dim_), - smiAdapter(NULL) + smiAdapter(NULL) { FUNCNAME("HL_SignedDistBornemann::HL_SignedDistBornemann"); // ===== Read parameters from init file. ===== - GET_PARAMETER(0,name + "->tolerance", "%f", &tol); - GET_PARAMETER(0,name + "->count_how_often_saved_in_list", "%d", &count_in_list); - GET_PARAMETER(0,name + "->save_in_list->the ..th", "%d", &print_in_list); + GET_PARAMETER(0, name + "->tolerance", "%f", &tol); + GET_PARAMETER(0, name + "->count_how_often_saved_in_list", "%d", &count_in_list); + GET_PARAMETER(0, name + "->save_in_list->the ..th", "%d", &print_in_list); TEST_EXIT(tol > 0)("illegal tolerance !\n"); - }; + } protected: /** diff --git a/AMDiS/Reinit/src/HL_SignedDistLevels.h b/AMDiS/Reinit/src/HL_SignedDistLevels.h index aa02a31d43f6cd7ce44616d67631ff3e0023ff85..7c3e5fba3cd7ada81e62a7420c6c3ea1dc68a861 100644 --- a/AMDiS/Reinit/src/HL_SignedDistLevels.h +++ b/AMDiS/Reinit/src/HL_SignedDistLevels.h @@ -1,23 +1,19 @@ #ifndef HL_SIGNEDDISTLEVELS #define HL_SIGNEDDISTLEVELS -#include "ElInfo.h" -#include "FixVec.h" -// #include "MemoryManager.h" -#include "Traverse.h" #include <queue> #include <vector> #include <time.h> - +#include "ElInfo.h" +#include "FixVec.h" +#include "Traverse.h" #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" #include "ElementUpdate.h" #include "ElementUpdate_2d.h" #include "ElementUpdate_3d.h" #include "HL_SignedDist.h" #include "VelocityExt.h" - #include "SMIAdapter.h" #include "smi.h" @@ -31,31 +27,26 @@ typedef struct class HL_SignedDistLevels : public HL_SignedDist { - public: - MEMORY_MANAGED(HL_SignedDistLevels); - - /** - * Constructor. - */ +public: HL_SignedDistLevels(const char *name_, int dim_, bool doVelocityExt = false, Flag velExtType_ = VEL_EXT) : HL_SignedDist(name_, dim_, doVelocityExt, velExtType_), - smiAdapter(NULL) + smiAdapter(NULL) { - FUNCNAME("HL_SignedDistLevels::HL_SignedDistLevels"); + FUNCNAME("HL_SignedDistLevels::HL_SignedDistLevels()"); // ===== Read parameters from init file. ===== - GET_PARAMETER(0,name + "->tolerance", "%f", &tol); - GET_PARAMETER(0,name + "->count_how_often_saved_in_list", "%d", &count_in_list); - GET_PARAMETER(0,name + "->save_in_list->the ..th", "%d", &print_in_list); - GET_PARAMETER(0,name + "->save_in_list->after the ..th traversing of the list", "%d", &print_in_list_2); - GET_PARAMETER(0,name + "->print_level", "%d", &chosen_level); - GET_PARAMETER(0,name + "->print_level_yes_no", "%d", &print_level); + GET_PARAMETER(0, name + "->tolerance", "%f", &tol); + GET_PARAMETER(0, name + "->count_how_often_saved_in_list", "%d", &count_in_list); + GET_PARAMETER(0, name + "->save_in_list->the ..th", "%d", &print_in_list); + GET_PARAMETER(0, name + "->save_in_list->after the ..th traversing of the list", "%d", &print_in_list_2); + GET_PARAMETER(0, name + "->print_level", "%d", &chosen_level); + GET_PARAMETER(0, name + "->print_level_yes_no", "%d", &print_level); TEST_EXIT(tol > 0)("illegal tolerance !\n"); - }; + } protected: /** diff --git a/AMDiS/Reinit/src/HL_SignedDistTraverse.cc b/AMDiS/Reinit/src/HL_SignedDistTraverse.cc index 534db935c3f8baf78f48159308bdaf3e22c7b73f..ace9bb587c0810defbd8d3d775062ba17f6c1497 100644 --- a/AMDiS/Reinit/src/HL_SignedDistTraverse.cc +++ b/AMDiS/Reinit/src/HL_SignedDistTraverse.cc @@ -1,9 +1,7 @@ #include "HL_SignedDistTraverse.h" - #include "VelocityExtFromVelocityField.h" -void -HL_SignedDistTraverse::initializeBoundary() +void HL_SignedDistTraverse::initializeBoundary() { FUNCNAME("HL_SignedDistTraverse::initializeBoundary()"); @@ -16,7 +14,7 @@ HL_SignedDistTraverse::initializeBoundary() int elStatus; const int nBasFcts = feSpace->getBasisFcts()->getNumber(); - DegreeOfFreedom *locInd = GET_MEMORY(DegreeOfFreedom, nBasFcts); + DegreeOfFreedom *locInd = new DegreeOfFreedom[nBasFcts]; ElInfo *elInfo; if (velExt && velExtType.isSet(VEL_EXT_FROM_VEL_FIELD)) { @@ -26,8 +24,7 @@ HL_SignedDistTraverse::initializeBoundary() Mesh::FILL_BOUND | Mesh::FILL_COORDS | Mesh::FILL_GRD_LAMBDA); - } - else { + } else { elInfo = stack.traverseFirst(feSpace->getMesh(), -1, Mesh::CALL_LEAF_EL | @@ -94,9 +91,7 @@ HL_SignedDistTraverse::initializeBoundary() elInfo = stack.traverseNext(elInfo); } // end of: mesh traverse - FREE_MEMORY(locInd, DegreeOfFreedom, nBasFcts); - - return; + delete [] locInd; } void @@ -165,7 +160,7 @@ HL_SignedDistTraverse::HL_elementUpdate(ElInfo *elInfo) // ===== Get global indices of vertices of element. ===== const int nBasFcts = feSpace->getBasisFcts()->getNumber(); - DegreeOfFreedom *locInd = GET_MEMORY(DegreeOfFreedom, nBasFcts); + DegreeOfFreedom *locInd = new DegreeOfFreedom[nBasFcts]; feSpace->getBasisFcts()->getLocalIndices( const_cast<Element *>(elInfo->getElement()), @@ -203,13 +198,12 @@ HL_SignedDistTraverse::HL_elementUpdate(ElInfo *elInfo) } } - FREE_MEMORY(locInd, DegreeOfFreedom, nBasFcts); + delete [] locInd; } -double -HL_SignedDistTraverse::calcElementUpdate(ElInfo *elInfo, - int nXh, - const DegreeOfFreedom *locInd) +double HL_SignedDistTraverse::calcElementUpdate(ElInfo *elInfo, + int nXh, + const DegreeOfFreedom *locInd) { double update; FixVec<WorldVector<double> *, VERTEX> elVert(dim, NO_INIT); @@ -257,8 +251,7 @@ HL_SignedDistTraverse::calcElementUpdate(ElInfo *elInfo, return update; } -bool -HL_SignedDistTraverse::checkTol() +bool HL_SignedDistTraverse::checkTol() { DOFVector<double>::Iterator it_sD(sD_DOF, USED_DOFS); DOFVector<double>::Iterator it_sDOld(sDOld_DOF, USED_DOFS); diff --git a/AMDiS/Reinit/src/HL_SignedDistTraverse.h b/AMDiS/Reinit/src/HL_SignedDistTraverse.h index 2f41d015f7f4997319c52df16e044b4dda250bc4..bb0fe7c443d7ca0d830fa349901cf0d735a9aabc 100644 --- a/AMDiS/Reinit/src/HL_SignedDistTraverse.h +++ b/AMDiS/Reinit/src/HL_SignedDistTraverse.h @@ -3,11 +3,8 @@ #include "ElInfo.h" #include "FixVec.h" -//#include "MemoryManager.h" #include "Traverse.h" - #include "ElementLevelSet.h" - #include "BoundaryElementDist.h" #include "ElementUpdate.h" #include "ElementUpdate_2d.h" @@ -19,20 +16,15 @@ using namespace AMDiS; class HL_SignedDistTraverse : public HL_SignedDist { - public: - MEMORY_MANAGED(HL_SignedDistTraverse); - - /** - * Constructor. - */ +public: HL_SignedDistTraverse(const char *name_, int dim_, bool doVelocityExt = false, Flag velExtType_ = VEL_EXT) : HL_SignedDist(name_, dim_, doVelocityExt, velExtType_), - sDOld_DOF(NULL), - update_DOF(NULL), - tol_reached(false) + sDOld_DOF(NULL), + update_DOF(NULL), + tol_reached(false) { FUNCNAME("HL_SignedDistTraverse::HL_SignedDistTraverse"); @@ -49,11 +41,8 @@ class HL_SignedDistTraverse : public HL_SignedDist calcUpdate_Cntr = 0; setUpdate_Cntr = 0; // ---> end: for test purposes - }; + } - /** - * Destructor. - */ ~HL_SignedDistTraverse() { if (sDOld_DOF) @@ -62,7 +51,7 @@ class HL_SignedDistTraverse : public HL_SignedDist // ---> for test purposes: print result of update counting printUpdateCntr(); // ---> end: for test purposes - }; + } protected: /** diff --git a/AMDiS/Reinit/src/NormEps.h b/AMDiS/Reinit/src/NormEps.h index 1e83953b05a5173601857fad63d1c6de3eda7b61..27d6a3fc99ab0775e83ed51b55ac4568014ce962 100644 --- a/AMDiS/Reinit/src/NormEps.h +++ b/AMDiS/Reinit/src/NormEps.h @@ -3,43 +3,33 @@ #include "FixVec.h" #include "Global.h" -// #include "MemoryManager.h" #include "Parameters.h" using namespace AMDiS; class NormEps { - public: - MEMORY_MANAGED(NormEps); +public: - /** - * Calculates regularized norm of WorldVector x. - */ + /// Calculates regularized norm of WorldVector x. static double calcNormEps(const WorldVector<double> &x); - /** - * Calculates square of regularized norm of WorldVector x. - */ + /// Calculates square of regularized norm of WorldVector x. static double calcNormEps2(const WorldVector<double> &x); - /** - * Set regularization epsilon eps. - */ - static inline void setEps() { + /// Set regularization epsilon eps. + static inline void setEps() + { FUNCNAME("NormEps::setEps()"); eps = 0.0; GET_PARAMETER(0, "NormEps->epsilon for norm regularization", "%f", &eps); - TEST_EXIT(eps > 1.e-15) - ("illegal eps for norm regularization !\n"); - }; + TEST_EXIT(eps > 1.e-15)("illegal eps for norm regularization !\n"); + } - protected: - /** - * Epsilon for regularization. - */ +protected: + /// Epsilon for regularization. static double eps; }; diff --git a/AMDiS/Reinit/src/VelocityExt.h b/AMDiS/Reinit/src/VelocityExt.h index 0260fca24015de04f0a2ee857258232fa478b16e..953e89c620184354b77b75c9c1abd1915eed59c2 100644 --- a/AMDiS/Reinit/src/VelocityExt.h +++ b/AMDiS/Reinit/src/VelocityExt.h @@ -6,36 +6,25 @@ #include "Element.h" #include "FileWriter.h" #include "FixVec.h" -// #include "MemoryManager.h" using namespace AMDiS; class VelocityExt { - public: - MEMORY_MANAGED(VelocityExt); - - /** - * Constructor. - */ +public: VelocityExt(int dim_) : dim(dim_), - nVelDOFs(0), - lamVec(dim_,dim_+1,NO_INIT), - permutation(dim_,NO_INIT) + nVelDOFs(0), + lamVec(dim_, dim_ + 1, NO_INIT), + permutation(dim_, NO_INIT) { indexFace = -1; - }; + } - /** - * Destructor. - */ virtual ~VelocityExt() - {}; + {} - /** - * Print velDOF to file. - */ + /// Print velDOF to file. void printVelDOF(AdaptInfo *adaptInfo, int i = 0) { FUNCNAME("VelocityExt::printVelDOF()"); @@ -49,11 +38,9 @@ class VelocityExt fileWriter->writeFiles(adaptInfo, false); DELETE fileWriter; - }; + } - /** - * Print origVelDOF to file. - */ + /// Print origVelDOF to file. void printOrigVelDOF(AdaptInfo *adaptInfo, int i = 0) { FUNCNAME("VelocityExt::printOrigVelDOF()"); diff --git a/AMDiS/Reinit/src/VelocityExtFromVelocityField.h b/AMDiS/Reinit/src/VelocityExtFromVelocityField.h index d606a982b1468ac5897277140d95551be9f8203f..9d90b0b17d0db3a3fadfd2852a301d042628bf5a 100644 --- a/AMDiS/Reinit/src/VelocityExtFromVelocityField.h +++ b/AMDiS/Reinit/src/VelocityExtFromVelocityField.h @@ -5,107 +5,36 @@ #include "BasisFunction.h" #include "DOFVector.h" #include "ElInfo.h" -// #include "MemoryManager.h" - #include "NormEps.h" #include "VelocityExt.h" using namespace AMDiS; - -/* #ifndef NORMEPS_H */ -/* #define NORMEPS_H */ -/* ///////////////////////////////////////////////////////////////////////////// */ -/* // c l a s s N o r m E p s // */ -/* ///////////////////////////////////////////////////////////////////////////// */ -/* class NormEps */ -/* { */ -/* public: */ -/* MEMORY_MANAGED(NormEps); */ - -/* /\** */ -/* * Calculates regularized norm of WorldVector x. */ -/* *\/ */ -/* static double calcNormEps(const WorldVector<double> &x) */ -/* { */ -/* FUNCNAME("NormEps::calcNormEps()"); */ - -/* static double result; */ - -/* TEST_EXIT(eps > 1.e-15)("illegal eps for norm regularization !\n"); */ -/* result = x*x + eps*eps; */ - -/* return(std::sqrt(result)); */ -/* }; */ - -/* /\** */ -/* * Calculates square of regularized norm of WorldVector x. */ -/* *\/ */ -/* static double calcNormEps2(const WorldVector<double> &x) */ -/* { */ -/* FUNCNAME("NormEps::calcNormEps2()"); */ - -/* static double result; */ - -/* TEST_EXIT(eps > 1.e-15)("illegal eps for norm regularization !\n"); */ -/* result = x*x + eps*eps; */ - -/* return(result); */ -/* }; */ - -/* /\** */ -/* * Set regularization epsilon eps. */ -/* *\/ */ -/* static inline void setEps() { */ -/* FUNCNAME("NormEps::setEps()"); */ - -/* eps = 0.0; */ -/* GET_PARAMETER(0, "NormEps->epsilon for norm regularization", "%f", &eps); */ - -/* TEST_EXIT(eps > 1.e-15) */ -/* ("illegal eps for norm regularization !\n"); */ -/* }; */ - -/* protected: */ -/* /\** */ -/* * Epsilon for regularization. */ -/* *\/ */ -/* static double eps; */ -/* }; */ -/* #endif // NORMEPS_H */ - ///////////////////////////////////////////////////////////////////////////// // c l a s s V e l o c i t y E x t F r o m V e l o c i t y F i e l d // ///////////////////////////////////////////////////////////////////////////// class VelocityExtFromVelocityField : public VelocityExt { - public: - MEMORY_MANAGED(VelocityExtFromVelocityField); +public: - /** - * Constructor. - */ VelocityExtFromVelocityField(int dim_) : VelocityExt(dim_), - lSFct(NULL), - elNormalVel(dim_, NO_INIT), - basFcts(NULL) + lSFct(NULL), + elNormalVel(dim_, NO_INIT), + basFcts(NULL) { FUNCNAME("VelocityExtFromVelocityField::VelocityExtFromVelocityField()"); - - lSFctVal = GET_MEMORY(double, dim+1); - + + lSFctVal = new double[dim + 1]; + // ===== set epsilon for norm regularization ===== NormEps::setEps(); - }; + } - /** - * Destructor. - */ ~VelocityExtFromVelocityField() { - FREE_MEMORY(lSFctVal, double, dim+1); - }; + delete [] lSFctVal; + } /** * Set velocity field. @@ -141,9 +70,10 @@ class VelocityExtFromVelocityField : public VelocityExt /** * Sets elInfo. */ - inline void setElInfo(ElInfo *elInfo_) { + inline void setElInfo(ElInfo *elInfo_) + { elInfo = elInfo_; - }; + } protected: /**