[petsc-users] matrix-free KSPLSQR and matrix norm

[Matthew Knepley]

On Mon, Jun 27, 2011 at 9:34 AM, Alexander Grayver
<agrayver at gfz-potsdam.de>wrote:

> Thanks Barry,
>
> Now I get this:


Please send the entire error message (this is missing the stack), but I am
guessing you
are not changed the default preconditioner (ILU) which cannot act on a shell
matrix. Use
-ksp_view to see what solver is begin used.

   Matt


> [0]PETSC ERROR: --------------------- Error Message
> ----------------------------
> --------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Mat type shell!
>
> [0]PETSC ERROR: ------------------------------**
> ------------------------------**----
> --------
> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 7, Mon Dec 20 14:26:37
> CST 20
> 10
> [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: test.exe on a cygwin
> -cx named MT0. by Unknown Mon Jun 27 16:30:42 2011
>
> [0]PETSC ERROR: Libraries linked from /cygdrive/d/dev/petsc-3.1-p7/**
> cygwin-cxx-de
> bug/lib
> [0]PETSC ERROR: Configure run at Fri Mar  4 10:10:06 2011
> [0]PETSC ERROR: Configure options --with-cc="win32fe cl" --with-fc="win32fe
> ifor
> t" --with-cxx="win32fe cl" --with-windows-graphics=0
> --download-f-blas-lapack=1
> --with-precision=double --with-scalar-type=complex --with-clanguage=cxx
> --with-m
> pi-include=/cygdrive/d/Dev/**HPC_2008_SDK/Include
> --with-mpi-lib="[/cygdrive/d/**Dev
> /HPC_2008_SDK/Lib/amd64/msmpi.**lib,/cygdrive/d/Dev/HPC_2008_**
> SDK/Lib/amd64/msmpife
> c.lib]" --useThreads=0 --useThreads=0
> [0]PETSC ERROR: ------------------------------**
> ------------------------------**----
> --------
> [0]PETSC ERROR: MatNorm() line 4512 in src/mat/interface/D:\dev\**
> PETSC-3.1-P\src\mat\INTERF~1\**matrix.c
>
>
> On 27.06.2011 16:32, Barry Smith wrote:
>
>>   You should pass H as third argument to KSPSetOperators() as well as the
>> second.
>>
>>   Barry
>>
>> On Jun 27, 2011, at 9:23 AM, Alexander Grayver wrote:
>>
>>  Hello!
>>>
>>> I'm using matrix-free formulation and KSPLSQR to solve normal equations
>>> system:
>>>
>>> (A'*A)b = -g
>>>
>>> I managed to write all necessary code:
>>>
>>>  call MatCreateShell(MPI_COMM_WORLD,**m,N,N,N,PETSC_NULL_INTEGER,H,**
>>> ierr)
>>>  call MatShellSetOperation(H,MATOP_**MULT,MFMatMult,ierr)
>>>
>>>  call KSPCreate(MPI_COMM_WORLD,ksp,**ierr);CHKERRQ(ierr)
>>>  call KSPSetType(ksp,KSPLSQR,ierr);**CHKERRQ(ierr)
>>>  call KSPGetPC(ksp,pc,ierr);CHKERRQ(**ierr)
>>>  call PCSetType(pc,PCNONE,ierr);**CHKERRQ(ierr)
>>>  call KSPSetOperators(ksp,H,PETSC_**NULL_OBJECT,SAME_**
>>> PRECONDITIONER,ierr);CHKERRQ(**ierr)
>>>  call KSPSolve(ksp,g,b,ierr);**CHKERRQ(ierr)
>>>
>>>
>>> but fail in KSPSolve stage with that:
>>> [0]PETSC ERROR: --------------------- Error Message
>>> ----------------------------
>>> --------
>>> [0]PETSC ERROR: Object is in wrong state!
>>> [0]PETSC ERROR: Mat object's type is not set: Argument # 1!
>>> [0]PETSC ERROR: ------------------------------**
>>> ------------------------------**----
>>> --------
>>> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 7, Mon Dec 20 14:26:37
>>> CST 20
>>> 10
>>> [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: test.exe on a cygwin
>>>
>>> -cx named MT0. by Unknown Mon Jun 27 16:18:43 2011
>>> [0]PETSC ERROR: Libraries linked from /cygdrive/d/dev/petsc-3.1-p7/**
>>> cygwin-cxx-de
>>> bug/lib
>>> [0]PETSC ERROR: Configure run at Fri Mar  4 10:10:06 2011
>>> [0]PETSC ERROR: Configure options --with-cc="win32fe cl"
>>> --with-fc="win32fe ifor
>>> t" --with-cxx="win32fe cl" --with-windows-graphics=0
>>> --download-f-blas-lapack=1
>>> --with-precision=double --with-scalar-type=complex --with-clanguage=cxx
>>> --with-m
>>> pi-include=/cygdrive/d/Dev/**HPC_2008_SDK/Include
>>> --with-mpi-lib="[/cygdrive/d/**Dev
>>> /HPC_2008_SDK/Lib/amd64/msmpi.**lib,/cygdrive/d/Dev/HPC_2008_**
>>> SDK/Lib/amd64/msmpife
>>> c.lib]" --useThreads=0 --useThreads=0
>>> [0]PETSC ERROR: ------------------------------**
>>> ------------------------------**----
>>> --------
>>> [0]PETSC ERROR: MatNorm() line 4507 in src/mat/interface/D:\dev\**
>>> PETSC-3.1-P\src\mat\INTERF~1\**matrix.c
>>> [0]PETSC ERROR: KSPSolve_LSQR() line 90 in src/ksp/ksp/impls/lsqr/D:\dev\
>>> **PETSC-3.1-P\src\ksp\ksp\impls\**lsqr\lsqr.c
>>>
>>> Any suggestions? Do I have to set MATOP_NORM operation? Then how?
>>>
>>> Regards,
>>> Alexander
>>>
>>
>


-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110627/ec0057e6/attachment.htm>