ILU0Solver Class Reference

#include <ILU0Solver.h>

Inheritance diagram for ILU0Solver:

Collaboration diagram for ILU0Solver:

Public Member Functions
	ILU0Solver ()
virtual	~ILU0Solver ()
virtual MFRPtr	solve (LinearSystem &ls)
virtual void	cleanup ()
virtual void	smooth (LinearSystem &ls)
	DEFINE_TYPENAME ("ILU0Solver")
Public Member Functions inherited from LinearSolver
	LinearSolver ()

Private Member Functions
void	doSweeps (LinearSystem &ls, const int nSweeps)
	ILU0Solver (const ILU0Solver &)

Additional Inherited Members
Public Attributes inherited from LinearSolver
int	nMaxIterations
int	verbosity
double	relativeTolerance
double	absoluteTolerance

Detailed Description

Solve a linear system using Jacobi iterations

Definition at line 20 of file ILU0Solver.h.

Constructor & Destructor Documentation

ILU0Solver::ILU0Solver

(

)

Definition at line 15 of file ILU0Solver.cpp.

References logCtor.

{
  logCtor();
  
}

ILU0Solver::~ILU0Solver ( )

virtual

Definition at line 21 of file ILU0Solver.cpp.

References logDtor.

{
  logDtor();
}

ILU0Solver::ILU0Solver ( const ILU0Solver &  )

private

Member Function Documentation

void ILU0Solver::cleanup ( )

virtual

Implements LinearSolver.

Definition at line 42 of file ILU0Solver.cpp.

{
}

ILU0Solver::DEFINE_TYPENAME

(

"ILU0Solver"

)

void ILU0Solver::doSweeps ( LinearSystem &  ls, const int  nSweeps  )

private

Definition at line 27 of file ILU0Solver.cpp.

References LinearSystem::getB(), LinearSystem::getDelta(), LinearSystem::getMatrix(), LinearSystem::getResidual(), and MultiFieldMatrix::iluSolve().

Referenced by smooth(), and solve().

{
  const MultiFieldMatrix& m = ls.getMatrix();
  MultiField& delta = ls.getDelta();
  const MultiField& b = ls.getB();
  MultiField& r = ls.getResidual();

  for(int i=0; i<nSweeps; i++)
  {
      m.iluSolve(delta,b,r);
  }
}

void ILU0Solver::smooth ( LinearSystem &  ls )

virtual

Implements LinearSolver.

Definition at line 98 of file ILU0Solver.cpp.

References doSweeps().

{
  doSweeps(ls,1);
}

MFRPtr ILU0Solver::solve ( LinearSystem &  ls )

virtual

Implements LinearSolver.

Definition at line 47 of file ILU0Solver.cpp.

References LinearSolver::absoluteTolerance, MultiFieldMatrix::computeResidual(), doSweeps(), LinearSystem::getB(), LinearSystem::getDelta(), LinearSystem::getMatrix(), MultiField::getOneNorm(), LinearSystem::getResidual(), LinearSolver::nMaxIterations, LinearSolver::relativeTolerance, and LinearSolver::verbosity.

{
  ls.getMatrix().computeResidual(ls.getDelta(),
                                 ls.getB(),
                                 ls.getResidual());
  MFRPtr rNorm0(ls.getResidual().getOneNorm());

#ifdef FVM_PARALLEL
  if (verbosity >0 && MPI::COMM_WORLD.Get_rank() == 0 )
    cout << "0: " << *rNorm0 << "procID = " << MPI::COMM_WORLD.Get_rank() << endl;
#endif

#ifndef FVM_PARALLEL
   if ( verbosity > 0 )
      cout << "0: " << *rNorm0 << endl;
#endif

  if (*rNorm0 < absoluteTolerance )
    return rNorm0;

  for(int i=1; i<nMaxIterations; i++)
  {
      doSweeps(ls,1);

      ls.getMatrix().computeResidual(ls.getDelta(),
                                     ls.getB(),
                                     ls.getResidual());
      MFRPtr rNorm(ls.getResidual().getOneNorm());
      MFRPtr normRatio((*rNorm)/(*rNorm0));

#ifndef FVM_PARALLEL
      if (verbosity >0  )
        cout << i << ": " <<  *rNorm << endl;
#endif


#ifdef FVM_PARALLEL
     if (*rNorm < absoluteTolerance || *normRatio < relativeTolerance || i == nMaxIterations-1)
        if (verbosity >0 && MPI::COMM_WORLD.Get_rank() == 0  )
        cout <<i << ": " << "procID = " << MPI::COMM_WORLD.Get_rank() <<  *rNorm << endl;
#endif

      if (*rNorm < absoluteTolerance || *normRatio < relativeTolerance)
         break;
      
  }
  return rNorm0;
}

---

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

• ILU0Solver.h
• ILU0Solver.cpp