Пример #1
0
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);
}
Пример #2
0
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);
    }
}