42 assert(this->mpLinearSystem->rGetLhsMatrix() != NULL);
43 assert(this->mpLinearSystem->rGetRhsVector() != NULL);
50 mpMonodomainAssembler->SetMatrixToAssemble(this->mpLinearSystem->rGetLhsMatrix());
51 mpMonodomainAssembler->AssembleMatrix();
57 this->mpLinearSystem->FinaliseLhsMatrix();
76 index!= dist_vec_matrix_based.
End();
79 double V = distributed_current_solution[index];
87 dist_vec_matrix_based.
Restore();
92 MatMult(mMassMatrix, mVecForConstructingRhs, this->mpLinearSystem->rGetRhsVector());
101 mpNeumannSurfaceTermsAssembler->SetVectorToAssemble(this->mpLinearSystem->rGetRhsVector(),
false);
102 mpNeumannSurfaceTermsAssembler->AssembleVector();
105 this->mpLinearSystem->FinaliseRhsVector();
128 if (this->mpLinearSystem != NULL)
150 this->mpLinearSystem->SetMatrixIsSymmetric(
true);
155 Vec& r_template = this->mpLinearSystem->rGetRhsVector();
156 VecDuplicate(r_template, &mVecForConstructingRhs);
159 VecGetOwnershipRange(r_template, &ownership_range_lo, &ownership_range_hi);
160 PetscInt local_size = ownership_range_hi - ownership_range_lo;
162 this->mpMesh->CalculateMaximumNodeConnectivityPerProcess(),
163 local_size, local_size);
void SetupLinearSystem(Vec currentSolution, bool computeMatrix)
NaturalNeumannSurfaceTermAssembler< ELEMENT_DIM, SPACE_DIM, 1 > * mpNeumannSurfaceTermsAssembler
OperatorSplittingMonodomainSolver(AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > *pMesh, MonodomainTissue< ELEMENT_DIM, SPACE_DIM > *pTissue, BoundaryConditionsContainer< ELEMENT_DIM, SPACE_DIM, 1 > *pBoundaryConditions)
void InitialiseForSolve(Vec initialSolution)
void FollowingSolveLinearSystem(Vec currentSolution)
void PrepareForSetupLinearSystem(Vec currentSolution)
MonodomainAssembler< ELEMENT_DIM, SPACE_DIM > * mpMonodomainAssembler
MonodomainTissue< ELEMENT_DIM, SPACE_DIM > * mpMonodomainTissue