def __init__(self, lhs, solution, rhs, bcs=None):
self.solution = solution
self._init_lhs(lhs)
self._init_rhs(rhs)
self._apply_bcs(bcs)
preserve_solution_attributes(self.lhs, self.solution, self.rhs)
self.monitor = None
def __init__(self, residual_eval, solution, bcs, jacobian_eval):
self.residual_eval_callback = residual_eval
self.solution = solution
self.jacobian_eval_callback = jacobian_eval
# Preserve solution auxiliary attributes
self.residual_vector = self.residual_eval(solution)
self.jacobian_matrix = self.jacobian_eval(solution)
preserve_solution_attributes(self.jacobian_matrix, self.solution, self.residual_vector)
# Initialize BCs
self._init_bcs(bcs)
def __init__(self, lhs, solution, rhs, bcs=None):
self.lhs = lhs
self.solution = solution
self.rhs = rhs
self._apply_bcs(bcs)
preserve_solution_attributes(self.lhs, self.solution, self.rhs)
