mpspack
mpspack copied to clipboard
2D Helmholtz scattering and eigenvalue problems via particular solutions and integral equations
use the correct eigenvector in evp.gsvdtension to output modes. Use 'ms' type search for gsvd. A lot of cleaning up to do here. This includes adding Dirichlet tension methods to...
``` take the single airdom and use its scattered field rep to compute far field at a given single angle, or list of angles. Make a user-friendly scattering.showfarfield which automates...
periodic case: spectral diff non-periodic segment: high-order interpolation. accuracy will suffer a bit, but saves ppl a bunch of annoyance w/ Z, Zp, Zpp for Kress.
``` No way to supress the output figure of [uj,x,y]=p.showmode() I want to take only the data. In addition to that the matrix uj is not a square matrix. How...
``` Since 'fd' is the premier method for low eigenvalues, should identify degeneracies closer than sqrt(eps) (?) and handle separately using local search ('fd' or 'ms'?) to try and regain...
``` Running examples/twodieldiscspointsrcscat.m but with the point source lying inside a dielectric fails because u_inc is only nonzero in the airdom(s). Need to check for this and change u_inc to...
``` Docs need: review of various options, eg bb, dx in all the methods calling problem.gridsolution EVP examples in Tutorial. ``` Original issue reported on code.google.com by `[email protected]` on 21...
``` in evp make a method for evaluating phi_j on the bdry then using its L2 norm to return absolute error in k eigenfrequency. Only do for Dirichlet BC for...
``` 'ms' mode: if degeneracy detected via multiple small singvals in gridminfit, should also make solvemodescoeffs handle this first in iter=0 mode (easy: SVD), and then in iter=1 mode (need...
``` The Laplace cases need to be added here, so Alpert can work. ``` Original issue reported on code.google.com by `[email protected]` on 2 May 2011 at 9:49