David:
I'll patch petsc-maint (v3.7), then merge it to petsc-master.
It might take 1-2 days for regression tests.
Hong
On Tue, Sep 20, 2016 at 2:01 PM, Hong wrote:
>
>> David :
>> This is a bug in PETSc. A change
>> $ git diff ../../../pc/impls/factor/cholesky/cholesky.c
>>
On Tue, Sep 20, 2016 at 2:01 PM, Hong wrote:
> David :
> This is a bug in PETSc. A change
> $ git diff ../../../pc/impls/factor/cholesky/cholesky.c
> diff --git a/src/ksp/pc/impls/factor/cholesky/cholesky.c
> b/src/ksp/pc/impls/factor/cholesky/cholesky.c
> index
David :
This is a bug in PETSc. A change
$ git diff ../../../pc/impls/factor/cholesky/cholesky.c
diff --git a/src/ksp/pc/impls/factor/cholesky/cholesky.c
b/src/ksp/pc/impls/factor/cholesky/cholesky.c
index 953d551..cc28369 100644
--- a/src/ksp/pc/impls/factor/cholesky/cholesky.c
+++
On Mon, Sep 19, 2016 at 11:04 PM, Hong wrote:
> David :
> I did following:
>
> PC pc;
> Mat F;
> ierr = KSPGetPC(ksp,);CHKERRQ(ierr);
> ierr = PCReset(pc);CHKERRQ(ierr);
> ierr = KSPSetOperators(ksp,A,A);CHKERRQ(ierr);
> ierr =
On Tue, Sep 20, 2016 at 6:44 AM, Florian Lindner
wrote:
> Hello,
>
> to compute / approximate the condition number of a MATSBAIJ, I put
>
> -ksp_view
>
> # Conditon number estimate
> -pc_type none
> -ksp_type gmres
> -ksp_monitor_singular_value
> -ksp_gmres_restart 1000
>
>
Hello,
to compute / approximate the condition number of a MATSBAIJ, I put
-ksp_view
# Conditon number estimate
-pc_type none
-ksp_type gmres
-ksp_monitor_singular_value
-ksp_gmres_restart 1000
in my .petscrc
Output is like:
[...]
566 KSP Residual norm 1.241765807317e-07 % max