| /petsc/src/dm/impls/plex/transform/impls/extrude/ |
| H A D | plextrcohesive.c | 32 DMPlexTransform_Cohesive *ex = (DMPlexTransform_Cohesive *)tr->data; in DMPlexTransformView_Cohesive() local 44 …PetscCall(PetscViewerASCIIPrintf(viewer, " create tensor cells: %s\n", ex->useTensor ? "YES" : "N… in DMPlexTransformView_Cohesive() 53 DMPlexTransform_Cohesive *ex = (DMPlexTransform_Cohesive *)tr->data; in DMPlexTransformSetFromOptions_Cohesive() local 59 …PetscCall(PetscOptionsBool("-dm_plex_transform_extrude_use_tensor", "Create tensor cells", "", ex-… in DMPlexTransformSetFromOptions_Cohesive() 61 …sReal("-dm_plex_transform_cohesive_width", "Width of a cohesive cell", "", ex->width, &width, &flg… in DMPlexTransformSetFromOptions_Cohesive() 63 …sInt("-dm_plex_transform_cohesive_debug", "Det debugging level", "", ex->debug, &ex->debug, NULL)); in DMPlexTransformSetFromOptions_Cohesive() 144 static PetscErrorCode DMPlexTransformCohesiveExtrudeSetUp_Point(DMPlexTransform_Cohesive *ex) in DMPlexTransformCohesiveExtrudeSetUp_Point() argument 151 ex->Nt[rt] = 2; in DMPlexTransformCohesiveExtrudeSetUp_Point() 154 …PetscCall(PetscMalloc4(ex->Nt[rt], &ex->target[rt], ex->Nt[rt], &ex->size[rt], Nc, &ex->cone[rt], … in DMPlexTransformCohesiveExtrudeSetUp_Point() 155 ex->target[rt][0] = DM_POLYTOPE_POINT; in DMPlexTransformCohesiveExtrudeSetUp_Point() [all …]
|
| H A D | plextrextrude.c | 7 DMPlexTransform_Extrude *ex = (DMPlexTransform_Extrude *)tr->data; in DMPlexTransformView_Extrude() local 26 … PetscCall(PetscViewerASCIIPrintf(viewer, " number of layers: %" PetscInt_FMT "\n", ex->layers)); in DMPlexTransformView_Extrude() 27 …PetscCall(PetscViewerASCIIPrintf(viewer, " create tensor cells: %s\n", ex->useTensor ? "YES" : "N… in DMPlexTransformView_Extrude() 28 if (ex->periodic) PetscCall(PetscViewerASCIIPrintf(viewer, " periodic\n")); in DMPlexTransformView_Extrude() 29 …PetscCall(PetscViewerASCIIPrintf(viewer, " normal algorithm: %s\n", PlexNormalAlgs[ex->normalAlg]… in DMPlexTransformView_Extrude() 30 …if (ex->normalFunc) PetscCall(PetscViewerASCIIPrintf(viewer, " normal modified by user function\n… in DMPlexTransformView_Extrude() 39 DMPlexTransform_Extrude *ex = (DMPlexTransform_Extrude *)tr->data; in DMPlexTransformSetFromOptions_Extrude() local 47 …nt("-dm_plex_transform_extrude_layers", "Number of layers to extrude", "", ex->layers, &nl, &flg, … in DMPlexTransformSetFromOptions_Extrude() 49 …ex_transform_extrude_thickness", "Total thickness of extruded layers", "", ex->thickness, &th, &fl… in DMPlexTransformSetFromOptions_Extrude() 51 …PetscCall(PetscOptionsBool("-dm_plex_transform_extrude_use_tensor", "Create tensor cells", "", ex-… in DMPlexTransformSetFromOptions_Extrude() [all …]
|
| /petsc/src/dm/impls/da/hypre/ |
| H A D | mhyp.c | 30 Mat_HYPREStruct *ex = (Mat_HYPREStruct *)mat->data; in MatSetValuesLocal_HYPREStruct_3d() local 43 } else if (stencil == -ex->gnx) { in MatSetValuesLocal_HYPREStruct_3d() 45 } else if (stencil == ex->gnx) { in MatSetValuesLocal_HYPREStruct_3d() 47 } else if (stencil == -ex->gnxgny) { in MatSetValuesLocal_HYPREStruct_3d() 49 } else if (stencil == ex->gnxgny) { in MatSetValuesLocal_HYPREStruct_3d() 53 row = ex->gindices[irow[i]] - ex->rstart; in MatSetValuesLocal_HYPREStruct_3d() 54 index[0] = (HYPRE_Int)(ex->xs + (row % ex->nx)); in MatSetValuesLocal_HYPREStruct_3d() 55 index[1] = (HYPRE_Int)(ex->ys + ((row / ex->nx) % ex->ny)); in MatSetValuesLocal_HYPREStruct_3d() 56 index[2] = (HYPRE_Int)(ex->zs + (row / (ex->nxny))); in MatSetValuesLocal_HYPREStruct_3d() 57 …if (addv == ADD_VALUES) PetscCallHYPRE(HYPRE_StructMatrixAddToValues(ex->hmat, index, (HYPRE_Int)n… in MatSetValuesLocal_HYPREStruct_3d() [all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex8.c | 20 PetscInt startx, starty, startz, nx, ny, nz, ex, ey, ez, nExtrax, nExtray, nExtraz; in main() local 54 for (ex = startx; ex < startx + nx + nExtrax; ++ex) { in main() 55 arrSol[ez][ey][ex][iux] = 2 * PRESSURE_CONST; in main() 56 arrRHS[ez][ey][ex][iux] = 0.0; in main() 57 arrSol[ez][ey][ex][iuy] = 2 * PRESSURE_CONST; in main() 58 arrRHS[ez][ey][ex][iuy] = 0.0; in main() 59 arrSol[ez][ey][ex][iuz] = 2 * PRESSURE_CONST; in main() 60 arrRHS[ez][ey][ex][iuz] = 0.0; in main() 61 if (ex < startx + nx && ey < starty + ny && ez < startz + nz) { in main() 62 arrSol[ez][ey][ex][ip] = PRESSURE_CONST; in main() [all …]
|
| H A D | ex30.c | 66 PetscInt startx, starty, startz, N[3], nx, ny, nz, ex, ey, ez, d; in CreateMat() local 99 for (ex = startx; ex < startx + nx; ++ex) { in CreateMat() 100 if (ex == N[0] - 1) { in CreateMat() 104 row.i = ex; in CreateMat() 115 row.i = ex; in CreateMat() 126 row.i = ex; in CreateMat() 135 if (ex == 0) { in CreateMat() 139 row.i = ex; in CreateMat() 151 row.i = ex; in CreateMat() 159 col[0].i = ex; in CreateMat() [all …]
|
| H A D | ex11.c | 20 PetscInt startx, starty, nx, ny, ex, ey, nExtrax, nExtray; in main() local 52 for (ex = startx; ex < startx + nx + nExtrax; ++ex) { in main() 53 arrSol[ey][ex][iux] = 2 * PRESSURE_CONST; in main() 54 arrRHS[ey][ex][iux] = 0.0; in main() 55 arrSol[ey][ex][iuy] = 2 * PRESSURE_CONST; in main() 56 arrRHS[ey][ex][iuy] = 0.0; in main() 57 if (ex < startx + nx && ey < starty + ny) { in main() 58 arrSol[ey][ex][ip] = PRESSURE_CONST; in main() 59 arrRHS[ey][ex][ip] = PRESSURE_CONST; in main() 120 …PetscInt startx, starty, nx, ny, nExtrax, nExtray, ex, ey, idxP, idxUx, idxUy, startGhostx, … in ApplyOperator() local [all …]
|
| H A D | ex23.c | 9 PetscInt ex, ey, ez, n[3], start[3], nExtra[3], iNext, iPrev, iCenter, d, round; in main() local 29 for (ex = start[0]; ex < start[0] + n[0]; ++ex) { in main() 30 cArrX[ex][iPrev] *= 1.1; in main() 31 cArrX[ex][iNext] = cArrX[ex][iPrev] + 0.1; in main() 32 cArrX[ex][iCenter] = 0.5 * (cArrX[ex][iPrev] + cArrX[ex][iNext]); in main() 46 for (ex = start[0]; ex < start[0] + n[0]; ++ex) { in main() 47 cArrX[ex][iPrev] = 2.0; in main() 48 cArrX[ex][iNext] = 2.0; in main() 49 cArrX[ex][iCenter] = 2.0; in main()
|
| H A D | ex15.c | 244 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem2d() local 259 for (ex = startx; ex < startx + nx; ++ex) { in CreateSystem2d() 260 if (ex == N[0] - 1) { in CreateSystem2d() 265 row.i = ex; in CreateSystem2d() 278 row.i = ex; in CreateSystem2d() 292 row.i = ex; in CreateSystem2d() 305 row.i = ex; in CreateSystem2d() 309 if (ex == 0) { in CreateSystem2d() 311 col[0].i = ex; in CreateSystem2d() 316 col[1].i = ex; in CreateSystem2d() [all …]
|
| H A D | ex18.c | 252 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem() local 292 for (ex = startx; ex < startx + nx; ++ex) { in CreateSystem() 293 if (ex == N[0] - 1) { in CreateSystem() 299 row.i = ex; in CreateSystem() 305 valRhs = bogusScale * uxRef(cArrX[ex][inext], cArrY[ey][icenter]); in CreateSystem() 314 row.i = ex; in CreateSystem() 320 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][inext]); in CreateSystem() 330 row.i = ex; in CreateSystem() 336 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem() 345 row.i = ex; in CreateSystem() [all …]
|
| H A D | ex40.c | 359 for (PetscInt ex = start[0]; ex < start[0] + n[0] + n_extra[0]; ++ex) { in FormFunction2DNoCoupling() local 364 row.i = ex; in FormFunction2DNoCoupling() 372 if (ex < N[0]) { in FormFunction2DNoCoupling() 377 row.i = ex; in FormFunction2DNoCoupling() 391 row.i = ex; in FormFunction2DNoCoupling() 400 if (ex < N[0] && ey < N[1]) { in FormFunction2DNoCoupling() 405 row.i = ex; in FormFunction2DNoCoupling() 437 for (PetscInt ex = start[0]; ex < start[0] + n[0] + n_extra[0]; ++ex) { in FormJacobian2DNoCoupling() local 442 row.i = ex; in FormJacobian2DNoCoupling() 450 if (ex < N[0]) { in FormJacobian2DNoCoupling() [all …]
|
| H A D | ex5.c | 114 PetscInt start, n, nExtra, ex, idxP, idxU, startGhost, nGhost; in ApplyOperator() local 140 for (ex = start; ex < start + n + nExtra; ++ex) { in ApplyOperator() 141 if (ex < start + n) { /* Don't compute pressure outside domain */ in ApplyOperator() 142 arrOut[ex][idxP] = arrIn[ex][idxP]; in ApplyOperator() 144 arrOut[ex][idxU] = arrIn[ex][idxP] + arrIn[ex - 1][idxP] - arrIn[ex][idxU]; in ApplyOperator()
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex3.c | 170 PetscInt startx, starty, startz, N[3], nx, ny, nz, ex, ey, ez, d; in CreateSystem() local 205 for (ex = startx; ex < startx + nx; ++ex) { in CreateSystem() 206 if (ex == N[0] - 1) { in CreateSystem() 211 row.i = ex; in CreateSystem() 217 …valRhs = uxRef(arrCoord[ez][ey][ex][icux_right[0]], arrCoord[ez][ey][ex][icux_right[1]], arrCoord[… in CreateSystem() 225 row.i = ex; in CreateSystem() 231 …valRhs = uyRef(arrCoord[ez][ey][ex][icuy_up[0]], arrCoord[ez][ey][ex][icuy_up[1]], arrCoord[ez][ey… in CreateSystem() 239 row.i = ex; in CreateSystem() 245 …valRhs = uzRef(arrCoord[ez][ey][ex][icuz_front[0]], arrCoord[ez][ey][ex][icuz_front[1]], arrCoord[… in CreateSystem() 250 if (ex == 0) { in CreateSystem() [all …]
|
| H A D | ex6.c | 228 PetscInt N[3], ex, ey, ez, startx, starty, startz, nx, ny, nz, extrax, extray, extraz; in CreateLame() local 237 for (ex = startx; ex < startx + nx; ++ex) { in CreateLame() 240 pos.i = ex; in CreateLame() 245 pos.i = ex; in CreateLame() 254 for (ex = startx; ex < startx + nx + extrax; ++ex) { in CreateLame() 257 pos.i = ex; in CreateLame() 268 for (ex = startx; ex < startx + nx; ++ex) { in CreateLame() 271 pos.i = ex; in CreateLame() 277 pos.i = ex; in CreateLame() 289 for (ex = startx; ex < startx + nx + extrax; ++ex) { in CreateLame() [all …]
|
| H A D | ex2.c | 148 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem() local 181 for (ex = startx; ex < startx + nx; ++ex) { in CreateSystem() 182 if (ex == N[0] - 1) { in CreateSystem() 187 row.i = ex; in CreateSystem() 192 valRhs = uxRef(cArrX[ex][inext], cArrY[ey][icenter]); in CreateSystem() 200 row.i = ex; in CreateSystem() 205 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][inext]); in CreateSystem() 214 row.i = ex; in CreateSystem() 219 valRhs = uyRef(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem() 227 row.i = ex; in CreateSystem() [all …]
|
| H A D | ex4.c | 576 PetscInt ex, ey, startx, starty, nx, ny; in CreateSystem2d() local 616 for (ex = startx; ex < startx + nx; ++ex) { in CreateSystem2d() 617 const PetscBool left_boundary = (PetscBool)(ex == 0); in CreateSystem2d() 618 const PetscBool right_boundary = (PetscBool)(ex == N[0] - 1); in CreateSystem2d() 628 row.i = ex; in CreateSystem2d() 642 row.i = ex; in CreateSystem2d() 659 row.i = ex; in CreateSystem2d() 665 rhoPoint[0].i = ex; in CreateSystem2d() 669 rhoPoint[1].i = ex; in CreateSystem2d() 677 etaPoint[0].i = ex; in CreateSystem2d() [all …]
|
| /petsc/ |
| H A D | .gitignore | 51 ex[0-9] 52 ex[0-9]f 53 ex[0-9]dm 54 ex[0-9]d* 55 ex[0-9]f90 56 ex[0-9]f90t 57 ex[0-9][0-9] 58 ex[0-9][0-9]f 59 ex[0-9][0-9]dm 60 ex[0-9][0-9]d* [all …]
|
| /petsc/src/dm/impls/stag/ |
| H A D | stagda.c | 330 PetscInt ex; in DMStagTransferCoordinatesToDMDA() local 338 …for (ex = start[0]; ex < start[0] + n[0] + extraPoint[0]; ++ex) cArrDa[ex][0] = cArrStag[ex][slot]; in DMStagTransferCoordinatesToDMDA() 343 for (ex = start[0]; ex < start[0] + n[0] + extraPoint[0]; ++ex) cArrDa[ex][0] = cArrX[ex][0]; in DMStagTransferCoordinatesToDMDA() 348 PetscInt ex, ey; in DMStagTransferCoordinatesToDMDA() local 357 for (ex = start[0]; ex < start[0] + n[0] + extraPoint[0]; ++ex) { in DMStagTransferCoordinatesToDMDA() 358 for (d = 0; d < 2; ++d) cArrDa[ey][ex][d] = cArrStag[ey][ex][slot + d]; in DMStagTransferCoordinatesToDMDA() 366 for (ex = start[0]; ex < start[0] + n[0] + extraPoint[0]; ++ex) { in DMStagTransferCoordinatesToDMDA() 367 cArrDa[ey][ex][0] = cArrX[ex][0]; in DMStagTransferCoordinatesToDMDA() 368 cArrDa[ey][ex][1] = cArrY[ey][0]; in DMStagTransferCoordinatesToDMDA() 375 PetscInt ex, ey, ez; in DMStagTransferCoordinatesToDMDA() local [all …]
|
| /petsc/lib/petsc/conf/ |
| H A D | rules | 16 ex[0-9] ex[0-9][0-9] ex[0-9][0-9][0-9] \ 17 ex[0-9]f ex[0-9][0-9]f ex[0-9][0-9][0-9]f \ 18 ex[0-9]k ex[0-9][0-9]k ex[0-9][0-9][0-9]k \ 19 ex[0-9]f90 ex[0-9][0-9]f90 ex[0-9][0-9][0-9]f90 \ 20 ex[0-9]cu ex[0-9][0-9]cu ex[0-9][0-9][0-9]cu \ 21 ex[0-9].exe ex[0-9][0-9].exe ex[0-9][0-9][0-9].exe \ 22 ex[0-9]f.exe ex[0-9][0-9]f.exe ex[0-9][0-9][0-9]f.exe \ 23 ex[0-9]f90.exe ex[0-9][0-9]f90.exe ex[0-9][0-9][0-9]f90.exe \ 24 ex[0-9]cu.exe ex[0-9][0-9]cu.exe ex[0-9][0-9][0-9]cu.exe \ 25 ex[0-9]hip ex[0-9][0-9]hip ex[0-9][0-9][0-9]hip \ [all …]
|
| /petsc/src/vec/vec/impls/seq/seqviennacl/ |
| H A D | vecviennacl.cxx | 102 } catch (std::exception const &ex) { in PetscViennaCLInit() local 103 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex.what()); in PetscViennaCLInit() 120 } catch (std::exception const &ex) { in PetscViennaCLInit() local 121 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex.what()); in PetscViennaCLInit() 198 } catch (std::exception const &ex) { in VecViennaCLAllocateCheck() local 199 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex.what()); in VecViennaCLAllocateCheck() 218 } catch (std::exception const &ex) { in VecViennaCLCopyToGPU() local 219 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex.what()); in VecViennaCLCopyToGPU() 243 } catch (std::exception const &ex) { in VecViennaCLCopyFromGPU() local 244 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex.what()); in VecViennaCLCopyFromGPU() [all …]
|
| /petsc/src/ksp/pc/impls/hypre/ |
| H A D | hypre.c | 2737 PC_PFMG *ex = (PC_PFMG *)pc->data; in PCDestroy_PFMG() local 2740 if (ex->hsolver) PetscCallHYPRE(HYPRE_StructPFMGDestroy(ex->hsolver)); in PCDestroy_PFMG() 2741 PetscCall(PetscCommRestoreComm(PetscObjectComm((PetscObject)pc), &ex->hcomm)); in PCDestroy_PFMG() 2752 PC_PFMG *ex = (PC_PFMG *)pc->data; in PCView_PFMG() local 2758 PetscCall(PetscViewerASCIIPrintf(viewer, " max iterations %" PetscInt_FMT "\n", ex->its)); in PCView_PFMG() 2759 PetscCall(PetscViewerASCIIPrintf(viewer, " tolerance %g\n", ex->tol)); in PCView_PFMG() 2760 PetscCall(PetscViewerASCIIPrintf(viewer, " relax type %s\n", PFMGRelaxType[ex->relax_type])); in PCView_PFMG() 2761 PetscCall(PetscViewerASCIIPrintf(viewer, " RAP type %s\n", PFMGRAPType[ex->rap_type])); in PCView_PFMG() 2762 …mber pre-relax %" PetscInt_FMT " post-relax %" PetscInt_FMT "\n", ex->num_pre_relax, ex->num_post_… in PCView_PFMG() 2763 PetscCall(PetscViewerASCIIPrintf(viewer, " max levels %" PetscInt_FMT "\n", ex->max_levels)); in PCView_PFMG() [all …]
|
| /petsc/src/ksp/pc/impls/chowiluviennacl/ |
| H A D | chowiluviennacl.cxx | 46 } catch (char *ex) { in PCSetUp_CHOWILUVIENNACL() local 47 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_CHOWILUVIENNACL() 62 } catch (char *ex) { in PCSetUp_CHOWILUVIENNACL() local 63 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_CHOWILUVIENNACL() 103 } catch (char *ex) { in PCApply_CHOWILUVIENNACL() local 104 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCApply_CHOWILUVIENNACL() 129 } catch (char *ex) { in PCDestroy_CHOWILUVIENNACL() local 130 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCDestroy_CHOWILUVIENNACL()
|
| /petsc/src/ksp/pc/impls/rowscalingviennacl/ |
| H A D | rowscalingviennacl.cxx | 46 } catch (char *ex) { in PCSetUp_ROWSCALINGVIENNACL() local 47 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_ROWSCALINGVIENNACL() 62 } catch (char *ex) { in PCSetUp_ROWSCALINGVIENNACL() local 63 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_ROWSCALINGVIENNACL() 103 } catch (char *ex) { in PCApply_ROWSCALINGVIENNACL() local 104 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCApply_ROWSCALINGVIENNACL() 129 } catch (char *ex) { in PCDestroy_ROWSCALINGVIENNACL() local 130 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCDestroy_ROWSCALINGVIENNACL()
|
| /petsc/src/ksp/pc/impls/saviennacl/ |
| H A D | saviennacl.cxx | 45 } catch (char *ex) { in PCSetUp_SAVIENNACL() local 46 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_SAVIENNACL() 64 } catch (char *ex) { in PCSetUp_SAVIENNACL() local 65 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCSetUp_SAVIENNACL() 102 } catch (char *ex) { in PCApply_SAVIENNACL() local 103 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCApply_SAVIENNACL() 128 } catch (char *ex) { in PCDestroy_SAVIENNACL() local 129 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_LIB, "ViennaCL error: %s", ex); in PCDestroy_SAVIENNACL()
|
| /petsc/src/ts/tutorials/output/ |
| H A D | fig.tex.out | 23 \draw (0.5,0.9) node [rectangle, draw, align=center, inner sep=1ex] {TS: pseudo}; 24 \draw (0.5,0.83) node [rectangle, draw, align=center, inner sep=1ex] {SNES: newtonls}; 25 \draw (0.5,0.76) node [rectangle, draw, align=center, inner sep=1ex] {KSP: gmres}; 26 \draw (0.5,0.69) node [rectangle, draw, align=center, inner sep=1ex] {PC: ilu (16)};
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 282 void DeformationGradient(Field *ex, PetscInt qi, PetscInt qj, PetscInt qk, PetscScalar *invJ, Petsc… in DeformationGradient() argument 297 F[0] += lgrad[0] * ex[idx][0]; in DeformationGradient() 298 F[1] += lgrad[1] * ex[idx][0]; in DeformationGradient() 299 F[2] += lgrad[2] * ex[idx][0]; in DeformationGradient() 300 F[3] += lgrad[0] * ex[idx][1]; in DeformationGradient() 301 F[4] += lgrad[1] * ex[idx][1]; in DeformationGradient() 302 F[5] += lgrad[2] * ex[idx][1]; in DeformationGradient() 303 F[6] += lgrad[0] * ex[idx][2]; in DeformationGradient() 304 F[7] += lgrad[1] * ex[idx][2]; in DeformationGradient() 305 F[8] += lgrad[2] * ex[idx][2]; in DeformationGradient() [all …]
|