Go to the documentation of this file.
45 const scalar coeffVal,
50 source_(mSize, sourceVal)
113 <<
"Different size matrices"
117 if (source_.size() != m.source_.size())
120 <<
"Different size source vectors"
124 scalarSquareMatrix::operator=(m);
142 m1.source_ + m2.source_
150 const simpleMatrix<Type>& m1,
151 const simpleMatrix<Type>& m2
154 return simpleMatrix<Type>
158 m1.source_ - m2.source_
167 const simpleMatrix<Type>& m
170 return simpleMatrix<Type>(
s*m.matrix_,
s*m.source_);
180 const simpleMatrix<Type>& m
183 os << static_cast<const scalarSquareMatrix&>(m) <<
nl << m.source_;
void operator=(const simpleMatrix< Type > &)
void LUsolve(scalarSquareMatrix &matrix, List< Type > &source)
A simple square matrix solver with scalar coefficients.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
Field< Type > LUsolve() const
Solve the matrix using LU decomposition with pivoting.
Generic templated field type.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
SolverPerformance< Type > solve(faMatrix< Type > &, Istream &)
Solve returning the solution statistics given convergence tolerance.
simpleMatrix(const label)
Construct given size.
SquareMatrix< scalar > scalarSquareMatrix
errorManip< error > abort(error &err)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Field< Type > solve() const
Solve the matrix using Gaussian elimination with pivoting.