search

guagua-pamcn / a-dda

Automatically exported from code.google.com/p/a-dda
0 stars 0 forks source link

Rayleigh limit #36

Open GoogleCodeExporter opened 8 years ago

GoogleCodeExporter commented 8 years ago 
Implement special Rayleigh mode, which enables some simplifications (like
using CG iterative solver and calculating only total polarizability
tensor). Use additional symmetries of the particle (rotation around y and x
axes by 90 degrees). See ideas of Budko & Samokhin (2005) for iterative
solvers. Also think, how to calculate the total static polarizability
tensor of the particle, using three different incident polarizations.

Original issue reported on code.google.com by yurkin on 24 Dec 2008 at 7:10

GoogleCodeExporter commented 8 years ago 
The general overrelaxation method proposed by Budko & Samokhin (2005) should 
not be better than the default QMR.

However there is another relevant idea, which is described by issue 121.

Original comment by yurkin on 6 Jan 2011 at 7:34

GoogleCodeExporter commented 8 years ago 
Currently, numerical instability happens inside ADDA when decreasing the size 
of the scatterer. So the results are generally reliable down to d/lambda = 1e-5 
for all iterative solvers except CSYM. For the latter problems start to appear 
already at 1e-3 (reason unknown).

Original comment by yurkin on 29 Mar 2013 at 7:23