void dsrtomo_init(int dim /* model dimension */, int *n /* model size */, float *d /* model sampling */) /*< initialize >*/ { int i; nn = n; nt = 1; for (i=0; i < dim; i++) { nt = nt*n[i]; } upg = upgrad_init(dim,n,d); temp = sf_floatalloc(nt); pstk = sf_floatalloc(nt); }
void fatomo_init(int dim /* model dimension */, int *n /* model size */, float *d /* model sampling */, int nshot /* number of shots */) /*< initialize >*/ { int i, is; nt = 1; for (i=0; i < dim; i++) { nt = nt*n[i]; } ns = nshot; tempt = sf_floatalloc(nt); tempx = sf_floatalloc(nt); upglist = (upgrad *)malloc(ns*sizeof(upgrad)); for (is=0; is < ns; is++) { upglist[is] = upgrad_init(dim,n,d); } }