void PCR::setup(){ set_abd_current(); setup_X(); setup_beta(); setup_aliases(); set_observers(); }
static void init_predictions(PRED_AT w) { int i; DPOINT *bp; DATA **d = NULL; #ifdef WITH_SPIRAL DATA_GRIDMAP *grid; #endif est = (double *) emalloc(get_n_outfile() * sizeof(double)); bp = get_block_p(); d = get_gstat_data(); switch (w) { case AT_POINTS: if (o_filename == NULL) ErrMsg(ER_VARNOTSET, "please specify output file"); write_points(o_filename, val_data, NULL, NULL, get_n_outfile()); if (bp->x == -1.0) { /* set default */ bp->x = bp->y = 1.0; pr_warning("default block size set to: dx=1, dy=1"); } break; case AT_GRIDMAP: /* open mask files: */ get_maskX(NULL, NULL, 0, 0); /* re-initializes static arrays */ masks = (GRIDMAP **) emalloc(get_n_masks() * sizeof(GRIDMAP *)); for (i = 0; i < get_n_masks(); i++) masks[i] = check_open(get_mask_name(i), i); /* read as float */ if (n_pred_locs > 0) strata_min = floor(masks[0]->cellmin); outmap = (GRIDMAP **) emalloc(get_n_outfile() * sizeof(GRIDMAP *)); printlog("initializing maps "); for (i = 0; i < get_n_outfile(); i++) { if (get_outfile_namei(i) != NULL) { printlog("."); /* creating maps ..... */ if (get_method() == ISI) masks[0]->celltype = CT_UINT8; outmap[i] = map_dup(get_outfile_namei(i), masks[0]); } else outmap[i] = NULL; } printlog("\n"); if (bp->x == -1.0) { /* set default to map cellsize */ bp->x = masks[0]->cellsizex; bp->y = masks[0]->cellsizey; pr_warning("default block size set to dx=%g, dy=%g", bp->x, bp->y); } for (i = 0; i < get_n_vars(); i++) { if (d[i]->dummy) { d[i]->minX = masks[0]->x_ul + 0.5 * masks[0]->cellsizex; d[i]->maxX = masks[0]->x_ul + masks[0]->cellsizex * (masks[0]->cols - 0.5); d[i]->maxY = masks[0]->y_ul - 0.5 * masks[0]->cellsizey; d[i]->minY = masks[0]->y_ul - masks[0]->cellsizey * (masks[0]->rows - 0.5); d[i]->minZ = d[i]->maxZ = 0.0; } if (d[i]->togrid) datagrid_rebuild(d[i], 1); } break; } if (gl_nsim > 1) init_simulations(d); if (is_simulation(get_method()) && get_n_beta_set() != get_n_vars()) setup_beta(d, get_n_vars(), gl_nsim); } /* init_predictions() */