42 const scalar coeffVal,
47 source_(mSize, sourceVal)
105 <<
"Attempted assignment to self" 112 <<
"Different size matrices" 116 if (source_.size() != m.source_.size())
119 <<
"Different size source vectors" 140 + static_cast<const scalarSquareMatrix&>(m2),
141 m1.source_ + m2.source_
156 - static_cast<const scalarSquareMatrix&>(m2),
157 m1.source_ - m2.source_
182 os << static_cast<const scalarSquareMatrix&>(
m) <<
nl <<
m.source_;
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
void operator=(const simpleMatrix< Type > &)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
void operator=(const zero)
Assignment of all elements to zero.
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().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.
Pre-declare SubField and related Field type.
A simple square matrix solver with scalar coefficients.
errorManip< error > abort(error &err)
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.
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
label m() const
Return the number of rows.
simpleMatrix(const label)
Construct given size.
SquareMatrix< scalar > scalarSquareMatrix
void LUsolve(scalarSquareMatrix &matrix, List< Type > &source)
Solve the matrix using LU decomposition with pivoting.