On Sun, Apr 10, 2011 at 10:14 AM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im">Matt,<br>
<br>
>><br>
>> The reference counting is ASA must be wrong.<br>
><br>
> I have fixed this, but the output does not seem to match. Did this work<br>
> before?<br>
</div>I got same thing - runs, but gives different output.<br>
I saw the crash from this Friday's nightly test, which did not occur<br>
in the previous weeks.<br>
I do not not know if previously the output differ.<br>
<br>
The residual of new output is in the same order as output/ex40_1.out<br>
and is slight better. Do you want update output/ex40_1.out?</blockquote><div><br></div><div>I guess it should be updated. I do not know who added it.</div><div><br></div><div>Does anyone want to claim responsibility for KSP ex40?</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font color="#888888"><br>
Hong<br>
</font><div><div></div><div class="h5"><br>
><br>
>><br>
>>    Matt<br>
>><br>
>>><br>
>>> petsc-dev/src/ksp/ksp/examples/tutorials/ex40.c (runex40) crashes<br>
>>> 0]PETSC ERROR: --------------------- Error Message<br>
>>> ------------------------------------<br>
>>> [0]PETSC ERROR: Invalid argument!<br>
>>> [0]PETSC ERROR: Invalid object classid 0<br>
>>> This often happens if you compile with PETSC_USE_DYNAMIC_LIBRARIES,<br>
>>> but link with static libraries.!<br>
>>> ...<br>
>>> [0]PETSC ERROR: ClassRegLogGetClass() line 282 in<br>
>>> src/sys/plog/utils/classLog.c<br>
>>> [0]PETSC ERROR: PetscLogObjDestroyDefault() line 383 in<br>
>>> src/sys/plog/utils/classLog.c<br>
>>> [0]PETSC ERROR: DMDestroy() line 164 in src/dm/interface/dm.c<br>
>>> [0]PETSC ERROR: PCDestroyLevel_ASA() line 300 in<br>
>>> src/ksp/pc/impls/asa/asa.c<br>
>>> [0]PETSC ERROR: PCDestroy_ASA() line 1928 in src/ksp/pc/impls/asa/asa.c<br>
>>> [0]PETSC ERROR: PCDestroy() line 116 in src/ksp/pc/interface/precon.c<br>
>>> [0]PETSC ERROR: KSPDestroy_() line 757 in src/ksp/ksp/interface/itfunc.c<br>
>>> [0]PETSC ERROR: main() line 234 in src/ksp/ksp/examples/tutorials/ex40.c<br>
>>> application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<br>
>>><br>
>>> Could it result of your recent change?  Hong<br>
>>><br>
>>> changeset:   18673:75eb7f4d3ab7<br>
>>> user:        Matthew G. Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
>>> date:        Sat Apr 02 18:46:15 2011 -0500<br>
>>><br>
>>> description:<br>
>>> SNES ex10 now runs in serial<br>
>>> - Added MatSetValuesTopology()<br>
>>> - DMGetMatrix_Mesh() now sets the local to global mapping<br>
>>> - We compose the Mesh with global/local vectors<br>
>>> - Added DMCreateLocalToGlobalMapping() to Mesh<br>
>>> - Added matrix filling to DMMeshCreateMatrix()<br>
>>> - Fixed DMDestroy() for all subclasses<br>
>>> - Moved DMDASetMatrPreallocateOnly() to DMSetMatrixPreallocateOnly()<br>
>>>  and changed option to -dm_preallocate_only<br>
>>> - Fixed test output<br>
>><br>
>><br>
>><br>
>> --<br>
>> What most experimenters take for granted before they begin their<br>
>> experiments is infinitely more interesting than any results to which their<br>
>> experiments lead.<br>
>> -- Norbert Wiener<br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments<br>
> is infinitely more interesting than any results to which their experiments<br>
> lead.<br>
> -- Norbert Wiener<br>
><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>