[Getdp] Question regarding the "MicroWaveElectricField" example on Wiki.
Christophe Geuzaine
cag32 at case.edu
Sat May 13 23:29:08 CEST 2006
Andrew Isaverdian wrote:
> In MWE.pro file, when I try
>
>
>
> Name MW_e_2D_PVP;
>
> System {
>
> { Name A; NameOfFormulation MW_e_2D; }
>
> }
>
> Operation {
>
> GenerateSeparate[A]; Lanczos[A, 19, {1:5}, 0]; SaveSolutions[A] ;
>
> }
>
>
>
> GetDP calculates only 2 first eigenvalues and stops at the third one
> with “BAD Eigenvalue…. Error = 10^99, etc.)
>
>
>
> If I try
>
>
>
> Name MW_e_2D_PVP;
>
> System {
>
> { Name A; NameOfFormulation MW_e_2D; Type ComplexValue; } // or
> Type Complex
>
> }
>
> Operation {
>
> GenerateSeparate[A]; EigenSolve[A, 5, 0, 0]; SaveSolutions[A] ;
>
> }
>
>
>
> then the solver stops with a message “ … Error: LinAlg_CopyMatrix not
> yet implemented.”
>
Indeed, it's not implemented with the default (Sparskit-based) solvers.
It should work if you compile getdp with the PETSc solvers ("./configure
--enable-petsc").
>
>
> Is there a way around it? For example, EigenSolve works fine with
> beam.pro, Elasticity_3D.pro.
>
>
>
> Thank you,
>
>
>
> Andrew
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
--
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine