Lines Matching refs:ierr
32 PetscErrorCode :: ierr
98 PetscCallA(PetscInitialize(PETSC_NULL_CHARACTER, ierr))
99 if (ierr /= 0) then
104 PetscCallMPIA(MPI_Comm_rank(PETSC_COMM_WORLD, rank, ierr))
105 PetscCallMPIA(MPI_Comm_size(PETSC_COMM_WORLD, numProc, ierr))
106 …cCallA(PetscOptionsGetString(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-i', ifilename, flg, ierr))
108 …cCallA(PetscOptionsGetString(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-o', ofilename, flg, ierr))
110 …etscCallA(PetscOptionsGetInt(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-order', order, flg, ierr))
117 …allA(DMPlexCreateFromFile(PETSC_COMM_WORLD, ifilename, PETSC_NULL_CHARACTER, PETSC_TRUE, dm, ierr))
118 PetscCallA(DMPlexDistributeSetDefault(dm, PETSC_FALSE, ierr))
119 PetscCallA(DMSetFromOptions(dm, ierr))
120 PetscCallA(DMGetDimension(dm, sdim, ierr))
121 PetscCallA(DMViewFromOptions(dm, PETSC_NULL_OBJECT, '-dm_view', ierr))
126 PetscCallA(exopts(EXVRBS + EXDEBG, ierr))
128 PetscCallA(PetscViewerExodusIIOpen(PETSC_COMM_WORLD, ofilename, FILE_MODE_WRITE, viewer, ierr))
137 PetscCallA(PetscViewerExodusIISetOrder(viewer, order, ierr))
138 PetscCallA(PetscViewerView(viewer, PETSC_VIEWER_STDOUT_WORLD, ierr))
146 PetscCallA(DMView(dm, viewer, ierr))
147 PetscCallA(PetscViewerView(viewer, PETSC_VIEWER_STDOUT_WORLD, ierr))
156 PetscCallA(PetscViewerExodusIISetZonalVariable(viewer, numZonalVar, ierr))
157 PetscCallA(PetscViewerExodusIISetNodalVariable(viewer, numNodalVar, ierr))
159 PetscCallA(PetscViewerExodusIISetZonalVariableName(viewer, i - 1, zonalVarName2D(i), ierr))
162 PetscCallA(PetscViewerExodusIISetNodalVariableName(viewer, i - 1, nodalVarName2D(i), ierr))
167 PetscCallA(PetscViewerExodusIISetZonalVariable(viewer, numZonalVar, ierr))
168 PetscCallA(PetscViewerExodusIISetNodalVariable(viewer, numNodalVar, ierr))
170 PetscCallA(PetscViewerExodusIISetZonalVariableName(viewer, i - 1, zonalVarName3D(i), ierr))
173 PetscCallA(PetscViewerExodusIISetNodalVariableName(viewer, i - 1, nodalVarName3D(i), ierr))
178 PetscCallA(PetscViewerView(viewer, PETSC_VIEWER_STDOUT_WORLD, ierr))
182 PetscCallA(PetscViewerExodusIIGetId(viewer, exoid, ierr))
183 PetscCallA(exinq(exoid, EX_INQ_ELEM_BLK, numCS, PETSC_NULL_REAL, sjunk, ierr))
186 PetscCallA(expvtt(exoid, numCS, numZonalVar, truthtable, ierr))
191 PetscCallA(exptim(exoid, step, real(step, kind=kPR), ierr))
194 PetscCallA(PetscObjectGetComm(dm, comm, ierr))
195 PetscCallA(PetscSectionCreate(comm, section, ierr))
196 PetscCallA(PetscSectionSetNumFields(section, 3_kPI, ierr))
197 PetscCallA(PetscSectionSetFieldName(section, fieldU, 'U', ierr))
198 PetscCallA(PetscSectionSetFieldName(section, fieldA, 'Alpha', ierr))
199 PetscCallA(PetscSectionSetFieldName(section, fieldS, 'Sigma', ierr))
200 PetscCallA(DMPlexGetChart(dm, pStart, pEnd, ierr))
201 PetscCallA(PetscSectionSetChart(section, pStart, pEnd, ierr))
206 PetscCallA(DMPlexGetDepthStratum(dm, d - 1, pStartDepth(d), pEndDepth(d), ierr))
210 PetscCallA(PetscSectionSetFieldComponents(section, fieldU, sdim, ierr))
211 PetscCallA(PetscSectionSetFieldComponents(section, fieldA, 1_kPI, ierr))
212 PetscCallA(PetscSectionSetFieldComponents(section, fieldS, sdim*(sdim + 1)/2, ierr))
215 PetscCallA(DMGetLabelSize(dm, 'Cell Sets', numCS, ierr))
216 PetscCallA(DMGetLabelIdIS(dm, 'Cell Sets', csIS, ierr))
217 PetscCallA(ISGetIndices(csIS, csID, ierr))
223 PetscCallA(DMGetStratumSize(dm, 'Cell Sets', csID(set), numCells, ierr))
224 PetscCallA(DMGetStratumIS(dm, 'Cell Sets', csID(set), cellIS, ierr))
238 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
240 …tscCallA(DMPlexGetTransitiveClosure(dm, cellID(1), PETSC_TRUE, PETSC_NULL_INTEGER, closureA, ierr))
278 …allA(DMPlexRestoreTransitiveClosure(dm, cellID(1), PETSC_TRUE, PETSC_NULL_INTEGER, closureA, ierr))
281 …cCallA(DMPlexGetTransitiveClosure(dm, cellID(cell), PETSC_TRUE, PETSC_NULL_INTEGER, closure, ierr))
286 PetscCallA(PetscSectionSetDof(section, closure(p), dofU(d) + dofA(d) + dofS(d), ierr))
287 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldU, dofU(d), ierr))
288 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldA, dofA(d), ierr))
289 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldS, dofS(d), ierr))
293 …lA(DMPlexRestoreTransitiveClosure(dm, cellID(cell), PETSC_TRUE, PETSC_NULL_INTEGER, closure, ierr))
295 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
296 PetscCallA(ISDestroy(cellIS, ierr))
299 PetscCallA(ISRestoreIndices(csIS, csID, ierr))
300 PetscCallA(ISDestroy(csIS, ierr))
301 PetscCallA(PetscSectionSetUp(section, ierr))
302 PetscCallA(DMSetLocalSection(dm, section, ierr))
303 PetscCallA(PetscObjectViewFromOptions(section, PETSC_NULL_OBJECT, '-dm_section_view', ierr))
304 PetscCallA(PetscSectionGetValueLayout(PETSC_COMM_WORLD, section, layout, ierr))
305 PetscCallA(PetscLayoutDestroy(layout, ierr))
306 PetscCallA(PetscSectionDestroy(section, ierr))
308 PetscCallA(DMSetUseNatural(dm, PETSC_TRUE, ierr))
309 PetscCallA(DMPlexGetPartitioner(dm, part, ierr))
310 PetscCallA(PetscPartitionerSetFromOptions(part, ierr))
311 PetscCallA(DMPlexDistribute(dm, 0_kPI, migrationSF, pdm, ierr))
314 PetscCallA(DMPlexSetMigrationSF(pdm, migrationSF, ierr))
315 PetscCallA(PetscSFDestroy(migrationSF, ierr))
319 PetscCallA(DMViewFromOptions(pdm, PETSC_NULL_OBJECT, '-dm_view', ierr))
322 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldU], isU, dmU, ierr))
323 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldA], isA, dmA, ierr))
324 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldS], isS, dmS, ierr))
325 PetscCallA(DMCreateSubDM(pdm, 2_kPI, fieldUA, isUA, dmUA, ierr))
331 PetscCallA(DMCreateSuperDM(dmList, 2_kPI, PETSC_NULL_IS_POINTER, dmUA2, ierr))
334 PetscCallA(DMGetGlobalVector(pdm, X, ierr))
335 PetscCallA(DMGetGlobalVector(dmU, U, ierr))
336 PetscCallA(DMGetGlobalVector(dmA, A, ierr))
337 PetscCallA(DMGetGlobalVector(dmS, S, ierr))
338 PetscCallA(DMGetGlobalVector(dmUA, UA, ierr))
339 PetscCallA(DMGetGlobalVector(dmUA2, UA2, ierr))
341 PetscCallA(PetscObjectSetName(U, 'U', ierr))
342 PetscCallA(PetscObjectSetName(A, 'Alpha', ierr))
343 PetscCallA(PetscObjectSetName(S, 'Sigma', ierr))
344 PetscCallA(PetscObjectSetName(UA, 'UAlpha', ierr))
345 PetscCallA(PetscObjectSetName(UA2, 'UAlpha2', ierr))
346 PetscCallA(VecSet(X, -111.0_kPR, ierr))
349 PetscCallA(DMGetLocalSection(dmUA, sectionUA, ierr))
350 PetscCallA(DMGetLocalVector(dmUA, UALoc, ierr))
351 PetscCallA(VecGetArray(UALoc, cval, ierr))
352 PetscCallA(DMGetCoordinateSection(dmUA, coordSection, ierr))
353 PetscCallA(DMGetCoordinatesLocal(dmUA, coord, ierr))
354 PetscCallA(DMPlexGetChart(dmUA, pStart, pEnd, ierr))
357 PetscCallA(PetscSectionGetDof(sectionUA, p, dofUA, ierr))
359 PetscCallA(PetscSectionGetOffset(sectionUA, p, offUA, ierr))
360 PetscCallA(DMPlexVecGetClosure(dmUA, coordSection, coord, p, PETSC_NULL_INTEGER, xyz, ierr))
369 … PetscCallA(DMPlexVecRestoreClosure(dmUA, coordSection, coord, p, PETSC_NULL_INTEGER, xyz, ierr))
373 PetscCallA(VecRestoreArray(UALoc, cval, ierr))
374 PetscCallA(DMLocalToGlobalBegin(dmUA, UALoc, INSERT_VALUES, UA, ierr))
375 PetscCallA(DMLocalToGlobalEnd(dmUA, UALoc, INSERT_VALUES, UA, ierr))
376 PetscCallA(DMRestoreLocalVector(dmUA, UALoc, ierr))
379 PetscCallA(VecISCopy(X, isUA, SCATTER_FORWARD, UA, ierr))
381 PetscCallA(VecISCopy(X, isU, SCATTER_REVERSE, U, ierr))
382 PetscCallA(VecISCopy(X, isA, SCATTER_REVERSE, A, ierr))
383 PetscCallA(VecViewFromOptions(UA, PETSC_NULL_OBJECT, '-ua_vec_view', ierr))
384 PetscCallA(VecViewFromOptions(U, PETSC_NULL_OBJECT, '-u_vec_view', ierr))
385 PetscCallA(VecViewFromOptions(A, PETSC_NULL_OBJECT, '-a_vec_view', ierr))
387 PetscCallA(VecISCopy(X, isUA, SCATTER_REVERSE, UA2, ierr))
388 PetscCallA(VecViewFromOptions(UA2, PETSC_NULL_OBJECT, '-ua2_vec_view', ierr))
391 PetscCallA(DMSetOutputSequenceNumber(dmU, 0_kPI, time, ierr))
392 PetscCallA(DMSetOutputSequenceNumber(dmA, 0_kPI, time, ierr))
394 PetscCallA(VecView(U, viewer, ierr))
395 PetscCallA(VecView(A, viewer, ierr))
401 PetscCallA(DMSetOutputSequenceNumber(dmUA, 1_kPI, time, ierr))
402 PetscCallA(DMGetGlobalVector(dmUA, tmpVec, ierr))
403 PetscCallA(VecCopy(UA, tmpVec, ierr))
404 PetscCallA(PetscObjectSetName(tmpVec, 'U', ierr))
405 PetscCallA(VecView(tmpVec, viewer, ierr))
408 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
409 PetscCallA(VecLoad(tmpVec, viewer, ierr))
410 PetscCallA(VecAXPY(UA, -1.0_kPR, tmpVec, ierr))
411 PetscCallA(VecNorm(UA, NORM_INFINITY, norm, ierr))
415 PetscCallA(DMRestoreGlobalVector(dmUA, tmpVec, ierr))
418 PetscCallA(DMGetGlobalVector(dmUA2, tmpVec, ierr))
419 PetscCallA(VecCopy(UA2, tmpVec, ierr))
420 PetscCallA(PetscObjectSetName(tmpVec, 'U', ierr))
421 PetscCallA(DMSetOutputSequenceNumber(dmUA2, 2_kPI, time, ierr))
422 PetscCallA(VecView(tmpVec, viewer, ierr))
425 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
426 PetscCallA(VecLoad(tmpVec, viewer, ierr))
427 PetscCallA(VecAXPY(UA2, -1.0_kPR, tmpVec, ierr))
428 PetscCallA(VecNorm(UA2, NORM_INFINITY, norm, ierr))
432 PetscCallA(DMRestoreGlobalVector(dmUA2, tmpVec, ierr))
438 PetscCallA(DMGetCoordinateSection(dmS, coordSection, ierr))
439 PetscCallA(DMGetCoordinatesLocal(dmS, coord, ierr))
440 PetscCallA(DMGetLabelIdIS(dmS, 'Cell Sets', csIS, ierr))
441 PetscCallA(DMGetLabelSize(dmS, 'Cell Sets', numCS, ierr))
442 PetscCallA(ISGetIndices(csIS, csID, ierr))
445 PetscCallA(DMGetStratumIS(dmS, 'Cell Sets', csID(set), cellIS, ierr))
446 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
447 PetscCallA(ISGetSize(cellIS, numCells, ierr))
449 …allA(DMPlexVecGetClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), PETSC_NULL_INTEGER, cval, ierr))
450 …scCallA(DMPlexVecGetClosure(dmS, coordSection, coord, cellID(cell), PETSC_NULL_INTEGER, xyz, ierr))
458 …CallA(DMPlexVecSetClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), cval, INSERT_ALL_VALUES, ierr))
459 …(DMPlexVecRestoreClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), PETSC_NULL_INTEGER, cval, ierr))
460 …llA(DMPlexVecRestoreClosure(dmS, coordSection, coord, cellID(cell), PETSC_NULL_INTEGER, xyz, ierr))
462 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
463 PetscCallA(ISDestroy(cellIS, ierr))
465 PetscCallA(ISRestoreIndices(csIS, csID, ierr))
466 PetscCallA(ISDestroy(csIS, ierr))
467 PetscCallA(VecViewFromOptions(S, PETSC_NULL_OBJECT, '-s_vec_view', ierr))
470 PetscCallA(DMSetOutputSequenceNumber(dmS, 0_kPI, time, ierr))
471 PetscCallA(VecView(S, viewer, ierr))
474 PetscCallA(DMGetGlobalVector(dmS, tmpVec, ierr))
475 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
476 PetscCallA(PetscObjectSetName(tmpVec, 'Sigma', ierr))
477 PetscCallA(VecLoad(tmpVec, viewer, ierr))
478 PetscCallA(VecAXPY(S, -1.0_kPR, tmpVec, ierr))
479 PetscCallA(VecNorm(S, NORM_INFINITY, norm, ierr))
483 PetscCallA(DMRestoreGlobalVector(dmS, tmpVec, ierr))
485 PetscCallA(DMRestoreGlobalVector(dmUA2, UA2, ierr))
486 PetscCallA(DMRestoreGlobalVector(dmUA, UA, ierr))
487 PetscCallA(DMRestoreGlobalVector(dmS, S, ierr))
488 PetscCallA(DMRestoreGlobalVector(dmA, A, ierr))
489 PetscCallA(DMRestoreGlobalVector(dmU, U, ierr))
490 PetscCallA(DMRestoreGlobalVector(pdm, X, ierr))
491 PetscCallA(DMDestroy(dmU, ierr))
492 PetscCallA(ISDestroy(isU, ierr))
493 PetscCallA(DMDestroy(dmA, ierr))
494 PetscCallA(ISDestroy(isA, ierr))
495 PetscCallA(DMDestroy(dmS, ierr))
496 PetscCallA(ISDestroy(isS, ierr))
497 PetscCallA(DMDestroy(dmUA, ierr))
498 PetscCallA(ISDestroy(isUA, ierr))
499 PetscCallA(DMDestroy(dmUA2, ierr))
500 PetscCallA(DMDestroy(pdm, ierr))
502 PetscCallA(DMDestroy(dm, ierr))
508 PetscCallA(PetscViewerDestroy(viewer, ierr))
509 PetscCallA(PetscFinalize(ierr))