Lines Matching refs:lsys
604 linearSystemPETSc<double> *lsys = new linearSystemPETSc<double>; in computeParametrization() local
613 lsys->setParameter("petsc_solver_options", options); in computeParametrization()
614 lsys->setParameter("matrix_reuse", "same_matrix"); in computeParametrization()
616 linearSystemCSRGmm<double> *lsys = new linearSystemCSRGmm<double>; in computeParametrization() local
618 linearSystemFull<double> *lsys = new linearSystemFull<double>; in computeParametrization() local
621 lsys->allocate(nodes.size()); in computeParametrization()
627 lsys->insertInSparsityPattern(nodeIndex[it->first.getVertex(i)], in computeParametrization()
664 lsys->addToMatrix(index0, index1, -c); in computeParametrization()
665 lsys->addToMatrix(index0, index0, c); in computeParametrization()
670 lsys->addToMatrix(row, row, 1); in computeParametrization()
673 lsys->zeroRightHandSide(); in computeParametrization()
676 lsys->addToRightHandSide(row, u[row]); in computeParametrization()
678 lsys->systemSolve(); in computeParametrization()
679 for(std::size_t i = 0; i < nodes.size(); i++) lsys->getFromSolution(i, u[i]); in computeParametrization()
681 lsys->zeroRightHandSide(); in computeParametrization()
684 lsys->addToRightHandSide(row, v[row]); in computeParametrization()
686 lsys->systemSolve(); in computeParametrization()
687 for(std::size_t i = 0; i < nodes.size(); i++) lsys->getFromSolution(i, v[i]); in computeParametrization()
689 delete lsys; in computeParametrization()