NonlinearSolverTest Class Reference

#include <nonlin.h>

Inheritance diagram for NonlinearSolverTest:

Public Types
enum	Type { NONE, NORM, MINMAX, LASTNONLINEARSOLVERTEST }

Public Member Functions
virtual	~NonlinearSolverTest (void)
virtual doublereal	MakeTest (Solver *pS, const integer &Size, const VectorHandler &Vec, bool bResidual=false, doublereal dScaleAlgEqu=1., doublereal *pTestDiff=0)
virtual void	TestOne (doublereal &dRes, const VectorHandler &Vec, const integer &iIndex, doublereal dCoef) const =0
virtual void	TestMerge (doublereal &dResCurr, const doublereal &dResNew) const =0
virtual doublereal	TestPost (const doublereal &dRes) const
virtual const doublereal &	dScaleCoef (const integer &iIndex) const

Detailed Description

Definition at line 74 of file nonlin.h.

Member Enumeration Documentation

enum NonlinearSolverTest::Type

Enumerator
NONE
NORM
MINMAX
LASTNONLINEARSOLVERTEST

Definition at line 76 of file nonlin.h.

                  {
                NONE,
                NORM,
                MINMAX,

                LASTNONLINEARSOLVERTEST
        };

Constructor & Destructor Documentation

NonlinearSolverTest::~NonlinearSolverTest ( void  )

virtual

Definition at line 57 of file nonlin.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

const doublereal & NonlinearSolverTest::dScaleCoef ( const integer &  iIndex ) const

virtual

Reimplemented in NonlinearSolverTestRange, NonlinearSolverTestScaleMinMax, NonlinearSolverTestScaleNorm, and NonlinearSolverTestScale.

Definition at line 167 of file nonlin.cc.

References dOne.

Referenced by NonlinearSolverTestRange::dScaleCoef(), ImplicitStepIntegrator::TestScale(), and InverseDynamicsStepSolver::TestScale().

{
        return ::dOne;
}

doublereal NonlinearSolverTest::MakeTest ( Solver *  pS, const integer &  Size, const VectorHandler &  Vec, bool  bResidual = false, doublereal  dScaleAlgEqu = 1., doublereal *  pTestDiff = 0  )

virtual

Definition at line 63 of file nonlin.cc.

References ASSERT, SchurSolutionManager::ComplExchIntRes(), SchurSolutionManager::ComplExchIntSol(), DEBUGCOUTFNAME, DofOrder::DIFFERENTIAL, SchurDataManager::GetDofsList(), SchurDataManager::HowManyDofs(), VectorHandler::iGetSize(), SchurDataManager::LOCAL, order, Solver::pGetDataManager(), Solver::pGetSolutionManager(), SchurSolutionManager::StartExchIntRes(), SchurSolutionManager::StartExchIntSol(), TestOne(), and TestPost().

Referenced by NonlinearSolver::MakeResTest(), and NonlinearSolver::MakeSolTest().

{
        DEBUGCOUTFNAME("NonlinearSolverTestNorm::MakeTest");

        doublereal dTest = 0.;

        if (pTestDiff) {
                *pTestDiff = 0.;
        }

#ifdef USE_SCHUR
        /* Only residual test is parallelized; the master node
         * always knows the entire solution */
        ASSERT(pS != NULL);
        SchurSolutionManager *pSSM = dynamic_cast<SchurSolutionManager *>(pS->pGetSolutionManager());
        if (pSSM) {
                SchurDataManager *pSDM = 
                        dynamic_cast<SchurDataManager *>(pS->pGetDataManager());
                ASSERT(pSDM);

                integer iNumLocDofs = pSDM->HowManyDofs(SchurDataManager::LOCAL);
                integer *pLocDofs = pSDM->GetDofsList(SchurDataManager::LOCAL);

#if 0
                silent_cout("NonlinearSolverTest::MakeTest("
                                << MBDynComm.Get_rank() << ") "
                                "iNumLocDofs=" << iNumLocDofs << std::endl);
#endif
        
                if (bResidual) {
                        /*
                         * Chiama la routine di comunicazione per la trasmissione 
                         * del residuo delle interfacce
                         */
                        pSSM->StartExchIntRes();

                        /* calcola il test per i dofs locali */
                        for (int iCnt = 0; iCnt < iNumLocDofs; iCnt++) {
                                TestOne(dTest, Vec, pLocDofs[iCnt]);
                        }

                        /* collect contributions from other nodes,
                         * plus that of the interface; merge them according
                         * to the NonlinearSolverTest type */
                        pSSM->ComplExchIntRes(dTest, this);

                } else {
                        /*
                         * Chiama la routine di comunicazione per la trasmissione 
                         * del residuo delle interfacce
                         */
                        pSSM->StartExchIntSol();

                        /* calcola il test per i dofs locali */
                        for (int iCnt = 0; iCnt < iNumLocDofs; iCnt++) {
                                TestOne(dTest, Vec, pLocDofs[iCnt]);
                        }

                        /* collect contributions from other nodes,
                         * plus that of the interface; merge them according
                         * to the NonlinearSolverTest type */
                        pSSM->ComplExchIntSol(dTest, this);
                }

        } else
#endif // USE_SCHUR
        {
                ASSERT(Vec.iGetSize() == Size);
                const DataManager* const pDM = pS->pGetDataManager();

                for (int iCntp1 = 1; iCntp1 <= Size; iCntp1++) {
                        const DofOrder::Order order = pDM->GetEqType(iCntp1);
                        const doublereal dCoef = order == DofOrder::DIFFERENTIAL ? 1. : dScaleAlgEqu;

                        TestOne(dTest, Vec, iCntp1, dCoef);

                        if (pTestDiff && order == DofOrder::DIFFERENTIAL) {
                                TestOne(*pTestDiff, Vec, iCntp1, dCoef);
                        }
                }
        }

        if (pTestDiff) {
                *pTestDiff = TestPost(*pTestDiff);
        }

        return TestPost(dTest);
}

Here is the call graph for this function:

virtual void NonlinearSolverTest::TestMerge ( doublereal &  dResCurr, const doublereal &  dResNew  ) const

pure virtual

Implemented in NonlinearSolverTestRange, NonlinearSolverTestScaleMinMax, NonlinearSolverTestScaleNorm, NonlinearSolverTestMinMax, NonlinearSolverTestNorm, and NonlinearSolverTestNone.

Referenced by NonlinearSolverTestRange::TestMerge().

virtual void NonlinearSolverTest::TestOne ( doublereal &  dRes, const VectorHandler &  Vec, const integer &  iIndex, doublereal  dCoef  ) const

pure virtual

Implemented in NonlinearSolverTestRange, NonlinearSolverTestScaleMinMax, NonlinearSolverTestScaleNorm, NonlinearSolverTestMinMax, NonlinearSolverTestNorm, and NonlinearSolverTestNone.

Referenced by MakeTest(), and NonlinearSolverTestRange::TestOne().

doublereal NonlinearSolverTest::TestPost ( const doublereal &  dRes ) const

virtual

Reimplemented in NonlinearSolverTestRange, and NonlinearSolverTestNorm.

Definition at line 155 of file nonlin.cc.

References MBDYN_EXCEPT_ARGS.

Referenced by MakeTest(), and NonlinearSolverTestRange::TestPost().

{
        if (!std::isfinite(dRes)) {
                throw NonlinearSolver::ErrSimulationDiverged(MBDYN_EXCEPT_ARGS);
        }

        return dRes;
}

---

The documentation for this class was generated from the following files:

• nonlin.h
• nonlin.cc