Merge branch 'master' of https://github.com/ahay/src

api/c/files.c

@@ -78,7 +78,7 @@ Outputs the number of dimensions dim and a dimension array n[dim] >*/
     return dim;
 }
 
-int sf_memsize()
+int sf_memsize(void)
 /*< Returns memory size by:
   1. checking RSFMEMSIZE environmental variable
   2. using hard-coded "def" constant

book/tccs/micro/SConstruct

@@ -1,3 +1,3 @@
 from rsf.tex import *
 
-End(color='ALL',hires='data sov location0-clean location0-noisy location0-hyb stage0-1 stage0-2 stage0-3 stage-1 stage-2 stage-3',use='amsmath,listings,hyperref')
+End(color='ALL',hires='data sov location0-clean location0-noisy location0-hyb stage0-1 stage0-2 stage0-3 stage-1 stage-2 stage-3',use='hyperref')

book/tccs/micro/paper.tex

@@ -1,7 +1,7 @@
 \published{SEG Expanded Abstracts, 2485-2490, (2015)}
 
 \title{Investigating the possibility of locating microseismic sources using distributed sensor networks}
-\author{Junzhe Sun$^{*1}$, Tieyuan Zhu$^1$, Sergey Fomel$^1$ and Wen-Zhan Song$^{2}$; $^1$The University of Texas at Austin; $^2$Georgia State University}
+\author{Junzhe Sun, Tieyuan Zhu, Sergey Fomel, The University of Texas at Austin, and Wen-Zhan Song, Georgia State University}
 
 \maketitle
 
@@ -88,8 +88,8 @@ \section{Numerical Examples}
 %   \label{threelayer}
 %\end{figure*}
 \inputdir{cross}
-\multiplot*{3}{sov,datatrace-clean,datatrace-noisy}{height=0.145\textheight}{(a) Microseismic source locations overlaid on a three-layer velocity model. (b) Noise-Free microseismic data from three-layer model. (c) Noisy microseismic data from three-layer model. The traces displayed in the right panel correspond to $X=1050\;m$.}
-\multiplot*{3}{location0-clean,location0-noisy,location0-hyb}{height=0.145\textheight}{(a) Microseismic source locations imaged by the cross-correlation imaging condition using noise-free data. (b) Microseismic source locations imaged by the cross-correlation imaging condition using noisy data. (c) Microseismic source locations imaged by the hybrid imaging condition using noisy data.}
+\multiplot{3}{sov,datatrace-clean,datatrace-noisy}{height=0.145\textheight}{(a) Microseismic source locations overlaid on a three-layer velocity model. (b) Noise-Free microseismic data from three-layer model. (c) Noisy microseismic data from three-layer model. The traces displayed in the right panel correspond to $X=1050\;m$.}
+\multiplot{3}{location0-clean,location0-noisy,location0-hyb}{height=0.145\textheight}{(a) Microseismic source locations imaged by the cross-correlation imaging condition using noise-free data. (b) Microseismic source locations imaged by the cross-correlation imaging condition using noisy data. (c) Microseismic source locations imaged by the hybrid imaging condition using noisy data.}
 
 \subsection{Three-layer model}
 
@@ -105,7 +105,7 @@ \subsection{Three-layer model}
 \plot{sov}{height=0.16\textheight}{Microseismic source locations overlaid on Marmousi velocity model. A total of three stages of microseismicity are assumed, each with six hypocenters. \label{fig:marmsov}}
 \plot{datatrace}{height=0.16\textheight}{Noisy microseismic data from Marmousi model. The trace in the right panel corresponds to $X=1.2\;km$.\label{fig:datamarm}}
 %\multiplot{2}{location0,location}{width=0.4\textwidth}{Imaged microseismic source locations using hybrid imaging condition of: (a) noisy data from surface array; (b) noisy data from both surface array and downhole sensors. \label{fig:marmloc}}
-\multiplot*{6}{stage0-1,stage0-2,stage0-3,stage-1,stage-2,stage-3}{height=0.15\textheight}{Accumulated microseismicity calculated by the hybrid imaging condition using: (a)-(c) noisy data from surface array; (d)-(f) noisy data from both surface array and downhole sensors. Each column is produced after one stage of perforation. \label{fig:marmloc}}
+\multiplot{6}{stage0-1,stage0-2,stage0-3,stage-1,stage-2,stage-3}{height=0.15\textheight}{Accumulated microseismicity calculated by the hybrid imaging condition using: (a)-(c) noisy data from surface array; (d)-(f) noisy data from both surface array and downhole sensors. Each column is produced after one stage of perforation. \label{fig:marmloc}}
 
 \subsection{Marmousi model}
 

book/tccs/owe/SConstruct

@@ -0,0 +1,3 @@
+from rsf.tex import *
+
+End(use='amsmath')

book/tccs/owe/Sage/ostolt.sage

@@ -0,0 +1,8 @@
+k = var('k')
+def ost(t, w):
+    return plot(w*cos(t) - k*sin(t), k, -2, 2)
+def env(w):
+    return plot(sqrt(w^2 + k^2), k, -2, 2, thickness=5)
+lines = [ost(t*pi/20,1) for t in range(-10, 11)]
+plot = sum(lines)+env(1) 
+plot.save(frame=True, axes=False, axes_labels=['Wavenumber (rad)', 'Frequency (rad)'], filename='junk_sage.pdf')

book/tccs/owe/owe.bib

@@ -0,0 +1,74 @@
+@Article{pnas,
+  author = 	 {S[ergey] Fo\-mel and J[ames] A[] Sethian},
+  title = 	 {Fast phase-space computation of multiple
+  arrivals},
+  journal = 	 {Proc. of the Nat. Acad. of Sciences},
+  year = 	 2002,
+  volume =	 99,
+  pages =	 {7329-7334}
+}
+
+@Article{jcp,
+  author = 	 {S[ergey] Fo\-mel and J[ames] A[] Sethian},
+  title = 	 {Fast Algorithms for Computing 
+Multiple-Arrival Multiple-Source Non-Viscosity
+Solutions of Static {Hamilton-Jacobi} Equations},
+  journal = 	 {Journal of Computational Physics, submitted},
+  year = 	 2003
+}
+
+@Book{maslov,
+  author =	 {V[] P[] Maslov and M[] V[] Fedoriuk},
+  title = 	 {Semi-classical approximations in
+quantum mechanics},
+  publisher = 	 {Reidel},
+  year = 	 1981,
+  address =	 {Dordrecht}
+}
+
+@Book{norm,
+  author =	 {N[] Bleisten and J[] K[] Cohen and J[] W[] Stockwell},
+  title = 	 {Mathematics of multidimensional seismic imaging, migration, and inversion},
+  publisher = 	 {Springer},
+  year = 	 2001
+}
+
+@Book{iei,
+  author =	 {J[on] F[] Claerbout},
+  title = 	 {Imaging the {Earth's} Interior},
+  publisher = 	 {Blackwell},
+  year = 	 1985
+}
+
+@Book{santo,
+  author =	 {J[] A[] DeSanto},
+  title = 	 {Scalar wave theory: {Green's} functions and applications},
+  publisher = 	 {Springer-Verlag},
+  year = 	 1992
+}
+
+@Article{curvelet,
+  author = 	 {E[mmanuel] J[] Cand\'{e}s and L[aurent] Demanet},
+  title = 	 {The curvelet representation of wave propagators is optimally sparse},
+  journal = 	 {Communications on Pure and Applied Mathematics},
+  year = 	 2005,
+  volume =	 58,
+  pages =	 {1472-1528}
+}
+
+@inproceedings{seislet,
+   author = {S[] Fomel},
+   booktitle = {76rd Ann. Internat. Mtg.},
+   pages = {2847-2850},
+   publisher = {Soc. of Expl. Geophys.},
+   title = {Towards the seislet transform},
+   year = {2006}
+}
+
+@Book{sethian,
+  author =	 {J[] A[] Sethian},
+  title = 	 {Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science},
+  publisher = 	 {Cambridge Univ. Press},
+  year = 	 1999
+}
+