[petsc-users] having issues with nullspace

Thu Mar 6 18:08:06 CST 2014

Mat,

here's sor. svd still running!

  0 KSP preconditioned resid norm 2.481182137607e+00 true resid norm
1.027599467617e+00 ||r(i)||/||b|| 1.000000000000e+00

  1 KSP preconditioned resid norm 1.428369058648e+00 true resid norm
4.775799486483e-01 ||r(i)||/||b|| 4.647530129184e-01

  2 KSP preconditioned resid norm 4.540176572297e-01 true resid norm
1.868438043300e-01 ||r(i)||/||b|| 1.818255168654e-01

  3 KSP preconditioned resid norm 1.999321637574e-01 true resid norm
1.185660060102e-01 ||r(i)||/||b|| 1.153815370158e-01

  4 KSP preconditioned resid norm 1.134340150403e-01 true resid norm
8.909919316428e-02 ||r(i)||/||b|| 8.670614959631e-02

  5 KSP preconditioned resid norm 7.566523445491e-02 true resid norm
7.766809539565e-02 ||r(i)||/||b|| 7.558207048878e-02

  6 KSP preconditioned resid norm 5.537203663154e-02 true resid norm
7.485665618568e-02 ||r(i)||/||b|| 7.284614146334e-02

  7 KSP preconditioned resid norm 4.390492285958e-02 true resid norm
7.302891311152e-02 ||r(i)||/||b|| 7.106748827039e-02

  8 KSP preconditioned resid norm 3.621562502013e-02 true resid norm
7.279539496786e-02 ||r(i)||/||b|| 7.084024200273e-02

  9 KSP preconditioned resid norm 3.309171702542e-02 true resid norm
7.139006972953e-02 ||r(i)||/||b|| 6.947266126469e-02

 10 KSP preconditioned resid norm 3.651904851324e-02 true resid norm
7.153789926206e-02 ||r(i)||/||b|| 6.961652036271e-02

 11 KSP preconditioned resid norm 5.973622761183e-02 true resid norm
7.347037429585e-02 ||r(i)||/||b|| 7.149709260383e-02

 12 KSP preconditioned resid norm 9.586867231154e-02 true resid norm
7.608732859604e-02 ||r(i)||/||b|| 7.404376023325e-02

 13 KSP preconditioned resid norm 1.635382207023e-02 true resid norm
7.118213033182e-02 ||r(i)||/||b|| 6.927030674402e-02

 14 KSP preconditioned resid norm 9.385509745433e-03 true resid norm
7.101728347294e-02 ||r(i)||/||b|| 6.910988737434e-02

 15 KSP preconditioned resid norm 7.423535884473e-03 true resid norm
7.104945342327e-02 ||r(i)||/||b|| 6.914119329785e-02

 16 KSP preconditioned resid norm 6.323417551714e-03 true resid norm
7.104173274422e-02 ||r(i)||/||b|| 6.913367998231e-02

 17 KSP preconditioned resid norm 5.970492489076e-03 true resid norm
7.099818043932e-02 ||r(i)||/||b|| 6.909129741373e-02

 18 KSP preconditioned resid norm 6.777117918864e-03 true resid norm
7.098851590105e-02 ||r(i)||/||b|| 6.908189244753e-02

 19 KSP preconditioned resid norm 1.556240130595e-02 true resid norm
7.103123679068e-02 ||r(i)||/||b|| 6.912346593114e-02

 20 KSP preconditioned resid norm 4.049582768289e-03 true resid norm
7.120382212555e-02 ||r(i)||/||b|| 6.929141593532e-02

 21 KSP preconditioned resid norm 6.862248582464e-04 true resid norm
7.111495778110e-02 ||r(i)||/||b|| 6.920493832682e-02

 22 KSP preconditioned resid norm 3.280255317757e-04 true resid norm
7.109751049924e-02 ||r(i)||/||b|| 6.918795964747e-02

 23 KSP preconditioned resid norm 1.943287835633e-04 true resid norm
7.109218692978e-02 ||r(i)||/||b|| 6.918277905948e-02

 24 KSP preconditioned resid norm 1.265015512889e-04 true resid norm
7.108929821126e-02 ||r(i)||/||b|| 6.917996792673e-02

 25 KSP preconditioned resid norm 1.264202759318e-04 true resid norm
7.108882409509e-02 ||r(i)||/||b|| 6.917950654446e-02

 26 KSP preconditioned resid norm 1.941274487871e-04 true resid norm
7.108881499846e-02 ||r(i)||/||b|| 6.917949769215e-02

 27 KSP preconditioned resid norm 1.237782236805e-04 true resid norm
7.108539247506e-02 ||r(i)||/||b|| 6.917616709155e-02

 28 KSP preconditioned resid norm 2.567842373425e-05 true resid norm
7.108693575066e-02 ||r(i)||/||b|| 6.917766891756e-02

 29 KSP preconditioned resid norm 1.624250348737e-05 true resid norm
7.108694548270e-02 ||r(i)||/||b|| 6.917767838821e-02

 30 KSP preconditioned resid norm 1.400709805305e-05 true resid norm
7.108704978977e-02 ||r(i)||/||b|| 6.917777989378e-02

 31 KSP preconditioned resid norm 1.258524583639e-05 true resid norm
7.108709377947e-02 ||r(i)||/||b|| 6.917782270200e-02

 32 KSP preconditioned resid norm 1.214563778665e-05 true resid norm
7.108699539645e-02 ||r(i)||/||b|| 6.917772696137e-02

 33 KSP preconditioned resid norm 1.175376138755e-05 true resid norm
7.108687362476e-02 ||r(i)||/||b|| 6.917760846025e-02

 34 KSP preconditioned resid norm 1.218910479225e-05 true resid norm
7.108702667143e-02 ||r(i)||/||b|| 6.917775739636e-02

 35 KSP preconditioned resid norm 1.786145773473e-05 true resid norm
7.108711007365e-02 ||r(i)||/||b|| 6.917783855855e-02

 36 KSP preconditioned resid norm 2.293730220946e-05 true resid norm
7.108631479962e-02 ||r(i)||/||b|| 6.917706464414e-02

 37 KSP preconditioned resid norm 5.223812358404e-06 true resid norm
7.108670686965e-02 ||r(i)||/||b|| 6.917744618388e-02

 38 KSP preconditioned resid norm 3.754603897641e-06 true resid norm
7.108673794435e-02 ||r(i)||/||b|| 6.917747642397e-02

 39 KSP preconditioned resid norm 3.359131707993e-06 true resid norm
7.108673993370e-02 ||r(i)||/||b|| 6.917747835988e-02

 40 KSP preconditioned resid norm 3.150846888622e-06 true resid norm
7.108674979713e-02 ||r(i)||/||b|| 6.917748795841e-02

 41 KSP preconditioned resid norm 2.872188456387e-06 true resid norm
7.108671888049e-02 ||r(i)||/||b|| 6.917745787212e-02

 42 KSP preconditioned resid norm 2.644889812492e-06 true resid norm
7.108676255906e-02 ||r(i)||/||b|| 6.917750037757e-02

 43 KSP preconditioned resid norm 2.702539684936e-06 true resid norm
7.108676352647e-02 ||r(i)||/||b|| 6.917750131900e-02

 44 KSP preconditioned resid norm 4.022429227783e-06 true resid norm
7.108676431469e-02 ||r(i)||/||b|| 6.917750208604e-02

 45 KSP preconditioned resid norm 1.579236965984e-06 true resid norm
7.108675412426e-02 ||r(i)||/||b|| 6.917749216931e-02

 46 KSP preconditioned resid norm 1.831765135201e-06 true resid norm
7.108674882093e-02 ||r(i)||/||b|| 6.917748700842e-02

 47 KSP preconditioned resid norm 1.846944675698e-06 true resid norm
7.108675619194e-02 ||r(i)||/||b|| 6.917749418146e-02

 48 KSP preconditioned resid norm 1.820006170791e-06 true resid norm
7.108675215517e-02 ||r(i)||/||b|| 6.917749025311e-02

 49 KSP preconditioned resid norm 1.770579229536e-06 true resid norm
7.108675662247e-02 ||r(i)||/||b|| 6.917749460042e-02

 50 KSP preconditioned resid norm 1.715571966759e-06 true resid norm
7.108675399614e-02 ||r(i)||/||b|| 6.917749204463e-02

 51 KSP preconditioned resid norm 1.646182320252e-06 true resid norm
7.108675922712e-02 ||r(i)||/||b|| 6.917749713512e-02

 52 KSP preconditioned resid norm 1.483705350928e-06 true resid norm
7.108675882398e-02 ||r(i)||/||b|| 6.917749674280e-02

 53 KSP preconditioned resid norm 1.116320649371e-06 true resid norm
7.108676123954e-02 ||r(i)||/||b|| 6.917749909349e-02

 54 KSP preconditioned resid norm 1.657360936891e-07 true resid norm
7.108676285301e-02 ||r(i)||/||b|| 6.917750066362e-02

 55 KSP preconditioned resid norm 9.060420145627e-07 true resid norm
7.108677215316e-02 ||r(i)||/||b|| 6.917750971399e-02

 56 KSP preconditioned resid norm 8.701510106194e-07 true resid norm
7.108675863210e-02 ||r(i)||/||b|| 6.917749655608e-02

 57 KSP preconditioned resid norm 6.250673067310e-07 true resid norm
7.108675858904e-02 ||r(i)||/||b|| 6.917749651417e-02

 58 KSP preconditioned resid norm 8.606839230868e-07 true resid norm
7.108676660388e-02 ||r(i)||/||b|| 6.917750431376e-02

 59 KSP preconditioned resid norm 2.619366455426e-07 true resid norm
7.108676234484e-02 ||r(i)||/||b|| 6.917750016910e-02

 60 KSP preconditioned resid norm 4.200943714852e-08 true resid norm
7.108676245753e-02 ||r(i)||/||b|| 6.917750027877e-02

 61 KSP preconditioned resid norm 7.666291884631e-09 true resid norm
7.108676261926e-02 ||r(i)||/||b|| 6.917750043615e-02

 62 KSP preconditioned resid norm 6.267989202765e-09 true resid norm
7.108676267362e-02 ||r(i)||/||b|| 6.917750048905e-02

 63 KSP preconditioned resid norm 6.322065776573e-09 true resid norm
7.108676266915e-02 ||r(i)||/||b|| 6.917750048470e-02

 64 KSP preconditioned resid norm 5.419388728381e-09 true resid norm
7.108676275176e-02 ||r(i)||/||b|| 6.917750056510e-02

 65 KSP preconditioned resid norm 4.864432437470e-09 true resid norm
7.108676262299e-02 ||r(i)||/||b|| 6.917750043978e-02

 66 KSP preconditioned resid norm 5.889407294927e-09 true resid norm
7.108676258555e-02 ||r(i)||/||b|| 6.917750040335e-02

 67 KSP preconditioned resid norm 2.329377620253e-09 true resid norm
7.108676262702e-02 ||r(i)||/||b|| 6.917750044371e-02

 68 KSP preconditioned resid norm 6.789865008325e-10 true resid norm
7.108676263968e-02 ||r(i)||/||b|| 6.917750045602e-02

 69 KSP preconditioned resid norm 9.951873144471e-10 true resid norm
7.108676262735e-02 ||r(i)||/||b|| 6.917750044402e-02

 70 KSP preconditioned resid norm 3.345568224138e-10 true resid norm
7.108676262789e-02 ||r(i)||/||b|| 6.917750044455e-02

 71 KSP preconditioned resid norm 2.723475696116e-10 true resid norm
7.108676263029e-02 ||r(i)||/||b|| 6.917750044688e-02

 72 KSP preconditioned resid norm 2.270307705471e-10 true resid norm
7.108676263411e-02 ||r(i)||/||b|| 6.917750045060e-02

 73 KSP preconditioned resid norm 3.952312974708e-10 true resid norm
7.108676262930e-02 ||r(i)||/||b|| 6.917750044592e-02

 74 KSP preconditioned resid norm 2.215148179871e-09 true resid norm
7.108676260334e-02 ||r(i)||/||b|| 6.917750042066e-02

 75 KSP preconditioned resid norm 6.854898442912e-11 true resid norm
7.108676262831e-02 ||r(i)||/||b|| 6.917750044496e-02

 76 KSP preconditioned resid norm 8.140088634442e-12 true resid norm
7.108676262834e-02 ||r(i)||/||b|| 6.917750044499e-02

 77 KSP preconditioned resid norm 2.214472895697e-12 true resid norm
7.108676262855e-02 ||r(i)||/||b|| 6.917750044519e-02

Linear solve converged due to CONVERGED_RTOL iterations 77

KSP Object: 1 MPI processes

  type: bcgs

  maximum iterations=10000, initial guess is zero

  tolerances:  relative=1e-12, absolute=1e-50, divergence=10000

  left preconditioning

  has attached null space

  using PRECONDITIONED norm type for convergence test

PC Object: 1 MPI processes

  type: sor

    SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1

  linear system matrix = precond matrix:

  Mat Object:   1 MPI processes

    type: seqaij

    rows=4225, cols=4225

    total: nonzeros=19301, allocated nonzeros=20273

    total number of mallocs used during MatSetValues calls =0

      has attached null space

      not using I-node routines

On Thu, Mar 6, 2014 at 4:06 PM, Mohammad Mirzadeh <mirzadeh at gmail.com>wrote:

> hummm just tried it in Matlab and you are correct -- they are different.
> What does this mean for my system? Also what is the correct approach here?
>
>
> On Thu, Mar 6, 2014 at 4:02 PM, Jed Brown <jed at jedbrown.org> wrote:
>
>> Mohammad Mirzadeh <mirzadeh at gmail.com> writes:
>>
>> > Jed,
>> >
>> > No the matrix is actually non-symmetric due to grid adaptivity (hanging
>> > nodes of QuadTree). Anyway, what do you exactly mean by the system being
>> > inconsistent?
>>
>> Sounds like the right and left null spaces are different.  You can test
>> by checking the null space using the transpose of your system.
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140306/99274854/attachment-0001.html>