def _jacobian_matrix_assemble(self,
jacobian_form: ParametrizedTensorFactory,
petsc_jacobian: PETSc.Mat):
self.jacobian_matrix = PETScMatrix(petsc_jacobian)
evaluate(jacobian_form, tensor=self.jacobian_matrix)
def _residual_vector_assemble(self,
residual_form: ParametrizedTensorFactory,
petsc_residual: PETSc.Vec):
self.residual_vector = PETScVector(petsc_residual)
evaluate(residual_form, tensor=self.residual_vector)
def _jacobian_matrix_assemble(self,
jacobian_form: ParametrizedTensorFactory):
return evaluate(jacobian_form)
def _residual_vector_assemble(self,
residual_form: ParametrizedTensorFactory):
return evaluate(residual_form)
def _init_rhs(self, rhs, bcs):
self.rhs = evaluate(rhs)
def _init_lhs(self, lhs, bcs):
self.lhs = evaluate(lhs)
def _assemble_if_form(mat: ParametrizedTensorFactory):
return evaluate(mat)
def _init_rhs(self, rhs: ParametrizedTensorFactory,
bcs: (list_of(DirichletBC), ProductOutputDirichletBC,
dict_of(str, list_of(DirichletBC)),
dict_of(str, ProductOutputDirichletBC), None)):
self.rhs = evaluate(rhs)
