VecView behaviour

Andreas Grassl Andreas.Grassl at student.uibk.ac.at
Wed Jun 3 17:29:22 CDT 2009


Barry Smith schrieb:
> 
>   When properly running nn-cg (are you sure everything is symmetric?)
> should require 10-30 iterations (certainly for model problems)

ok, this was the number I expected.

> 
>> nn-cg on 2 nodes 229 iterations, condition 6285
>> nn-cg on 4 nodes 331 iterations, condition 13312
> 
>   Are you sure that your operator has the null space of only constants?

no, I didn't touch anything regarding the null space since I thought it would be
done inside the NN-preconditioner. Does this mean I have to set up a null space
of the size of the Schur complement system, i.e. the number of interface DOF's?

cheers,

ando

> 
>     Barry
> 
> 
> On Jun 3, 2009, at 5:29 AM, Andreas Grassl wrote:
> 
>> Barry Smith schrieb:
>>>   Hmm, it sounds like the difference between local "ghosted" vectors
>>> and the global parallel vectors. But I do not understand why any of the
>>> local vector entries would be zero.
>>> Doesn't the vector X that is passed into KSP (or SNES) have the global
>>> entries and uniquely define the solution? Why is viewing that not right?
>>>
>>
>> I still don't understand fully the underlying processes of the whole PCNN
>> solution procedure, but trying around I substituted
>>
>> MatCreateIS(commw, ind_length, ind_length, PETSC_DECIDE, PETSC_DECIDE,
>> gridmapping, &A);
>>
>> by
>>
>> MatCreateIS(commw, PETSC_DECIDE, PETSC_DECIDE, actdof, actdof,
>> gridmapping, &A);
>>
>> and received the needed results.
>>
>> Furthermore it seems, that the load balance is now better, although I
>> still
>> don't reach the expected values, e.g.
>> ilu-cg 320 iterations, condition 4601
>> cg only 1662 iterations, condition 84919
>>
>> nn-cg on 2 nodes 229 iterations, condition 6285
>> nn-cg on 4 nodes 331 iterations, condition 13312
>>
>> or is it not to expect, that nn-cg is faster than ilu-cg?
>>
>> cheers,
>>
>> ando
>>
>> -- 
>> /"\                               Grassl Andreas
>> \ /    ASCII Ribbon Campaign      Uni Innsbruck Institut f. Mathematik
>>  X      against HTML email        Technikerstr. 13 Zi 709
>> / \                               +43 (0)512 507 6091
> 

-- 
 /"\
 \ /      ASCII Ribbon
  X    against HTML email
 / \


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 315 bytes
Desc: OpenPGP digital signature
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090604/a2b62ef8/attachment.pgp>


More information about the petsc-users mailing list