Synopsis
		sfiwigrad < in.rsf > out.rsf source=source.rsf data=data.rsf imask=imask.rsf weight=weight.rsf precon=precon.rsf verb=n nh=0 load=n update=y uts=0 npml=10 prect1=5 prect2=1 prect3=5 pliter=20 plower=0.1 pupper=3. dorder=6 grect1=5 grect2=5 gliter=1 geps=0. gscale=0.5 miter=10 lower=1.5 upper=7.5 order= cost=

Parameters
		string cost= cost functional type (default classic DSO)   string data= auxiliary input file name   int dorder=6 image derivative accuracy order   float geps=0. regularization parameter for Gauss-Newton   int gliter=1 # of Gauss-Newton iterations   int grect1=5 gradient smoothing radius on axis 1   int grect2=5 gradient smoothing radius on axis 2   float gscale=0.5 gradient re-scale   string imask= auxiliary input file name   bool load=n [y/n] load LU   float lower=1.5 lower bound of feasible set   int miter=10 Nonlinear-CG maximum # of iterations   int nh=0 horizontal space-lag   int npml=10 PML width   string order= discretization scheme (default optimal 9-point)   int pliter=20 slope estimation # of linear iterations   float plower=0.1 slope thresholding lower limit   string precon= auxiliary input file name   int prect1=5 slope smoothing radius on axis 1   int prect2=1 slope smoothing radius on axis 2   int prect3=5 slope smoothing radius on axis 3   float pupper=3. slope thresholding upper limit   string source= auxiliary input file name   bool update=y [y/n] Non-linear CG update   float upper=7.5 upper bound of feasible set   int uts=0 number of OMP threads   bool verb=n [y/n] verbosity flag   string weight= auxiliary input file name