search

GeoMop / Genie

ERT 3D inversion tool
GNU General Public License v3.0
1 stars 0 forks source link

pg.geometricFactors() počítá geometrický faktor pro nás nevhodně #22

Closed radeksrb closed 3 years ago

radeksrb commented 3 years ago

Funkce https://github.com/gimli-org/gimli/blob/master/core/src/bert/bertMisc.cpp#L131 počíta geometrický faktor medou "3D mirrorspace", pro naše účely by byla lepší medota "3D fullspace".

Navrhované řešeni: Vytvořit vlastní funkci, která bude počítat tak jak chceme.

Touto chybou mohou být dotčeny výpočty 10A, 10B, 11A, 22A, 24A.

jbrezmorf commented 3 years ago

Jaká je motivace toto řešit?

Full space a mirror space by se měl lišit přenásobením 0.5. Jelikož by k tomu došlo pro všechna měření nemělo by to mít na inverzi vliv. Takže pokud tu není nějaká indicie toho, že by to vliv mít mělo, tak bych to neřešil.

radeksrb commented 3 years ago

Ve výsledku se požije vzorec https://github.com/gimli-org/gimli/blob/master/core/src/bert/bertMisc.cpp#L214 ovšem my pořebujeme https://github.com/gimli-org/gimli/blob/master/core/src/bert/bertMisc.cpp#L202 , který se nepožije protože funkce exactDCSolution je vždy volána s parametrem surfaceZ=0.0. Dle vzorců je patrné že nejde pouze o dvojnásobnou hodnotu, ale výsledek u "mirror space" není invariantní vůči posunutí ani natočení.

jbrezmorf commented 3 years ago

Jasný. Zdá se tedy, že exactDCSolution má podporu pro full 3D, ale 'geometricFactors' to nepoužívá. Bude tedy lépe si 'k' počítat sami podle potenciálu: '1.0 / (4.0 PI r)' umístěného v jednotlivých elektrodách.

radeksrb commented 3 years ago

Hotovo.