Skip to content
This repository has been archived by the owner on Jul 12, 2023. It is now read-only.

Difficulties to handle with waveform #48

Open
jc-barreto opened this issue Nov 11, 2020 · 0 comments
Open

Difficulties to handle with waveform #48

jc-barreto opened this issue Nov 11, 2020 · 0 comments

Comments

@jc-barreto
Copy link

Hello,

Due to a lack of documentation, I'm facing some difficulties in understanding the waveform

First, I have tried to reproduce the results from this example that Heagy went through in this video to compare with the results from SimpegEM1D but they are not similar.

Second, I have tried to understand why in some cases the time zero in the waveform is setting when the current is one, and in other situations, like in this example not. This is the point that makes me more confused

Finally, to achieve my goal that is to have a good forward that returns the Jacobian and the predicted data, I tried to compare the 3 waveforms: VTEM, TEMFAST, and trapezoidal.The real case I'm studying is a square loop performed with TEMFAST equipment and I don't understand why I can't predict well the data with the original waveform

Last question: why to extract the Jacobian I need to multiply by a vector?

J = []
for i in range(len(model)):
    ei = np.zeros_like(model)
    ei[i] = 1.0
    J.append(simulation.Jvec(model, ei, f=fields))

J = np.vstack(J)

Is this the Jacobian that is called several times during the inversion?

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant