[petsc-users] lsqr usage

[Ildar Rakhmanov]

Here is full message:

Matrix A after construction
row 0: (0, 1)  (1, 0)
row 1: (0, 0)  (1, 1)
row 2: (0, 1)  (1, 1)
u
1
1
b
0
0
0
b
1
1
2
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Invalid argument!
[0]PETSC ERROR: Must be square matrix, rows 3 columns 2!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48
CDT 2011
[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: ./lsqr on a linux-gnu named ildar by ildar Wed May 23
17:13:44 2012
[0]PETSC ERROR: Libraries linked from /home/ildar/pbsm/lib
[0]PETSC ERROR: Configure run at Sat Mar  3 12:01:22 2012
[0]PETSC ERROR: Configure options --prefix=/home/ildar/pbsm
--download-hypre=1 --download-spooles=1 --download-plapack=1
--download-spai=1 --download-blacs=1 --download-triangle=1
--download-f-blas-lapack=1 --download-umfpack=1 --download-sowing=1
--download-c2html=1 --download-superlu_dist=1 --download-parmetis=1
--download-scalapack=1 --download-superlu=1 --download-mumps=1
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatGetOrdering() line 223 in src/mat/order/sorder.c
[0]PETSC ERROR: PCSetUp_ILU() line 194 in src/ksp/pc/impls/factor/ilu/ilu.c
[0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: main() line 109 in "unknowndirectory/"lsqr.c
application called MPI_Abort(MPI_COMM_WORLD, 62) - process 0




On Wed, May 23, 2012 at 5:29 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>  You cut and threw away the most important part of the error message.
> Where is the stack trace? Are we to guess where the error message was
> triggered? Please send the entire error message.
>
>   Barry
>
> On May 23, 2012, at 4:15 PM, Ildar Rakhmanov wrote:
>
> > Hi,
> > I seems I am missing something. I want to use KSPLSQR.
> > I set up lsqr.c (see attachment)
> > it uses 3x2 matrix.
> >
> > The results is
> >
> > ildar at ildar:~/lsqr> ./lsqr
> > Matrix A after construction
> > row 0: (0, 1)  (1, 0)
> > row 1: (0, 0)  (1, 1)
> > row 2: (0, 1)  (1, 1)
> > u
> > 1
> > 1
> > b
> > 0
> > 0
> > 0
> > b
> > 1
> > 1
> > 2
> > [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> > [0]PETSC ERROR: Invalid argument!
> > [0]PETSC ERROR: Must be square matrix, rows 3 columns 2!
> > [0]PETSC ERROR:
> ------------------------------------------------------------------------
> >
> >
> >
> > Could any one point out how to correctly use LSQR solver?
> >
> > <lsqr.c>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120523/c26b5342/attachment.html>