tiny changes

fdelmodc/demo/benchmark.scr

@@ -15,7 +15,8 @@ makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1
 
 xsrc=-200
 
-time ../fdelmodc \
+set -x
+../fdelmodc \
     file_cp=syncl_cp.su ischeme=1 \
     file_den=syncl_ro.su \
     file_src=wave.su \
@@ -35,5 +36,5 @@ time ../fdelmodc \
     zrcv1=0 zrcv2=0 \
     xsrc=0 zsrc=1150 \
     ntaper=200 \
-    left=4 right=4 top=1 bottom=4
+    left=2 right=2 top=1 bottom=2
 

marchenko/demo/oneD/README

@@ -151,30 +151,31 @@ mergeGreenRef.eps 	=> Figure 7
 
 To compute the snapshots 
 
-==> run backpropf2.scr This will take about 2 minutes. The generated output files are
+==> run backpropf2.scr This will take about 1 minute. The generated output file is
 	- backpropf2_sp.su
-	- backpropf2detail_sp.su
 
-The postscript files of Figure 8  are generated with 
+The postscript files of Figure 8 are generated with 
 
 ==> run epsBackprop.scr
 
 -- Figure 8 column 1
 backpropf2_-0.30.eps
 backpropf2_-0.15.eps
-backpropf2_0.0.eps
-backpropf2_0.15.eps
-backpropf2_0.30.eps
--- Figure 8 column 2
 backpropf2_-0.03.eps
 backpropf2_-0.02.eps
 backpropf2_0.00.eps
-backpropf2_0.02.eps
+-- Figure 8 column 2
+backpropf2_0.30.eps
+backpropf2_0.15.eps
 backpropf2_0.03.eps
+backpropf2_0.02.eps
+backpropf2_0.00.eps
 -- Figure 8 column 3
-backpropf2sum_0.03.eps
-backpropf2sum_0.15.eps
 backpropf2sum_0.30.eps
+backpropf2sum_0.15.eps
+backpropf2sum_0.03.eps
+backpropf2sum_0.02.eps
+backpropf2_0.00.eps
 
 
 The figures in the appendix, to explain the different options in the programs, are reproduced by
@@ -190,3 +191,4 @@ noise_above4.eps
 
 -- Figure A-2
 iniFocus_shifts.eps
+

marchenko/demo/oneD/backpropf2.scr

@@ -60,31 +60,3 @@ fdelmodc \
     left=2 right=2 top=2 bottom=2
 
 
-#for detailed, middle column, 
-#backpropagate f2.su and collect snapshots
-fdelmodc \
-    file_cp=$file_cp ischeme=1 iorder=4 \
-    file_den=$file_ro \
-    file_src=pplusdt.su \
-    file_rcv=backprop_f2_z900.su \
-    grid_dir=0 \
-    src_type=1 \
-    src_injectionrate=1 \
-	src_orient=1 \
-    rec_type_vz=0 \
-    rec_type_p=1 \
-    rec_int_vz=2 \
-    dtrcv=0.004 \
-	rec_delay=0.0 \
-    verbose=2 \
-    tmod=3.10 \
-    dxrcv=5.0 \
-    xrcv1=-2250 xrcv2=2250 \
-    zrcv1=900 zrcv2=900 \
-	zsrc=0 xsrc=0 \
-    npml=101 \
-	file_snap=backpropf2detail.su tsnap1=2.0245 dtsnap=0.0005 tsnap2=2.0645 dxsnap=2.5 dzsnap=2.5 zsnap1=0 zsnap2=1250 xsnap1=-1000 xsnap2=1000 \
-    sna_type_vz=0 \
-    sna_type_p=1 \
-    left=2 right=2 top=2 bottom=2
-

marchenko/demo/oneD/epsBackprop.scr

@@ -54,7 +54,7 @@ done
 #add pos and negative times to get response of homogenoeus Green's function
 
 file_snap="backpropf2"
-for fldr in 71 86 98 101;
+for fldr in 71 86 98 99 101;
 do
   times=$(echo "scale=2; -0.01*(${fldr}-101)" | bc -l)
   atime=`printf "%4.2f" $times`

marchenko/demo/oneD/epsCompare.scr

@@ -28,5 +28,5 @@ suwind < referenceP_rp.su j=50 s=1 | \
     n1tic=2 d2=250 f1=0.0 x1beg=0 x1end=2.004 d1num=0.4 \
     f2=-2250 f2num=-2000 d2num=500 > ref.eps
 
-psmerge in=green.eps in=ref.eps > mergeGreenRef.eps
+psmerge in=ref.eps in=green.eps > mergeGreenRef.eps
 

marchenko/fmute.c

@@ -38,7 +38,7 @@ char *sdoc[] = {
 " Required parameters: ",
 " ",
 "   file_mute= ................ input file with event that defines the mute line",
-"   file_shot= ................ input data thats is muted",
+"   file_shot= ................ input data that is muted",
 " ",
 " Optional parameters: ",
 " ",

marchenko/marchenko.c

@@ -405,14 +405,14 @@ int main (int argc, char **argv)
                 j = 0;
                 Ni[l*nxs*nts+i*nts+j]    = -iRN[l*nxs*nts+i*nts+j];
                 pmin[l*nxs*nts+i*nts+j] += iRN[l*nxs*nts+i*nts+j];
-                energyNi = sqrt(iRN[l*nxs*nts+i*nts+j]*iRN[l*nxs*nts+i*nts+j]);
+                energyNi = iRN[l*nxs*nts+i*nts+j]*iRN[l*nxs*nts+i*nts+j];
                 for (j = 1; j < nts; j++) {
                     Ni[l*nxs*nts+i*nts+j]    = -iRN[l*nxs*nts+i*nts+nts-j];
                     pmin[l*nxs*nts+i*nts+j] += iRN[l*nxs*nts+i*nts+j];
-                    energyNi = sqrt(iRN[l*nxs*nts+i*nts+j]*iRN[l*nxs*nts+i*nts+j]);
+                    energyNi += iRN[l*nxs*nts+i*nts+j]*iRN[l*nxs*nts+i*nts+j];
                 }
             }
-            if (verbose >=2) vmess("    - operator %d at iteration %d has energy %e", l, iter, energyNi);
+            if (verbose >=2) vmess("    - operator %d at iteration %d has energy %e", l, iter, sqrt(energyNi));
 
         }
 
@@ -474,7 +474,7 @@ int main (int argc, char **argv)
         }
         else {
             /* next iterations  */
-            /* N_k(x,t) = -N_(k-1)(x,-t) */
+            /*  -N_(k-1)(x,-t) */
 
             /* update f2 */
             for (l = 0; l < Nsyn; l++) {