b'\n \n \n
 
\n  
sfautofocusing (4.0)
index
user/fbroggin/Mautofocusing.c
\n Marchenko-Wapenaar-Broggini iterative scheme\n

\n \n \n \n \n \n
 
\n Synopsis
       sfautofocusing < Fplus.rsf refl=FRefl.rsf > FGp.rsf Gm=FGm.rsf G=FG.rsf H=FH.rsf p=Fp.rsf q=Fq.rsf window=Ftwin.rsf verb=n conj=n twin=n pandq=n Gtot=n Htot=n niter=1 nshots=1 scale=1.0 eps=1e-4 shift=5

\nsfmarchenko < downgoing.rsf refl=REFL_000.rsf conj=y verb=n Gtot=y niter=21 nshots=401 scale=1 eps=1e-4 shift=5 Gm=Gm.rsf G=G.rsf> Gp.rsf
\n
\n======= INPUTS ============
\n
\np0plus.rsf: initial downgoing wavefield
\n
\nREFL_000.rsf: Fourier transform of the reflection response
\n
\n======= PARAMETERS ========
\n
\nconj = [y]/n\t- complex-conjugation of the first input (corresponds to time-reversal in time)
\nverb = y/[n]\t- verbosity flag
\ntwin = y/[n]\t- returns the timewindow as one of the outputs (window=window.rsf)
\npandq = y/[n]\t- pandq=true: returns p and q, pandq=false returns Gp and Gm
\nGtot = y/[n]\t- Gtot=true returns G=Gp+Gm
\nHtot = y/[n]\t- Htot=true returns H=Gp-Gm
\nniter = 1\t\t- number of iterations
\nnshots = 1\t\t- number of shots in the reflection response
\nscale = 1.0\t- scale factor (often due to resampling)
\neps = 1e-4\t\t- threshold for the timewindow
\nshift = 5\t\t- shift in samples for the timewindow
\n\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
file G=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Gm=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
bool Gtot=n [y/n]
\tGtot=true: returns G=Gp+Gm
\n
\n \n\n \n \n
\n  
file H=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
bool Htot=n [y/n]
\tHtot=true: returns H=Gp-Gm
\n
\n \n\n \n \n
\n  
bool conj=n [y/n]
\tcomplex conjugation (time-reversal) flag
\n
\n \n\n \n \n
\n  
float eps=1e-4
\tthreshold for the timewindow
\n
\n \n\n \n \n
\n  
int niter=1
\tnumber of iterations
\n
\n \n\n \n \n
\n  
int nshots=1
\tnumber of shots
\n
\n \n\n \n \n
\n  
file p=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
bool pandq=n [y/n]
\tpandq=true: returns p and q
\n
\n \n\n \n \n
\n  
file q=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
string refl=
\t000.rsf are 7 characters (auxiliary input file name)
\n
\n \n\n \n \n
\n  
float scale=1.0
\tscale factor
\n
\n \n\n \n \n
\n  
int shift=5
\tshift in samples for the timewindow
\n
\n \n\n \n \n
\n  
bool twin=n [y/n]
\treturns the timewindow as one of the outputs
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]
\tverbosity flag
\n
\n \n\n \n \n
\n  
file window=
\tauxiliary output file name
\n
\n \n
'