macos

Author: sfomel

api/c/lapack.h

@@ -5,9 +5,9 @@
 
 /* #define __VECLIB__ */
 
-#include <Accelerate/Accelerate.h>
+/* #include <Accelerate/Accelerate.h> */
 
-/* #include <vecLib/clapack.h> */
+#include <vecLib/clapack.h> 
 
 #define MKL_Complex8 __CLPK_complex 
 #define MKL_Complex16 __CLPK_doublecomplex 

user/fomels/Mcausinv.c

@@ -0,0 +1,52 @@
+/* Smooth derivative by regularized causint inversion */
+/*
+  Copyright (C) 2004 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 <rsf.h>
+
+int main(int argc, char* argv[])
+{
+  int n1, niter;
+  float eps;
+  bool verb;
+  float *x, *y, **xmov;
+  sf_file in, out, mov;
+
+  sf_init(argc,argv);
+  in = sf_input("in");
+  out = sf_output("out");
+
+  mov = sf_output("mov");
+
+  if (!sf_histint(in, "n1", &n1)) sf_error("No n1= in input");
+  if (!sf_getint("niter", &niter)) niter=1;
+  if (!sf_getfloat("eps", &eps)) eps=0.0;
+  if (!sf_getbool("verb", &verb)) verb=false;
+
+  x = sf_floatalloc(n1);
+  y = sf_floatalloc(n1);
+  xmov = sf_floatalloc2(n1,niter);
+
+  sf_floatread(y, n1, in);
+  
+  sf_solver_reg (sf_causint_lop, sf_cgstep, sf_igrad1_lop, n1, n1, n1, x, y, niter, eps, "verb", verb, "xmov", xmov, "end");
+
+  sf_floatwrite(x, n1, out);
+  sf_floatwrite(xmov[0], n1*niter, mov);
+
+  exit(0);
+}

user/fomels/SConstruct

@@ -4,25 +4,25 @@ import bldutil
 
 progs = '''
 abalance analytical angle angle2 approx arrival bdix beamform1 bil1
-bilat2 blur boxcascade cbeamform1 ccausint cchebyshevp cfftexpa-dev
-cfftwave1 cdivn cflow chain2dfft chebvc chebyshev chebyshevp clpf
-cltft cltftfft constperm constpermh constpermh1 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 semblancew shape shapeagc shapefill shearer shift1
-sh1 sh2 similarity similarity2 simenv slschain2d smoothderw smoothreg
-smspray stfchain stfchain2 stcontrib stltft stpf stphase stream
-streamiss sttimefreq taupfit tfchain tf2dprec thin timecont timefreq
-tomo tree trace2 tricascade tristack tristack2 twofreq2 upgrad var2
-velcon velinv vipmig0 vofzperm warp1 warpadd warpscan zero zmarch
-ztrace
+bilat2 blur boxcascade cbeamform1 causinv ccausint cchebyshevp
+cfftexpa-dev cfftwave1 cdivn cflow chain2dfft chebvc chebyshev
+chebyshevp clpf cltft cltftfft constperm constpermh constpermh1
+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 semblancew shape shapeagc
+shapefill shearer shift1 sh1 sh2 similarity similarity2 simenv
+slschain2d smoothderw smoothreg smspray stfchain stfchain2 stcontrib
+stltft stpf stphase stream streamiss sttimefreq taupfit tfchain
+tf2dprec thin timecont timefreq tomo tree trace2 tricascade tristack
+tristack2 twofreq2 upgrad var2 velcon velinv vipmig0 vofzperm warp1
+warpadd warpscan zero zmarch ztrace
 '''
 
 ccprogs = 'anisolr2 isolr2 isolr3 lrvc0 permlr1 permlr2 permlr3'

user/fomels/cfft2w.c

@@ -202,7 +202,7 @@ void icfft2(sf_complex *out /* [n1*n2] */,
     }
 }
 
-void cfft2_finalize()
+void cfft2_finalize(void)
 /*< clean up fftw >*/
 {
 /* make sure everything is back to its pristine state */