Roy Stogner writes:
>> Interesting. Can you include the full unpreconditioned output? And
>> please compare with and without -ksp_gmres_modifiedgramschmidt.
>
> From "./main-opt -ksp_monitor_true_residual -ksp_norm_type unpreconditioned":
>0 KSP unpreconditioned resid norm 1.137340546775e-01
On Thu, 8 Jan 2015, Jed Brown wrote:
> Roy Stogner writes:
>
>> By "neither seems to do much"?
>>
>> I mean that if I set either "-ksp_pc_side right" or "-ksp_norm_type
>> unpreconditioned", then instead of ending up falsely "converged" at
>> 10 KSP preconditioned resid norm 1.580464771307e-13
Roy Stogner writes:
> By "neither seems to do much"?
>
> I mean that if I set either "-ksp_pc_side right" or "-ksp_norm_type
> unpreconditioned", then instead of ending up falsely "converged" at
> 10 KSP preconditioned resid norm 1.580464771307e-13 true resid norm
> 3.084021824157e-04 ||r(i)||
On Thu, 8 Jan 2015, Jed Brown wrote:
> Roy Stogner writes:
>
>> On Thu, 8 Jan 2015, Jed Brown wrote:
>>
10 KSP preconditioned resid norm 1.580464771307e-13 true resid norm
3.084021824157e-04 ||r(i)||/||b|| 2.711608086867e-03
number of iterations to solve adjoint: 10 final resid
Roy Stogner writes:
> On Thu, 8 Jan 2015, Jed Brown wrote:
>
>>> 10 KSP preconditioned resid norm 1.580464771307e-13 true resid norm
>>> 3.084021824157e-04 ||r(i)||/||b|| 2.711608086867e-03
>>> number of iterations to solve adjoint: 10 final residual of adjoint solve:
>>> 1.58046e-13
>>
>> Th
On Thu, 8 Jan 2015, Jed Brown wrote:
>> 10 KSP preconditioned resid norm 1.580464771307e-13 true resid norm
>> 3.084021824157e-04 ||r(i)||/||b|| 2.711608086867e-03
>> number of iterations to solve adjoint: 10 final residual of adjoint solve:
>> 1.58046e-13
>
> This is usually caused by a sing
John Peterson writes:
> Are you definitely allowed to pass the same matrix for both "A" and
> the preconditioner when you call KSPSolveTranspose()?
Yes, at least with KSPs capable of transpose solves. Can you send
-ksp_view along with the output requested in my last email?
signature.asc
Descri
On Thu, Jan 8, 2015 at 1:31 PM, Roy Stogner wrote:
>
> Copying discussion from libmesh-users, both because it looks like a
> library- rather than a user-level problem and because I'm hoping one
> of our PETSc expert lurkers will chime in.
>
> Summary: running a libMesh adjoint_solve() on a particu
Roy Stogner writes:
> On Thu, 8 Jan 2015, Jed Brown wrote:
>
>> Can you provide output with -ksp_monitor_true_residual? Are you
>> comparing a preconditioned residual to an unpreconditioned residual?
>
> Yes, and yes:
>
>0 KSP preconditioned resid norm 5.893300515271e-01 true resid norm
> 1
On Thu, 8 Jan 2015, Jed Brown wrote:
> Can you provide output with -ksp_monitor_true_residual? Are you
> comparing a preconditioned residual to an unpreconditioned residual?
Yes, and yes:
0 KSP preconditioned resid norm 5.893300515271e-01 true resid norm
1.137340546775e-01 ||r(i)||/||b|| 1
Roy Stogner writes:
> Copying discussion from libmesh-users, both because it looks like a
> library- rather than a user-level problem and because I'm hoping one
> of our PETSc expert lurkers will chime in.
>
> Summary: running a libMesh adjoint_solve() on a particular coupled
> multiphysics syste
Copying discussion from libmesh-users, both because it looks like a
library- rather than a user-level problem and because I'm hoping one
of our PETSc expert lurkers will chime in.
Summary: running a libMesh adjoint_solve() on a particular coupled
multiphysics system on a small mesh claims to solv
12 matches
Mail list logo
