ptatin3d_energyfv.c:71:  energy->time = time;
ptatin3d_energyfv.c:72:  energy->dt = dt;
ptatin3d_energyfv.c:74:    energy->nsubdivision[0] = nsub[0];
ptatin3d_energyfv.c:75:    energy->nsubdivision[1] = nsub[1];
ptatin3d_energyfv.c:76:    energy->nsubdivision[2] = nsub[2];
ptatin3d_energyfv.c:78:    energy->npoints_macro = (energy->nsubdivision[0] + 1)*(energy->nsubdivision[1] + 1)*(energy->nsubdivision[2] + 1);
ptatin3d_energyfv.c:91:    if (n != 3) SETERRQ(energy->fv->comm,PETSC_ERR_USER,"Must provide 3 values for option -ptatin_energyfv_nsub");
ptatin3d_energyfv.c:92:    ierr = PhysCompEnergyFVSetParams(energy,energy->time,energy->dt,sub);CHKERRQ(ierr);
ptatin3d_energyfv.c:116:  energy->mi_parent[0] = q2_mi[0];
ptatin3d_energyfv.c:117:  energy->mi_parent[1] = q2_mi[1];
ptatin3d_energyfv.c:118:  energy->mi_parent[2] = q2_mi[2];
ptatin3d_energyfv.c:121:  fv_mi[0] = q2_mi[0] * energy->nsubdivision[0];
ptatin3d_energyfv.c:122:  fv_mi[1] = q2_mi[1] * energy->nsubdivision[1];
ptatin3d_energyfv.c:123:  fv_mi[2] = q2_mi[2] * energy->nsubdivision[2];
ptatin3d_energyfv.c:125:  ierr = fvgeometry_dmda3d_create_from_element_partition(energy->fv->comm,decomp,fv_mi,&energy->dmv);CHKERRQ(ierr);
ptatin3d_energyfv.c:127:  ierr = DMDAGetElementsSizes(energy->dmv,&mi[0],&mi[1],&mi[2]);CHKERRQ(ierr);
ptatin3d_energyfv.c:130:    if (mi[d] != fv_mi[d]) SETERRQ1(energy->fv->comm,PETSC_ERR_USER,"DMDA for FV has inconsistent number of elements (direction %D)",d);
ptatin3d_energyfv.c:134:  Mi[0] = pctx->mx * energy->nsubdivision[0];
ptatin3d_energyfv.c:135:  Mi[1] = pctx->my * energy->nsubdivision[1];
ptatin3d_energyfv.c:136:  Mi[2] = pctx->mz * energy->nsubdivision[2];
ptatin3d_energyfv.c:138:  ierr = FVDASetSizes(energy->fv,mi,Mi);CHKERRQ(ierr);
ptatin3d_energyfv.c:140:  ierr = FVDASetProblemType(energy->fv,PETSC_TRUE,FVDA_PARABOLIC,0,0);CHKERRQ(ierr);
ptatin3d_energyfv.c:143:  ierr = DMGetCoordinateDM(energy->dmv,&fv_dmgeom);CHKERRQ(ierr);
ptatin3d_energyfv.c:144:  ierr = FVDASetGeometryDM(energy->fv,fv_dmgeom);CHKERRQ(ierr);
ptatin3d_energyfv.c:148:    ierr = DMGetCoordinates(energy->dmv,&gcoor);CHKERRQ(ierr);
ptatin3d_energyfv.c:150:    energy->fv->vertex_coor_geometry = gcoor;
ptatin3d_energyfv.c:154:  ierr = FVDASetUp(energy->fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:156:  //ierr = FVDASetup_TimeDep(energy->fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:157:  ierr = FVDASetup_ALE(energy->fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:160:  ierr = FVDARegisterFaceProperty(energy->fv,"v",3);CHKERRQ(ierr);
ptatin3d_energyfv.c:161:  ierr = FVDARegisterFaceProperty(energy->fv,"xDot",3);CHKERRQ(ierr);
ptatin3d_energyfv.c:163:  ierr = FVDARegisterFaceProperty(energy->fv,"v.n",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:164:  ierr = FVDARegisterFaceProperty(energy->fv,"xDot.n",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:165:  ierr = FVDARegisterFaceProperty(energy->fv,"k",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:167:  ierr = FVDARegisterCellProperty(energy->fv,"rho.cp",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:168:  ierr = FVDARegisterCellProperty(energy->fv,"k",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:169:  ierr = FVDARegisterCellProperty(energy->fv,"H",1);CHKERRQ(ierr);
ptatin3d_energyfv.c:177:    ierr = PetscCalloc1(3*energy->npoints_macro,&energy->xi_macro);CHKERRQ(ierr);
ptatin3d_energyfv.c:179:      dxi[d] = 2.0 / ((PetscReal)energy->nsubdivision[d]);
ptatin3d_energyfv.c:181:    for (kk=0; kk<energy->nsubdivision[2]+1; kk++) {
ptatin3d_energyfv.c:182:      for (jj=0; jj<energy->nsubdivision[1]+1; jj++) {
ptatin3d_energyfv.c:183:        for (ii=0; ii<energy->nsubdivision[0]+1; ii++) {
ptatin3d_energyfv.c:184:          energy->xi_macro[3*cnt+0] = -1.0 + ii * dxi[0];
ptatin3d_energyfv.c:185:          energy->xi_macro[3*cnt+1] = -1.0 + jj * dxi[1];
ptatin3d_energyfv.c:186:          energy->xi_macro[3*cnt+2] = -1.0 + kk * dxi[2];
ptatin3d_energyfv.c:192:    ierr = PetscCalloc1(energy->npoints_macro,&energy->basis_macro);CHKERRQ(ierr);
ptatin3d_energyfv.c:193:    for (d=0; d<energy->npoints_macro; d++) {
ptatin3d_energyfv.c:194:      ierr = PetscCalloc1(Q2_NODES_PER_EL_3D,&energy->basis_macro[d]);CHKERRQ(ierr);
ptatin3d_energyfv.c:196:    for (d=0; d<energy->npoints_macro; d++) {
ptatin3d_energyfv.c:197:      P3D_ConstructNi_Q2_3D(&energy->xi_macro[3*d],energy->basis_macro[d]);CHKERRQ(ierr);
ptatin3d_energyfv.c:201:  ierr = DMCreateGlobalVector(energy->fv->dm_fv,&energy->T);CHKERRQ(ierr);
ptatin3d_energyfv.c:202:  ierr = PetscObjectSetName((PetscObject)energy->T,"T");CHKERRQ(ierr);
ptatin3d_energyfv.c:203:  ierr = DMCreateGlobalVector(energy->fv->dm_fv,&energy->Told);CHKERRQ(ierr);
ptatin3d_energyfv.c:204:  ierr = DMCreateGlobalVector(energy->fv->dm_fv,&energy->F);CHKERRQ(ierr);
ptatin3d_energyfv.c:205:  ierr = DMCreateGlobalVector(energy->fv->dm_fv,&energy->G);CHKERRQ(ierr);
ptatin3d_energyfv.c:206:  /*ierr = DMCreateMatrix(energy->fv->dm_fv,&energy->J);CHKERRQ(ierr);*/
ptatin3d_energyfv.c:207:  ierr = FVDACreateMatrix(energy->fv,DMDA_STENCIL_STAR,&energy->J);CHKERRQ(ierr);
ptatin3d_energyfv.c:209:    ierr = DMCreateGlobalVector(energy->dmv,&energy->velocity);CHKERRQ(ierr);
ptatin3d_energyfv.c:210:    ierr = DMCreateGlobalVector(fv_dmgeom,&energy->Xold);CHKERRQ(ierr);
ptatin3d_energyfv.c:214:  ierr = SNESCreate(energy->fv->comm,&energy->snes);CHKERRQ(ierr);
ptatin3d_energyfv.c:215:  ierr = SNESSetOptionsPrefix(energy->snes,"energyfv_");CHKERRQ(ierr);
ptatin3d_energyfv.c:216:  ierr = SNESSetDM(energy->snes,energy->fv->dm_fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:217:  ierr = SNESSetSolution(energy->snes,energy->T);CHKERRQ(ierr);
ptatin3d_energyfv.c:218:  ierr = SNESSetApplicationContext(energy->snes,(void*)energy->fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:219:  //ierr = SNESSetApplicationContext(energy->snes,(void*)energy);CHKERRQ(ierr);
ptatin3d_energyfv.c:221:  //ierr = SNESSetFunction(energy->snes,energy->F,          fvda_eval_F_timedep,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:222:  //ierr = SNESSetFunction(energy->snes,energy->F,          fvda_highres_eval_F_timedep,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:223:  //ierr = SNESSetJacobian(energy->snes,energy->J,energy->J,fvda_eval_J_timedep,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:226:  //ierr = SNESSetFunction(energy->snes,energy->F,          fvda_eval_F_forward_ale,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:227:  //ierr = SNESSetJacobian(energy->snes,energy->J,energy->J,fvda_eval_J_forward_ale,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:230:  ierr = SNESSetFunction(energy->snes,energy->F,          fvda_highres_eval_F_forward_ale,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:231:  ierr = SNESSetJacobian(energy->snes,energy->J,energy->J,fvda_eval_J_forward_ale,NULL);CHKERRQ(ierr);
ptatin3d_energyfv.c:234:  ierr = SNESSetFromOptions(energy->snes);CHKERRQ(ierr);
ptatin3d_energyfv.c:259:  dmv_fv = energy->dmv;
ptatin3d_energyfv.c:265:  dmc_fv = energy->fv->dm_geometry;
ptatin3d_energyfv.c:269:  ierr = PetscCalloc1(energy->npoints_macro,&fvindex);CHKERRQ(ierr);
ptatin3d_energyfv.c:281:    ierr = _cart_convert_index_to_ijk(e,(const PetscInt*)energy->mi_parent,q2_start);CHKERRQ(ierr);
ptatin3d_energyfv.c:284:    fv_start[0] = q2_start[0] * energy->nsubdivision[0];
ptatin3d_energyfv.c:285:    fv_start[1] = q2_start[1] * energy->nsubdivision[1];
ptatin3d_energyfv.c:286:    fv_start[2] = q2_start[2] * energy->nsubdivision[2];
ptatin3d_energyfv.c:288:    ierr = _cart_convert_ijk_to_index((const PetscInt*)fv_start,(const PetscInt*)energy->fv->mi,&fv_cell);CHKERRQ(ierr);
ptatin3d_energyfv.c:294:    for (kk=0; kk<energy->nsubdivision[2]+1; kk++) {
ptatin3d_energyfv.c:295:      for (jj=0; jj<energy->nsubdivision[1]+1; jj++) {
ptatin3d_energyfv.c:296:        for (ii=0; ii<energy->nsubdivision[0]+1; ii++) {
ptatin3d_energyfv.c:303:    for (m=0; m<energy->npoints_macro; m++) {
ptatin3d_energyfv.c:309:        x[0] += energy->basis_macro[m][k] * elcoords[3*k+0];
ptatin3d_energyfv.c:310:        x[1] += energy->basis_macro[m][k] * elcoords[3*k+1];
ptatin3d_energyfv.c:311:        x[2] += energy->basis_macro[m][k] * elcoords[3*k+2];
ptatin3d_energyfv.c:325:  ierr = DMLocalToGlobal(dmc_fv,coor_fv,INSERT_VALUES,energy->fv->vertex_coor_geometry);CHKERRQ(ierr);
ptatin3d_energyfv.c:327:  ierr = FVDAUpdateGeometry(energy->fv);CHKERRQ(ierr);
ptatin3d_energyfv.c:349:  if (bs > 3) SETERRQ(energy->fv->comm,PETSC_ERR_SUP,"Blocksize > 3 not supported - requires a small patch/fix to be applied (trivial)");
ptatin3d_energyfv.c:358:  ierr = PetscCalloc1(energy->npoints_macro,&fvindex);CHKERRQ(ierr);
ptatin3d_energyfv.c:375:    ierr = _cart_convert_index_to_ijk(e,(const PetscInt*)energy->mi_parent,q2_start);CHKERRQ(ierr);
ptatin3d_energyfv.c:378:    fv_start[0] = q2_start[0] * energy->nsubdivision[0];
ptatin3d_energyfv.c:379:    fv_start[1] = q2_start[1] * energy->nsubdivision[1];
ptatin3d_energyfv.c:380:    fv_start[2] = q2_start[2] * energy->nsubdivision[2];
ptatin3d_energyfv.c:382:    ierr = _cart_convert_ijk_to_index((const PetscInt*)fv_start,(const PetscInt*)energy->fv->mi,&fv_cell);CHKERRQ(ierr);
ptatin3d_energyfv.c:388:    for (kk=0; kk<energy->nsubdivision[2]+1; kk++) {
ptatin3d_energyfv.c:389:      for (jj=0; jj<energy->nsubdivision[1]+1; jj++) {
ptatin3d_energyfv.c:390:        for (ii=0; ii<energy->nsubdivision[0]+1; ii++) {
ptatin3d_energyfv.c:397:    for (m=0; m<energy->npoints_macro; m++) {
ptatin3d_energyfv.c:406:          field[b] += energy->basis_macro[m][k] * elfield[bs*k+b];
ptatin3d_energyfv.c:459:  fv = energy->fv;
ptatin3d_energyfv.c:460:  dm = energy->dmv;
ptatin3d_energyfv.c:522:  fv = energy->fv;
ptatin3d_energyfv.c:523:  dm = energy->dmv;
ptatin3d_energyfv.c:770:    //ierr = FVDAView_CellData(energy->fv,energy->T,PETSC_TRUE,"xcell");CHKERRQ(ierr);
ptatin3d_energyfv.c:989:        conductivity_mp = mpp_energy->diffusivity;
ptatin3d_energyfv.c:1029:          H_mp += mpp_energy->heat_source;
ptatin3d_energyfv.c:1520:  fv = energy->fv;
ptatin3d_energyfv.c:1521:  dm = energy->dmv;
ptatin3d_energyfv.c:1527:  ierr = DMGetLocalVector(energy->fv->dm_geometry,&geometry_coorl);CHKERRQ(ierr);
ptatin3d_energyfv.c:1528:  ierr = DMGlobalToLocal(energy->fv->dm_geometry,energy->fv->vertex_coor_geometry,INSERT_VALUES,geometry_coorl);CHKERRQ(ierr);
ptatin3d_energyfv.c:1531:  ierr = FVDAGetCellPropertyByNameArray(energy->fv,"k",&k_cell);CHKERRQ(ierr);
