qpms/BUGS.rst

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