convection

Author: sfomel

user/fomels/Mcnvdip.c

@@ -0,0 +1,67 @@
+/* Dip estimation from convection */
+/*
+  Copyright (C) 2025 University of Texas at Austin
+  
+  This program is free software; you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation; either version 2 of the License, or
+  (at your option) any later version.
+  
+  This program is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+  
+  You should have received a copy of the GNU General Public License
+  along with this program; if not, write to the Free Software
+  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+*/
+#include <math.h>
+
+#include <rsf.h>
+#include "predict.h"
+
+int main (int argc, char* argv[])
+{
+    int order, n1, n2, i1, i2, nw;
+    float **conv, *trace, eps;
+    sf_file cnv, dip;
+
+    sf_init(argc, argv);
+    cnv = sf_input("in");
+    dip = sf_output("out");
+
+    if (SF_FLOAT != sf_gettype(cnv)) sf_error("Need float input");
+    if (!sf_histint(cnv,"n1",&n1)) sf_error("Need n1= in input");
+    if (!sf_histint(cnv,"n2",&nw)) sf_error("Need n2= in input");    
+    if (!sf_histint(cnv,"n3",&n2)) sf_error("Need n3= in input");
+
+    order = nw/2;
+
+    sf_unshiftdim(cnv,dip,2);
+
+    if (!sf_getfloat("eps",&eps)) eps=0.01;
+    /* regularization */
+    
+    predict_init (n1, n2, eps*eps, order, 1, false);
+
+    conv = sf_floatalloc2(n1,nw);
+    trace = sf_floatalloc(n1);
+ 
+    for (i2=0; i2 < n2; i2++) {	
+	for (i1=0; i1 < n1; i1++) {
+	    trace[i1] = i1;
+	}
+	sf_floatread(conv[0],n1*nw,cnv);
+	predict_step(false,true,trace,conv);
+	for (i1=0; i1 < n1; i1++) {
+	    trace[i1] = i1 - trace[i1];
+	}
+	sf_floatwrite(trace,n1,dip);
+    }
+  
+    exit(0);
+}
+  
+
+  

user/fomels/Mcnvpaint.c

@@ -25,7 +25,7 @@ int main (int argc, char* argv[])
 {
   int order, n1, n2, n12, i1, i2, i0, nw;
     float **u, ***conv, *trace;
-    float o1, d1, o2, d2, eps, *time;
+    float o1, d1, eps, *time;
     sf_file cnv, out, seed;
 
     sf_init(argc, argv);
@@ -34,14 +34,10 @@ int main (int argc, char* argv[])
 
     if (SF_FLOAT != sf_gettype(cnv)) sf_error("Need float input");
     if (!sf_histint(cnv,"n1",&n1)) sf_error("Need n1= in input");
-    
+    if (!sf_histint(cnv,"n2",&nw)) sf_error("Need n2= in input");    
     if (!sf_histint(cnv,"n3",&n2)) sf_error("Need n3= in input");
-    if (!sf_histfloat(cnv,"o3",&o2)) o2=0.; 
-    if (!sf_histfloat(cnv,"d3",&d2)) d2=1.; 
 
     n12 = n1*n2;
-    
-    if (!sf_histint(cnv,"n2",&nw)) sf_error("Need n2= in input");
     order = nw/2;
 
     sf_unshiftdim(cnv,out,2);

user/fomels/SConstruct

@@ -6,14 +6,14 @@ progs = '''
 abalance analytical angle angle2 approx arrival bdix beamform1 bil1
 bilat2 blur boxcascade cbeamform1 causinv ccausint cchebyshevp
 cfftexpa-dev cfftwave1 cdivn cflow chain2dfft chebvc chebyshev
-chebyshevp clpf cltft cltftfft cnvd cnvpaint constperm constpermh
-constpermh1 convection cosftwave1 cpef cr ctf2dprec deblur diphase
-distance divn dix donut dpeiko eikonal eikonalvti eno2 erf exgr
-fedchain fedchain1 fedchain2 fedchain21 fft2 fftexp0 fftexp1 fftexp3
-fftexpa fftexp0a fftone ffttest fftwave1 fftwave2 fftwave3 findmin2
-focus fpow freqest gaussmooth gbeamform imray interf interp2 interpt
-iphase kdsort kdtree kron label legacy lfftexp0 llpf localskew locov
-lpf lrmig0 lsfit max2 median mffit mig3 miss3 morph nconv nnint
+chebyshevp clpf cltft cltftfft cnvd cnvdip cnvpaint constperm
+constpermh constpermh1 convection cosftwave1 cpef cr ctf2dprec deblur
+diphase distance divn dix donut dpeiko eikonal eikonalvti eno2 erf
+exgr fedchain fedchain1 fedchain2 fedchain21 fft2 fftexp0 fftexp1
+fftexp3 fftexpa fftexp0a fftone ffttest fftwave1 fftwave2 fftwave3
+findmin2 focus fpow freqest gaussmooth gbeamform imray interf interp2
+interpt iphase kdsort kdtree kron label legacy lfftexp0 llpf localskew
+locov lpf lrmig0 lsfit max2 median mffit mig3 miss3 morph nconv nnint
 nnshape nnshapet nsmooth nsmooth1 ocparcel octentwt ofilp ofsemb ortho
 patch pchain pchain1 phaserot pick pick2 pick3 plane poly polyfit
 psmig2 regr rdiv rect1 reshape riesz rsin sc seislet1 semblance