Skip to content
Snippets Groups Projects
Commit a7ad8d66 authored by Jan Thorbecke's avatar Jan Thorbecke
Browse files

solved bug in deconvolutiuon code

parent 67f11bde
No related branches found
No related tags found
No related merge requests found
......@@ -40,7 +40,7 @@ makewave fp=250e+6 dt=$dt file_out=wave.su nt=4096 t0=6e-9 verbose=1
file_ks=syncl_ro.su \
file_src=wave.su \
file_rcv=shot_zsrc1150.su \
src_type=7 \
src_type=1 \
src_orient=1 \
src_injectionrate=1 \
dtrcv=1e-10 \
......
......@@ -26,6 +26,7 @@ PRG = raytime
SRCC = $(PRG).c \
JespersRayTracer.c \
vidale.c \
Grid2Time1.c \
getParameters.c \
getWaveletInfo.c \
writeSrcRecPos.c \
......
#!/bin/bash
makewave nt=512 w=g2 file_out=wave.su w=g2 dt=0.004
makewave nt=512 file_out=wave.su w=g2 dt=0.004
green c=2000 nt=512 dt=0.004 zsrc1=500 xsrc1=0 file_src=wave.su file_out=green.su
fconv file_in1=green.su file_in2=wave.su mode=dec eps=0.0 reps=0.01 file_out=dec.su mode=dec verbose=1
../fconv file_in2=green.su file_in1=wave.su eps=0.01 reps=0.01 file_out=dec.su mode=dec verbose=2
......@@ -647,9 +647,11 @@ void deconv(float *data1, float *data2, float *decon, int nrec, int nsam,
qr = (float *) &cdec[0].r;
p1i = p1r + 1;
p2i = p2r + 1;
qi = qr + 1;
leps = reps*maxden+eps;
// fprintf(stderr,"eps=%e reps=%e max=%e => leps=%e\n", eps, reps, maxden, leps);
for (j = 0; j < n; j++) {
if (fabs(*p2r)>=fabs(*p2i)) {
*qr = (*p2r**p1r+*p2i**p1i)/(den[j]+leps);
*qi = (*p2r**p1i-*p2i**p1r)/(den[j]+leps);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment