QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Wrong radial part of longitudinal waves -> BUGS file 

Former-commit-id: 6144e79efefcb54b43e85b1b3c14ef9d907bb336
This commit is contained in:
Marek Nečada 2018-02-04 22:20:38 +00:00
parent 4dc151d5b9
commit b51b1dc2b5
1 changed files with 18 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
qpms/BUGS.rst
View File

@ -4,3 +4,21 @@ abort při určitých vstupech, např. -6 10 8 10
(ačkoliv fortran originál paří bez problémů)
vswf.c
======
Rozklad a znovusložení rovinné vlny vrací nesprávné výsledky pro radiální část podélných
vln, viz test_planewave_decomposition.c. Při reálné amplitudě (E) rovnoběžné s vlnovým
vektorem (k) to vypadá, že reálná část „radiální“ (ve vztahu k počátku souřadnic) složky
je správně, ale imaginární část neodpovídá. Pakliže amplituda E je komplexní, chyba je
i v reálné části (pro čistě imaginární amplitudu je chyba *právě* v reálné části).
Chybu jsem hledal především ve výpočtu kulových vln v bodě, tj. ve funkci
qpms_vswf_fill(), ale tam se mi ji zatím najít nepodařilo. Pochopitelně
také může být blbě přímo rozklad, tj. qpms_planewave2vswf_fill_sph().
Výše uvedené (chyba jen v imaginární části pro reál. amplitudu) snad časem napoví.
Neradiální část (podélných vln) je snad v pořádku. Rozklad čistě příčných vln (k.E = 0)
funguje dobře, což pro elektromagnetismus přirozeně stačí, ale stejně....
Alternativní funkce qpms_vswf_fill_alternative() *nevrací* stejné výsledky
jako qpms_vswf_fill, oproti qpms_vswf_fill() je nejspíš ještě chybovější.