Thanks a lot for your explanation, Barry,
This makes sense!
Fande,
On Fri, Feb 24, 2017 at 1:56 PM, Barry Smith wrote:
>
> Fande,
>
> Yes. Say one is doing a timestepping and using a direct solver. With
> time stepping we do not want to necessarily stop (in fact we
Fande,
Yes. Say one is doing a timestepping and using a direct solver. With time
stepping we do not want to necessarily stop (in fact we almost never) on a
failed solve (due to for example a failed factorization). We want the code to
continue up the stack until it gets to the time
Hi All,
In MatSolve(), there is a piece of code:
* if (mat->errortype) {ierr = PetscInfo1(mat,"MatFactorError
%D\n",mat->errortype);CHKERRQ(ierr);ierr = VecSetInf(x);CHKERRQ(ierr);
} else {ierr = (*mat->ops->solve)(mat,b,x);CHKERRQ(ierr); }*
If a direct solver such as LU or