simpegEM1D icon indicating copy to clipboard operation
simpegEM1D copied to clipboard

Difficulties to handle with waveform

Open jc-barreto opened this issue 5 years ago • 0 comments

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?

jc-barreto avatar Nov 11 '20 16:11 jc-barreto