[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