static WebData AllocUserData(void) { int i, ngrp = NGRP, ns = NS; WebData wdata; wdata = (WebData) malloc(sizeof *wdata); for(i=0; i < ngrp; i++) { (wdata->P)[i] = denalloc(ns); (wdata->pivot)[i] = denallocpiv(ns); } wdata->rewt = N_VNew_Serial(NEQ+1); return(wdata); }
static UserData AllocUserData(void) { int jx, jy; UserData data; data = (UserData) malloc(sizeof *data); for (jx=0; jx < MX; jx++) { for (jy=0; jy < MY; jy++) { (data->P)[jx][jy] = denalloc(NUM_SPECIES); (data->Jbd)[jx][jy] = denalloc(NUM_SPECIES); (data->pivot)[jx][jy] = denallocpiv(NUM_SPECIES); } } return(data); }
static UserData AllocUserData(void) { int jx, jz; UserData data; data = (UserData) malloc(sizeof *data); for (jx=0; jx < MX; jx++) { for (jz=0; jz < MZ; jz++) { (data->P)[jx][jz] = denalloc(NUM_SPECIES); (data->Jbd)[jx][jz] = denalloc(NUM_SPECIES); (data->pivot)[jx][jz] = denallocpiv(NUM_SPECIES); } } data->p = (realtype *) malloc(NP*sizeof(realtype)); return(data); }
static PreconData AllocPreconData(UserData fdata) { int lx, ly; PreconData pdata; pdata = (PreconData) malloc(sizeof *pdata); pdata->f_data = fdata; for (lx = 0; lx < MXSUB; lx++) { for (ly = 0; ly < MYSUB; ly++) { (pdata->P)[lx][ly] = denalloc(NVARS, NVARS); (pdata->Jbd)[lx][ly] = denalloc(NVARS, NVARS); (pdata->pivot)[lx][ly] = denallocpiv(NVARS); } } return(pdata); }
static UserData AllocUserData(void) { int jx, jy; UserData data; data = (UserData) malloc(sizeof *data); for (jx = 0; jx < MXSUB; jx++) { for (jy = 0; jy < MYSUB; jy++) { (data->P)[jx][jy] = denalloc(NUM_SPECIES, NUM_SPECIES); (data->pivot)[jx][jy] = denallocpiv(NUM_SPECIES); } } acoef = denalloc(NUM_SPECIES, NUM_SPECIES); bcoef = (realtype *)malloc(NUM_SPECIES * sizeof(realtype)); cox = (realtype *)malloc(NUM_SPECIES * sizeof(realtype)); coy = (realtype *)malloc(NUM_SPECIES * sizeof(realtype)); return(data); }
static UserData AllocUserData(MPI_Comm comm, long int local_N, long int SystemSize) { long int ix, jy; UserData webdata; webdata = (UserData) malloc(sizeof *webdata); webdata->rates = N_VNew_Parallel(comm, local_N, SystemSize); for (ix = 0; ix < MXSUB; ix++) { for (jy = 0; jy < MYSUB; jy++) { (webdata->PP)[ix][jy] = denalloc(NUM_SPECIES, NUM_SPECIES); (webdata->pivot)[ix][jy] = denallocpiv(NUM_SPECIES); } } webdata->acoef = denalloc(NUM_SPECIES, NUM_SPECIES); webdata->ewt = N_VNew_Parallel(comm, local_N, SystemSize); return(webdata); }