[petsc-users] EVP Solutions with Subspace

[Sonya Blade]

Dear All,

I'd like to solve the problem of Generalized EVP where A is positive definite
and where B has products only on its diagonal and some of them might be zeros.

I started implementation with modifying the existing example(namely example 2
in the official Slepc site). Instead of using the Matrix viewer object I used 
reading A and B matrices directly from text file. Program produces the following 
error message, this has been triggered right after ierr = EPSSolve(eps);CHKERRQ(ierr);
call. I supply with main program body and required matrices stored in files.

[0]PETSC ERROR: --------------------- Error Message ----------------------------
--------
[0]PETSC ERROR: Object is in wrong state!
[0]PETSC ERROR: Not for unassembled matrix!
[0]PETSC ERROR: ----------------------------------------------------------------
--------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 6, Mon Feb 11 12:26:34 CST 20
13

I appreciated your help, any other example of solving GEVP with subspace, arnoldi 
or Lanczos methods is also wellcome.

Regards, 		 	   		  
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: main.c
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130405/f975736a/attachment-0001.c>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: IN.txt
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130405/f975736a/attachment-0002.txt>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: mass.txt
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130405/f975736a/attachment-0003.txt>