https://github.com/cran/RandomFields
Revision 6fb75ab304ae785fdefa87c0017d8d39dbb762d2 authored by Martin Schlather on 19 December 2005, 00:00:00 UTC, committed by Gabor Csardi on 19 December 2005, 00:00:00 UTC
1 parent 4d4c43e
Tip revision: 6fb75ab304ae785fdefa87c0017d8d39dbb762d2 authored by Martin Schlather on 19 December 2005, 00:00:00 UTC
version 1.3.13
version 1.3.13
Tip revision: 6fb75ab
genu
extremes.cc: if (kc->simugrid) {
extremes.cc: if (kc->genuine_dim[d]) {
extremes.cc: if ((next = ((min[d] > kc->x[XENDD[d]]) ||
extremes.cc: (max[d] < kc->x[XSTARTD[d]])))) { break;}
extremes.cc: if (min[d] < kc->x[XSTARTD[d]]) {start[d]=0;}
extremes.cc: else start[d] = (int) ((min[d] - kc->x[XSTARTD[d]]) /
extremes.cc: kc->x[XSTEPD[d]]);
extremes.cc: end[d] = (int) ((max[d] - kc->x[XSTARTD[d]]) /
extremes.cc: kc->x[XSTEPD[d]]);
extremes.cc: coord[kc->idx[d]] = startcoord[d] =
extremes.cc: kc->x[XSTARTD[d]] +(double)start[d] * kc->x[XSTEPD[d]];
extremes.cc: coord[kc->idx[d]] += kc->x[XSTEPD[d]];
extremes.cc: coord[kc->idx[d]] = startcoord[d];
extremes.cc: coord[kc->idx[d]] += kc->x[XSTEPD[d]];
extremes.cc: for (d=0; d<dim; d++) y[d] = kc->x[j++];
Hyperfcts.cc: // printf("get: %s %d %f\n", CovList[kc->nr].name, instance, kc->param[KAPPA]);
Hyperfcts.cc: if (kc->nr == EXPONENTIAL) {
Hyperfcts.cc: } else if (kc->nr == CAUCHY) {
Hyperfcts.cc: if (kc->nr == GENERALISEDCAUCHY || kc->nr==STABLE) {
Hyperfcts.cc: } else if (kc->nr == WHITTLEMATERN) {
Hyperfcts.cc: if (kc->param[KAPPA] <= thres[0]) {
Hyperfcts.cc: } else if (kc->param[KAPPA] <= thres[1]) {
Hyperfcts.cc: param = kc->param;
Hyperfcts.cc: nsub= (int) (kc->param[HYPERNR]);
Hyperfcts.cc: store[v] = kc->aniso[0];
Hyperfcts.cc: kc->aniso[0] = d;
Hyperfcts.cc: if ((kc->nr == GENERALISEDCAUCHY && kc->param[KAPPA] <= 1.0) ||
Hyperfcts.cc: (kc->nr == STABLE && kc->param[KAPPA] <= 1.0) ||
Hyperfcts.cc: (kc->nr == EXPONENTIAL) ||
Hyperfcts.cc: (kc->nr == WHITTLEMATERN && kc->param[KAPPA] <= 0.5)) {
Hyperfcts.cc: } else if (kc->nr == BROWNIAN) {
Hyperfcts.cc: q[INTRINSIC_RAWR] = (kc->truetimespacedim <= 2
Hyperfcts.cc: ? ((kc->param[KAPPA] <= 1.5) ? 1.0 : 2.0)
Hyperfcts.cc: : ((kc->param[KAPPA] <= 1.0) ? 1.0 : 2.0));
Hyperfcts.cc: kc->param[INTRINSIC_RAWR] *= 2.0;
Hyperfcts.cc: param = kc->param;
Hyperfcts.cc: nsub= (int) (kc->param[HYPERNR]);
Hyperfcts.cc: store[v] = kc->aniso[0];
Hyperfcts.cc: kc->aniso[0] = d;
Hyperfcts.cc: if (v<nsub && CovList[kc->nr].variogram && kc->op)
Hyperfcts.cc: param = kc->param;
Hyperfcts.cc: if (kc->genuine_last_dimension)
Hyperfcts.cc: nsub= (int) (kc->param[HYPERNR]);
Hyperfcts.cc: endfor = ANISO + kc->dim * kc->dim - 1;
Hyperfcts.cc: if (kc->param[i]!=0.0) {
Hyperfcts.cc: kc->param[i], 1 - (ANISO - i), v); // i - ANISO + 1 gibt
MPP.cc: if ((kc->method==AdditiveMpp) && kc->left
MPP.cc: && (kc->param[VARIANCE]>0)) {
MPP.cc: assert(kc->nr>=0 && kc->nr<currentNrCov);
MPP.cc: cov = &(CovList[kc->nr]);
MPP.cc: if (key->Time && !kc->genuine_last_dimension &&
MPP.cc: if ((error=cov->check(kc->param, kc->truetimespacedim,
MPP.cc: if ((v<key->ncov-1 && kc->op)) {
MPP.cc: kc->left=false;
MPP.cc: timespacedim = kc->truetimespacedim;
MPP.cc: if (kc->simugrid) {
MPP.cc: s->min[kc->idx[d]] = kc->x[XSTARTD[d]];
MPP.cc: s->length[kc->idx[d]]= kc->x[XSTEPD[d]] * (double) (key->length[d] - 1);
MPP.cc: if (kc->x[ix+d] < s->min[d]) s->min[d] = kc->x[ix+d];
MPP.cc: if (kc->x[ix+d] > max[d]) max[d] = kc->x[ix+d];
MPP.cc: CovList[kc->nr].add_mpp_scl(s, timespacedim, kc->param);
MPP.cc: s->MppFct = CovList[kc->nr].add_mpp_rnd;
MPP.cc: factor = sqrt(kc->param[VARIANCE] /
MPP.cc: lambda = kc->param[VARIANCE] * s->effectivearea;
MPP.cc: timespacedim = kc->simugrid ? key->timespacedim : kc->truetimespacedim;
MPP.cc: if (kc->simugrid) {
MPP.cc: if (kc->genuinedim[d]) {
MPP.cc: if (min[d]< kc->x[XSTARTD[d]]) {start[d]=0;}
MPP.cc: else start[d] = (int) ((min[d] - kc->x[XSTARTD[d]]) /
MPP.cc: kc->x[XSTEPD[d]]);
MPP.cc: end[d] = 1 + (int) ((max[d] - kc->x[XSTARTD[d]]) /
MPP.cc: kc->x[XSTEPD[d]]);
MPP.cc: coord[kc->idx[d]] = startcoord[d] =
MPP.cc: kc->x[XSTARTD[d]] + (double) start[d] * kc->x[XSTEPD[d]];
MPP.cc: coord[kc->idx[d]] += kc->x[XSTEPD[d]];
MPP.cc: coord[kc->idx[d]] = startcoord[d];
MPP.cc: coord[kc->idx[d]] += kc->x[XSTEPD[d]];
MPP.cc: for (d=0; d<timespacedim; d++, j++) y[d] = kc->x[j];
RFcircembed.cc: dim = kc->truetimespacedim;
RFcircembed.cc: GetCornersOfGrid(key, dim, kc->aniso, sx);
RFcircembed.cc: dummy = key->x[d][XSTEP] * (double) (key->length[d] - 1) * kc->aniso[0];
RFcircembed.cc:// printf("** %d\n", kc->nr);
RFcircembed.cc: memcpy(&(sc->param[ANISO]), &(kc->param[ANISO]), sizeof(double) * n_aniso);
RFcircembed.cc: err = GetOrthogonalUnitExtensions(kc->aniso, key->timespacedim, grid_ext);
RFgetNset.cc: SET_VECTOR_ELT(submodel[i],subi++, mkString(METHODNAMES[kc->method]));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, ScalarInteger(kc->dim));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, ScalarInteger(kc->truetimespacedim));
RFgetNset.cc: Int(kc->length, kc->simugrid ? kc->dim : 0, MAX_INT));
RFgetNset.cc: Int(kc->idx, kc->simugrid ? kc->dim :0, MAX_INT));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, mkString(CovList[kc->nr].name));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, mkString(OP_SIGN[kc->op]));
RFgetNset.cc: Logi(kc->genuine_dim, kc->simugrid ? kc->dim : 0, MAX_INT));
RFgetNset.cc: ScalarLogical(kc->genuine_last_dimension));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, ScalarLogical(kc->simugrid));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, ScalarLogical(!kc->left));
RFgetNset.cc: : Num(kc->param, totalparam, MAX_INT));
RFgetNset.cc: SET_VECTOR_ELT(submodel[i], subi++, (kc->dim<=0) ? allocVector(VECSXP, 0)
RFgetNset.cc: : Mat(kc->aniso, kc->dim, // kc->truetimespacedim
RFgetNset.cc: kc->dim, MAX_INT)
RFgetNset.cc: if (kc->x == NULL)
RFgetNset.cc: Mat(kc->x, kc->simugrid ? 3 : totalpoints,
RFgetNset.cc: kc->truetimespacedim, MAX_INT));
RFhyperplan.cc: optdim=2; // falls dies gelockert wird, so kc->idx[d] nicht vergessen!
RFhyperplan.cc: if (kc->method==Hyperplane && kc->left) {
RFhyperplan.cc: cov = &(CovList[kc->nr]);
RFhyperplan.cc: assert(kc->nr>=0 && kc->nr<currentNrCov);
RFhyperplan.cc: assert(kc->param[VARIANCE] >= 0.0);
RFhyperplan.cc: if ((key->ncov>v+1 && kc->op) ||
RFhyperplan.cc: timespacedim = kc->truetimespacedim;
RFhyperplan.cc: v += (int) kc->param[HYPERNR];
RFhyperplan.cc: if ((error = cov->check(kc->param, timespacedim, Hyperplane))
RFhyperplan.cc: kc->left = false;
RFhyperplan.cc: timespacedim = kc->truetimespacedim;
RFhyperplan.cc: GetCenterAndDiameter(key, kc->simugrid, timespacedim, kc->truetimespacedim,
RFhyperplan.cc: kc->x, kc->aniso, s->center, s->rx, &(s->radius));
RFhyperplan.cc: for (d=0; d<kc->truetimespacedim; d++) s->rx[d] *= 0.5;
RFhyperplan.cc: timespacedim = kc->truetimespacedim;
RFhyperplan.cc: deltax = kc->x[XSTEPD[0]];
RFhyperplan.cc: deltay = kc->x[XSTEPD[1]];
RFhyperplan.cc: if (kc->simugrid) {
RFhyperplan.cc: for (gy=kc->x[XSTARTD[1]], resindex=j=0; j<key->length[1]; j++) {
RFhyperplan.cc: for (gx= kc->x[XSTARTD[0]], i=0; i<key->length[0]; i++,
RFhyperplan.cc: if ((cell=determine_cell(kc->x[j], kc->x[j+1], hx, hy, hr,
RFnugget.cc:// printf("nugget %d %d %d %s \n", v, kc->nr, key->ncov, CovList[kc->nr].name);
RFnugget.cc: assert((kc->nr >= 0) && (kc->nr < currentNrCov));
RFnugget.cc: assert(kc->param[VARIANCE] >= 0.0);
RFnugget.cc: cov = &(CovList[kc->nr]);
RFnugget.cc: if ((kc->method==Nugget) && (kc->left)) {
RFnugget.cc: quotient = kc->param[nonzero_pos]/first->param[nonzero_pos];
RFnugget.cc: equal &= fabs(first->param[w] * quotient - kc->param[w]) <
RFnugget.cc: (fabs(kc->param[w]) + 1.0*(kc->param[w]==0.0))
RFnugget.cc: kc->left = false;
RFnugget.cc: s->simple = key->timespacedim == kc->truetimespacedim;
RFnugget.cc: s->simugrid = kc->simugrid;
RFnugget.cc: if (kc->simugrid) {
RFnugget.cc: s->diag[d] = kc->param[ANISO + i];
RFnugget.cc: ordering(kc->x, key->totalpoints, kc->truetimespacedim, pos);
RFnugget.cc: if (equal(oldpos, pos[i], kc->x, kc->truetimespacedim))
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: // printf("%d %d %d %d\n", v, kc->method, Method, (int) kc->left);
RFsimu.cc: if (kc->method==Method && kc->left) {
RFsimu.cc: assert((kc->nr >= 0) && (kc->nr < currentNrCov));
RFsimu.cc: assert(kc->param[VARIANCE] >= 0.0);
RFsimu.cc: error = cov->check(kc->param, kc->truetimespacedim, Method);
RFsimu.cc: (int) kc->param[HYPERNR] : 0;
RFsimu.cc: kc->left = false;
RFsimu.cc: kc->left = false;
RFsimu.cc: if (kc->simugrid) {
RFsimu.cc: kc->length[d] = kc->genuine_dim[d] ? key->length[d] : 1;
RFsimu.cc:// kc->length[d], kc->genuine_dim[d], key->length[d]);
RFsimu.cc: return Transform2NoGrid(key, kc->aniso, kc->truetimespacedim,
RFsimu.cc: kc->simugrid, &(kc->x));
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: z = fabs(x[0] * kc->aniso[0]);
RFsimu.cc: fct[w1] = kc->param[VARIANCE] * cov->cov(&z, kc->param, dim);
RFsimu.cc: abl[w1] = kc->param[VARIANCE] * cov->derivative(&z, kc->param)
RFsimu.cc: * kc->aniso[0];
RFsimu.cc: snd[w1] = kc->param[VARIANCE] * cov->secondderivt(&z, kc->param)
RFsimu.cc: * kc->aniso[0] * kc->aniso[0];
RFsimu.cc:// cov->secondderivt(&z, kc->param), cov->name, cov->secondderivt);
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: z = fabs(x[0] * kc->aniso[0]);
RFsimu.cc: fct[w] = kc->param[VARIANCE] * cov->cov(&z, kc->param, dim);
RFsimu.cc: abl[w] = kc->param[VARIANCE] * cov->derivative(&z, kc->param)
RFsimu.cc: * kc->aniso[0];
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: var *= kc->param[VARIANCE];
RFsimu.cc: for (k=0; k<kc->truetimespacedim; k++) z[k]=0.0;
RFsimu.cc: for (ani=0, k=0; k<kc->truetimespacedim; k++) {
RFsimu.cc: z[k] += kc->aniso[ani++] * x[j];
RFsimu.cc: zw *= cov->cov(z, kc->param, kc->truetimespacedim);
RFsimu.cc: endfor = kc->truetimespacedim - 1;
RFsimu.cc: nsub = (int) kc->param[HYPERNR];
RFsimu.cc: zw *= cov->cov(z, kc->param, subdim);
RFsimu.cc:// z[0], z[1], z[2], kc->param[VARIANCE],
RFsimu.cc:// kc->param[KAPPA1], kc->param[KAPPA2], kc->param[KAPPA3],
RFsimu.cc:// subdim, cov->name, cov->cov(z, kc->param, subdim), zw);
RFsimu.cc: zw *= cov->cov(z, kc->param, 1 + (int)(cov_type==SPACEISOTROPIC));
RFsimu.cc: if ((v<ncovM1) && (kc->op==0)) break;
RFsimu.cc: //printf("covfct %d %d\n", v, kc->nr);
RFsimu.cc: cov = &(CovList[kc->nr]); /* cov needed in FORMULA for Vario */
RFsimu.cc: nsub = (int) kc->param[HYPERNR];
RFsimu.cc: z[0] = d * kc->aniso[0];
RFsimu.cc:// d, kc->aniso[0]);
RFsimu.cc: zw *= cov->cov(z, kc->param, subdim);
RFsimu.cc: zz = d * kc->aniso[0];
RFsimu.cc: var *= kc->param[VARIANCE];
RFsimu.cc: zw *= cov->cov(&zz, kc->param, 1);
RFsimu.cc: if ((v<ncovM1) && (kc->op==0)) break;
RFsimu.cc:// kc->param[KAPPA], v, covlist[v], kc->nr);
RFsimu.cc:// *x, z[0], z[1], result, d, kc->aniso[0],
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: *equal &= !memcmp(kc->param, ParamList, sizeof(double) * kappasP1);
RFsimu.cc: memcpy(kc->param, ParamList, sizeof(double)* kappasP1);
RFsimu.cc: if (kc->param[VARIANCE]<0.0) return ERRORNEGATIVEVAR;
RFsimu.cc: *equal &= kc->op == op[v];
RFsimu.cc: kc->op = op[v];
RFsimu.cc:// kc->op, false xor false, false xor true, true xor true);
RFsimu.cc: xor (kc->op == 2))
RFsimu.cc: GetNaturalScaling(&(kc->nr), &(kc->param[KAPPA]),
RFsimu.cc: *equal &= !memcmp(&(kc->param[ANISO]), &(param[ANISO]),
RFsimu.cc: memcpy(&(kc->param[ANISO]), &(param[ANISO]), sizeof(double) * pAniso);
RFsimu.cc: GetTrueDim(anisotropy, timespacedim, kc->param, cov->type,
RFsimu.cc: &(kc->genuine_last_dimension),
RFsimu.cc: &(kc->truetimespacedim), kc->aniso);
RFsimu.cc: kc->dim = timespacedim;
RFsimu.cc: if (cov->type==SPACEISOTROPIC && !kc->genuine_last_dimension)
RFsimu.cc: if (kc->truetimespacedim==0) return ERRORTRIVIAL;
RFsimu.cc: error = check_within_range(kc->param, cov, kc->truetimespacedim, NULL);
RFsimu.cc: cov->info(kc->param, &maxdim, &dummy);
RFsimu.cc:// printf("%d %d\n", maxdim, kc->truetimespacedim);
RFsimu.cc: if (maxdim < kc->truetimespacedim) return ERRORCOVNOTALLOWED;
RFsimu.cc: kc->simugrid = // kc->param[ANISO] !!
RFsimu.cc: grid && (!anisotropy || is_diag(&(kc->param[ANISO]), timespacedim));
RFsimu.cc: if (kc->simugrid) {
RFsimu.cc: if ((kc->genuine_dim[d] = param[ANISO + n] != 0.0)) {
RFsimu.cc: kc->idx[d] = ii++;
RFsimu.cc: kc->idx[d] = kc->truetimespacedim + d - ii;
RFsimu.cc: kc->genuine_dim[d] = true;
RFsimu.cc: kc->idx[d] = d;
RFsimu.cc: kc->simugrid = false; // not programmed yet -- save version !!!!!!
RFsimu.cc: error = cov->check(kc->param, kc->truetimespacedim, Nothing);
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc:// memcpy(&(kc->param[ANISO]), &(keycov[v+1].param[ANISO]),
RFsimu.cc: memcpy(&(kc->param[ANISO]), &(keycov[v+1].param[ANISO]),
RFsimu.cc: GetTrueDim(anisotropy, timespacedim, kc->param,
RFsimu.cc: &(kc->genuine_last_dimension), //
RFsimu.cc: &(kc->truetimespacedim), kc->aniso);//
RFsimu.cc: kc->aniso[0] = RF_NAN;
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: v += (int) kc->param[HYPERNR];
RFsimu.cc: if (CovList[kc->nr].cov == nugget && meth!=Nugget && ncov>1 &&
RFsimu.cc: kc->truetimespacedim < timespacedim)) {
RFsimu.cc: key->active &= kc->method == meth;
RFsimu.cc: cov = &(CovList[kc->nr]);
RFsimu.cc: error=cov->check(kc->param, 2, TBM2);
RFsimu.cc: kc->param[TBM2NUM] =
RFsimu.cc: if (GENERAL_PRINTLEVEL > 1 && (kc->param[TBM2NUM]!=0.0))
RFsimu.cc: if (kc->method >= (int) Nothing) {
RFsimu.cc:// printf("kc->method=%d %d %d \n", v, kc->method, (int) Nothing);
RFsimu.cc: if (CovList[kc->nr].type == ISOHYPERMODEL ||
RFsimu.cc: CovList[kc->nr].type == ANISOHYPERMODEL) {
RFsimu.cc: int i, nc = (int) kc->param[HYPERNR];
RFsimu.cc: if (kc->method != key->cov[i].method) {
RFsimu.cc: if (v>0 && op[v-1] && kc->method!=key->cov[v-1].method) {
RFsimu.cc: kc->method==CircEmbedIntrinsic) {
RFsimu.cc: key->ML[v][0] = kc->method;
RFsimu.cc: kc->left = left[v];
RFsimu.cc:// if (kc->left) {
RFsimu.cc:// if (kc->x != NULL) { free(kc->x); kc->x = NULL; }
RFsimu.cc: if (kc->left) {
RFsimu.cc: PRINTF("'%s' not initiated yet\n", CovList[kc->nr].name);
RFsimu.cc: Merr = kc->method;
RFsimu.cc: kc->left = left[w];
RFsimu.cc:// if (kc->left) {
RFsimu.cc: CovList[kc->nr].name);
RFspectral.cc: if ((kc->method==SpectralTBM) && (kc->left)) {
RFspectral.cc: timespacedim = kc->truetimespacedim;
RFspectral.cc: assert((kc->nr >= 0) && (kc->nr < currentNrCov));
RFspectral.cc: assert(kc->param[VARIANCE] >= 0.0);
RFspectral.cc: cov = &(CovList[kc->nr]);
RFspectral.cc: if ((Xerror=check_within_range(kc->param, cov, 2, " (spectral TBM)"))
RFspectral.cc: (Xerror=cov->check(kc->param, timespacedim, SpectralTBM))
RFspectral.cc: if (store_variance == kc->param[VARIANCE]) actcov++;
RFspectral.cc: kc->left=true;
RFspectral.cc: store_variance = kc->param[VARIANCE];
RFspectral.cc: kc->left = false;
RFspectral.cc: assert((kc->truetimespacedim>=1) && (kc->truetimespacedim<=2));
RFspectral.cc: sqrttwodivbyn = sqrt(kc->param[VARIANCE] * 2.0 / (double) SPECTRAL_LINES);
RFspectral.cc: Amp= (s->randomAmplitude[v])(kc->param);
RFspectral.cc: if (kc->simugrid) {
RFspectral.cc:// printf("%d %d \n", d, kc->length[d]);
RFspectral.cc: if (kc->length[d] != 1) {
RFspectral.cc: segt += kc->x[XSTARTD[d]] * cossin[idxcossind];
RFspectral.cc: inc[d] = kc->x[XSTEPD[d]] * cossin[idxcossind++];
RFspectral.cc: if (kc->truetimespacedim==1) {
RFspectral.cc: res[nx] += cos(cp * kc->x[nx] + VV);
RFspectral.cc: } else { // kc->truetimespacedim==2
RFspectral.cc: res[nx] += cos(cp * kc->x[twonx] + sp * kc->x[twonx+1] + VV);
RFtbm.cc: cov = &(CovList[kc->nr]);
RFtbm.cc: z[0] = kc->aniso[0] * x[0];
RFtbm.cc: z[1] = kc->aniso[TBMTIME] * x[1];
RFtbm.cc: if (z[0]==0.0) zw *= cov->cov(&(z[1]), kc->param, dim);
RFtbm.cc: zw *= TBM2integr(kc->nr, z[0], z[1], kc->param, 2, 2);
RFtbm.cc: zw *= kc->param[VARIANCE] * z[0] / r *
RFtbm.cc: (cov->cov_tbm2(&r, kc->param, dim) +
RFtbm.cc: TBM2integr_time(z[0] / r, r, cov->cov, kc->param, dim) );
RFtbm.cc: if (z[0]==0.0) zw *= cov->cov(z, kc->param, 2);
RFtbm.cc: else if (kc->param[TBM2NUM] == 0.0) {
RFtbm.cc: zw *= cov->cov_tbm2(z, kc->param);
RFtbm.cc: zw *= TBM2integr(kc->nr, z[0], z[1], kc->param, 2,
RFtbm.cc: zw *= kc->param[VARIANCE];
RFtbm.cc: cov = &(CovList[kc->nr]);
RFtbm.cc: z = x[0] * kc->aniso[0];
RFtbm.cc: else if (kc->param[TBM2NUM] == 0.0) {
RFtbm.cc: zw = cov->cov_tbm2(&z, kc->param);
RFtbm.cc: zw = TBM2integr(kc->nr, z, 0, kc->param, 2, cov->cov_tbm3 == NULL);
RFtbm.cc: result += kc->param[VARIANCE] * zw;
RFtbm.cc: cov = &(CovList[kc->nr]);
RFtbm.cc: z[0] = kc->aniso[0] * x[0];
RFtbm.cc: z[1] = kc->aniso[TBMTIME] * x[1];
RFtbm.cc: else abl[w] = kc->param[VARIANCE] * kc->aniso[0] * fabs(z[0]) /
RFtbm.cc: r * cov->derivative(&r, kc->param);
RFtbm.cc: zw *= (fct[w] = kc->param[VARIANCE] * cov->cov(&r, kc->param, dim));
RFtbm.cc: abl[w] = kc->param[VARIANCE] *
RFtbm.cc: kc->aniso[0] * cov->derivative(z, kc->param);
RFtbm.cc: zw *= (fct[w] = kc->param[VARIANCE] * cov->cov(z, kc->param, dim));
RFtbm.cc: cov = &(CovList[kc->nr]);
RFtbm.cc: z[0] = x[0] * kc->aniso[0];
RFtbm.cc: return kc->param[VARIANCE] * cov->cov_tbm3(z, kc->param);
RFtbm.cc: if ((kc->method==method) && (kc->left)) {
RFtbm.cc: ce_dim2 |= CovList[kc->nr].type==SPACEISOTROPIC ||
RFtbm.cc: kc->truetimespacedim == 1 + tbm_dim;
RFtbm.cc: while (!ce_dim2 && (++v)<key->ncov && kc->op) {
RFtbm.cc: ce_dim2 |= CovList[kc->nr].type==SPACEISOTROPIC;
RFtbm.cc: quot = kc->param[nonzero_pos] / first->param[nonzero_pos];
RFtbm.cc: if (fabs(first->param[w] * quot - kc->param[w]) >=
RFtbm.cc: (fabs(kc->param[w]) + (double)(kc->param[w]==0.0)) *
RFtbm.cc: if (kc->truetimespacedim > tbm_dim) {
RFtbm.cc: if ((kc->method==method) && (kc->left)) {
RFtbm.cc: cov = &(CovList[kc->nr]);
RFtbm.cc: assert((kc->nr>=0) && (kc->nr<currentNrCov));
RFtbm.cc: assert(kc->param[VARIANCE]>=0.0);
RFtbm.cc: if (v<key->ncov-1 && kc->op) {
RFtbm.cc: kc->left = true;
RFtbm.cc: if ((Xerror=check_within_range(kc->param, cov, tbm_dim + (int) ce_dim2,
RFtbm.cc: Xerror=cov->check(kc->param, tbm_dim + (int) ce_dim2, method);
RFtbm.cc: quot = kc->param[nonzero_pos] / first->param[nonzero_pos];
RFtbm.cc: equal = fabs(first->param[w] * quot - kc->param[w])
RFtbm.cc: < (fabs(kc->param[w]) + (double)(kc->param[w]==0.0)) *
RFtbm.cc: kc->left = true;
RFtbm.cc: memcpy(&(ParamList[cum_nParam[actcov]]), kc->param,
RFtbm.cc: covnr[actcov] = kc->nr;
RFtbm.cc: op[actcov] = kc->op;
RFtbm.cc: if (kc->param[w]!=0.0) {
RFtbm.cc: ParamList[Aniso + 3] = kc->param[lasttimecomponent];
RFtbm.cc: if (kc->param[w] != 0.0) {
RFtbm.cc: if (spacecomponent_found = (kc->param[w]!=0.0)) break;
RFtbm.cc: kc->left=false;
RFtbm.cc: } // kc->left
RFtbm.cc: newadditive = kc->op != 0;
Computing file changes ...