[petsc-users] issue with different versions of petsc-dev

[Jed Brown]

On Fri, Aug 17, 2012 at 2:20 PM, Mohammad Mirzadeh <mirzadeh at gmail.com>wrote:

> You mean columns within a row? Nope.
>

Okay, they are supposed to be. You can loop over your rows calling
PetscSortInt() on each row separately.

You didn't get the error before because debugging wasn't on, but you were
also very lucky because many matrix kernels would have behaved incorrectly
with unsorted data.


>
>
> On Fri, Aug 17, 2012 at 12:18 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> On Fri, Aug 17, 2012 at 2:11 PM, Mohammad Mirzadeh <mirzadeh at gmail.com>wrote:
>>
>>> Hi guys,
>>>
>>> I just realize my code breaks with new versions of petsc-dev. I use
>>> MatCreateMPIAIJWithSplitArrays() funtions and it works with this verison:
>>>
>>> Petsc Development HG revision: b27fc9457bdc002c5a5de2b07797ae1a93094370
>>>  HG Date: Thu May 03 13:42:15 2012 -0500
>>>
>>> but breaks with:
>>>
>>> Petsc Development HG revision: f9c6cac2d69c724a2258d4e0dc2f51b0825aa875
>>>  HG Date: Thu Aug 16 08:37:18 2012 -0700
>>>
>>> Here's the error code:
>>> [0]PETSC ERROR: --------------------- Error Message
>>> ------------------------------------
>>> [0]PETSC ERROR: Argument out of range!
>>> [0]PETSC ERROR: Column entry number 1 (actual colum 2) in row 3 is not
>>> sorted!
>>>
>>
>> Are your rows sorted?
>>
>>
>>>  [0]PETSC ERROR:
>>> ------------------------------------------------------------------------
>>> [0]PETSC ERROR: Petsc Development HG revision:
>>> f9c6cac2d69c724a2258d4e0dc2f51b0825aa875  HG Date: Thu Aug 16 08:37:18 2012
>>> -0700
>>> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
>>> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
>>> [0]PETSC ERROR: See docs/index.html for manual pages.
>>> [0]PETSC ERROR:
>>> ------------------------------------------------------------------------
>>> [0]PETSC ERROR: ./FVM on a arch-linu named mohammad-desktop by mohammad
>>> Fri Aug 17 12:06:15 2012
>>> [0]PETSC ERROR: Libraries linked from
>>> /home/mohammad/soft/petsc-dev_new/petsc-dev/arch-linux2-cxx-debug/lib
>>> [0]PETSC ERROR: Configure run at Fri Aug 17 11:45:31 2012
>>> [0]PETSC ERROR: Configure options --with-clanguage=cxx
>>> --with-mpi-dir=/usr/bin --download-f-blas-lapack=1 --download-hypre=1
>>> --download-parmetis=1 --download-metis=1 --download-superlu_dist=1
>>> [0]PETSC ERROR:
>>> ------------------------------------------------------------------------
>>> [0]PETSC ERROR: MatCreateSeqAIJWithArrays() line 4268 in
>>> /home/mohammad/soft/petsc-dev_new/petsc-dev/src/mat/impls/aij/seq/aij.c
>>> [0]PETSC ERROR: MatCreateMPIAIJWithSplitArrays() line 5731 in
>>> /home/mohammad/soft/petsc-dev_new/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c
>>> [0]PETSC ERROR: setLinearSystemMatrix() line 84 in
>>> "unknowndirectory/"../../../CASL/lib/algebra/petscLinearSolver.cpp
>>>
>>> Is there anything changed regarding numbering of diagonal and
>>> off-diagonal entries when using the said function between these two
>>> versions?
>>>
>>> Thanks,
>>> Mohammad
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120817/b359af75/attachment.html>