next up previous contents [pdf] index

Next: Input parameters Up: sf_solver Previous: Call

Definition

void sf_solver (sf_operator oper   /* linear operator */, 
                sf_solverstep solv /* stepping function */, 
                int nx             /* size of x */, 
                int ny             /* size of dat */, 
                float* x           /* estimated model */, 
                const float* dat   /* data */, 
                int niter          /* number of iterations */, 
               ...                /* variable number of arguments */)
/*< Generic linear solver.
  ---
  Solves
  oper{x}    =~ dat
  ---
  The last parameter in the call to this function should be "end".
  Example: 
  ---
  sf_solver (oper_lop,sf_cgstep,nx,ny,x,y,100,"x0",x0,"end");
  ---
  Parameters in ...:
  ---
  "wt":     float*:         weight      
  "wght":   sf_weight wght: weighting function
  "x0":     float*:         initial model
  "nloper": sf_operator:    nonlinear operator
  "mwt":    float*:         model weight
  "verb":   bool:           verbosity flag
  "known":  bool*:          known model mask
  "nmem":   int:            iteration memory
  "nfreq":  int:            periodic restart
  "xmov":   float**:        model iteration
  "rmov":   float**:        residual iteration
  "err":    float*:         final error
  "res":    float*:         final residual
  >*/ 
{
   ...
}




2011-07-02