mvillemarette / mpspack

Automatically exported from code.google.com/p/mpspack
Other
0 stars 0 forks source link

I cannot run the function showmodes(problem,options) #33

Open GoogleCodeExporter opened 8 years ago

GoogleCodeExporter commented 8 years ago
s = segment.smoothstar(900, 0.3, 5);
d = domain(s, 1); s.setbc(-1, 'D'); p = evp(d);
o.eps = 0.1; o.modes = 1; o.khat = 'r'; o.fhat = 's';
p.solvespectrum([1000 1000.1], 'ntd', o);
o = []; o.inds = 1; o.dx = 0.002; o.fmm = 1; o.col = 'bw';showmodes(p,o);

This is the code from Alex Barnett and Andrew Hassel's article "Fast 
Computation of the Dirichlet eigenmodes via the spectral flow of the interior 
Neumann-to-Dirichlet map".

I get the following error in Matlab

What version of the product are you using? On what operating system?

Please provide any additional information below.

??? Method 'hfmm2dparttarg' is not defined for class 'utils' or is removed from
MATLAB's search path.

Error in ==> layerpot.layerpot>layerpot.evalFMM at 338
        U = utils.hfmm2dparttarg(iprec,k,N,node,ifslp,b.a(1)*charge,...

Error in ==> problem.problem>problem.pointsolution at 334
              u(ii) = u(ii) + b.evalFMM(co, p.x(ii), opts); % opts needed, dom!

Error in ==> problem.problem>problem.gridsolution at 376
      [u di] = pr.pointsolution(pointset(zz(:)), o); % make zz a col vec

Error in ==> evp.evp>evp.showmodes at 571
        [u gx gy di] = pe.gridsolution(o); % either conventional or GRF evalu

Error in ==> trifoil_excited_state_mspack at 5
o = []; o.inds = 1; o.dx = 0.002; o.fmm = 1; o.col = 'bw';showmodes(p,o);

Original issue reported on code.google.com by sijothan...@gmail.com on 9 Oct 2014 at 10:46

GoogleCodeExporter commented 8 years ago
[deleted comment]
GoogleCodeExporter commented 8 years ago
[deleted comment]
GoogleCodeExporter commented 8 years ago
[deleted comment]
GoogleCodeExporter commented 8 years ago
[deleted comment]
GoogleCodeExporter commented 8 years ago
[deleted comment]
GoogleCodeExporter commented 8 years ago
I have solved the issue by installing FMM2D fortran package and followed the 
mspack manual instructions to link the fortran file location. Now the problem 
is solved and I can plot the high frequency modes. I hope this may help others.

Original comment by sijothan...@gmail.com on 10 Oct 2014 at 11:09