after the set functions.
PETSc needs to "assemble" the changes you made when calling set()
before you can do anything with the matrix.
--
Anders
> Now, when I'm trying to solve the problem with:
>
> SLEPcEigenvalueSolver esolver;
> esolver.solve(C);
>
> I get the following error message:
>
>
> [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 2.3.3, Patch 8, Fri Nov 16
> 17:03:40 CST 2007 HG revision:
414581156e67e55c761739b0deb119f7590d0f4b
> [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: ./mass on a linux-gnu named ... by witkowsk Tue Nov 27
> 16:17:51 2007
> [0]PETSC ERROR: Libraries linked from
> /u/witkowski/local/petsc-2.3.3-p8/lib/linux-gnu-cxx
> [0]PETSC ERROR: Configure run at Tue Nov 27 11:24:27 2007
> [0]PETSC ERROR: Configure options --with-mpi-dir=/u/witkowski/local
> --with-clanguage=cxx --download-hypre=yes --with-hypre=1
> --download-umfpack=yes --with-umfpack=1 --with-shared=0
> [0]PETSC ERROR:
>
------------------------------------------------------------------------
> [0]PETSC ERROR: MatMult() line 1616 in src/mat/interface/matrix.c
> [0]PETSC ERROR: STApply_Shift() line 30 in src/st/impls/shift/shift.c
> [0]PETSC ERROR: STApply() line 49 in src/st/interface/stsolve.c
> [0]PETSC ERROR: EPSBasicArnoldi() line 86 in
src/eps/impls/arnoldi/arnoldi.c
> [0]PETSC ERROR: EPSSolve_KRYLOVSCHUR() line 91 in
> src/eps/impls/krylovschur/krylovschur.c
> [0]PETSC ERROR: EPSSolve() line 71 in src/eps/interface/solve.c
> Eigenvalue solver (krylovschur) converged in 1 iterations.
>
>
> What does it meen? And what have I made wrong?
>
> Regards,
>
> Thomas
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
_______________________________________________
DOLFIN-dev mailing list
DOLFIN-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/dolfin-dev