Lines Matching refs:ierr
31 PetscErrorCode :: ierr
81 PetscCallA(PetscInitialize(PETSC_NULL_CHARACTER, ierr))
82 if (ierr /= 0) then
87 PetscCallA(MPI_Comm_rank(PETSC_COMM_WORLD, rank, ierr))
88 PetscCallA(MPI_Comm_size(PETSC_COMM_WORLD, numProc, ierr))
89 …cCallA(PetscOptionsGetString(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-i', ifilename, flg, ierr))
91 …cCallA(PetscOptionsGetString(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-o', ofilename, flg, ierr))
93 …etscCallA(PetscOptionsGetInt(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-order', order, flg, ierr))
100 …allA(DMPlexCreateFromFile(PETSC_COMM_WORLD, ifilename, PETSC_NULL_CHARACTER, PETSC_TRUE, dm, ierr))
101 PetscCallA(DMPlexDistributeSetDefault(dm, PETSC_FALSE, ierr))
102 PetscCallA(DMSetFromOptions(dm, ierr))
103 PetscCallA(DMGetDimension(dm, sdim, ierr))
104 PetscCallA(DMViewFromOptions(dm, PETSC_NULL_OBJECT, '-dm_view', ierr))
109 PetscCallA(exopts(EXVRBS + EXDEBG, ierr))
111 PetscCallA(PetscViewerExodusIIOpen(PETSC_COMM_WORLD, ofilename, FILE_MODE_WRITE, viewer, ierr))
120 PetscCallA(PetscViewerExodusIISetOrder(viewer, order, ierr))
121 PetscCallA(PetscViewerView(viewer, PETSC_VIEWER_STDOUT_WORLD, ierr))
129 PetscCallA(DMView(dm, viewer, ierr))
130 PetscCallA(PetscViewerView(viewer, PETSC_VIEWER_STDOUT_WORLD, ierr))
149 PetscCallA(PetscViewerExodusIIGetId(viewer, exoid, ierr))
150 PetscCallA(expvp(exoid, 'E', numZonalVar, ierr))
151 PetscCallA(expvan(exoid, 'E', numZonalVar, zonalVarName, ierr))
152 PetscCallA(expvp(exoid, 'N', numNodalVar, ierr))
153 PetscCallA(expvan(exoid, 'N', numNodalVar, nodalVarName, ierr))
154 PetscCallA(exinq(exoid, EX_INQ_ELEM_BLK, numCS, PETSC_NULL_REAL, sjunk, ierr))
160 PetscCallA(expvtt(exoid, numCS, numZonalVar, truthtable, ierr))
165 PetscCallA(exptim(exoid, step, real(step, kind=kPR), ierr))
168 PetscCallA(DMSetUseNatural(dm, PETSC_TRUE, ierr))
169 PetscCallA(DMPlexGetPartitioner(dm, part, ierr))
170 PetscCallA(PetscPartitionerSetFromOptions(part, ierr))
171 PetscCallA(DMPlexDistribute(dm, 0_kPI, migrationSF, pdm, ierr))
174 PetscCallA(DMPlexSetMigrationSF(pdm, migrationSF, ierr))
175 PetscCallA(PetscSFDestroy(migrationSF, ierr))
179 PetscCallA(DMViewFromOptions(pdm, PETSC_NULL_OBJECT, '-dm_view', ierr))
181 PetscCallA(PetscObjectGetComm(pdm, comm, ierr))
182 PetscCallA(PetscSectionCreate(comm, section, ierr))
183 PetscCallA(PetscSectionSetNumFields(section, 3_kPI, ierr))
184 PetscCallA(PetscSectionSetFieldName(section, fieldU, 'U', ierr))
185 PetscCallA(PetscSectionSetFieldName(section, fieldA, 'Alpha', ierr))
186 PetscCallA(PetscSectionSetFieldName(section, fieldS, 'Sigma', ierr))
187 PetscCallA(DMPlexGetChart(pdm, pStart, pEnd, ierr))
188 PetscCallA(PetscSectionSetChart(section, pStart, pEnd, ierr))
193 PetscCallA(DMPlexGetDepthStratum(pdm, d - 1, pStartDepth(d), pEndDepth(d), ierr))
197 PetscCallA(PetscSectionSetFieldComponents(section, fieldU, sdim, ierr))
198 PetscCallA(PetscSectionSetFieldComponents(section, fieldA, 1_kPI, ierr))
199 PetscCallA(PetscSectionSetFieldComponents(section, fieldS, sdim*(sdim + 1)/2, ierr))
202 PetscCallA(DMGetLabelSize(pdm, 'Cell Sets', numCS, ierr))
203 PetscCallA(DMGetLabelIdIS(pdm, 'Cell Sets', csIS, ierr))
204 PetscCallA(ISGetIndices(csIS, csID, ierr))
206 PetscCallA(DMGetStratumSize(pdm, 'Cell Sets', csID(set), numCells, ierr))
207 PetscCallA(DMGetStratumIS(pdm, 'Cell Sets', csID(set), cellIS, ierr))
221 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
223 …scCallA(DMPlexGetTransitiveClosure(pdm, cellID(1), PETSC_TRUE, PETSC_NULL_INTEGER, closureA, ierr))
261 …llA(DMPlexRestoreTransitiveClosure(pdm, cellID(1), PETSC_TRUE, PETSC_NULL_INTEGER, closureA, ierr))
264 …CallA(DMPlexGetTransitiveClosure(pdm, cellID(cell), PETSC_TRUE, PETSC_NULL_INTEGER, closure, ierr))
269 PetscCallA(PetscSectionSetDof(section, closure(p), dofU(d) + dofA(d) + dofS(d), ierr))
270 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldU, dofU(d), ierr))
271 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldA, dofA(d), ierr))
272 PetscCallA(PetscSectionSetFieldDof(section, closure(p), fieldS, dofS(d), ierr))
276 …A(DMPlexRestoreTransitiveClosure(pdm, cellID(cell), PETSC_TRUE, PETSC_NULL_INTEGER, closure, ierr))
278 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
279 PetscCallA(ISDestroy(cellIS, ierr))
282 PetscCallA(ISRestoreIndices(csIS, csID, ierr))
283 PetscCallA(ISDestroy(csIS, ierr))
284 PetscCallA(PetscSectionSetUp(section, ierr))
285 PetscCallA(DMSetLocalSection(pdm, section, ierr))
286 PetscCallA(PetscObjectViewFromOptions(section, PETSC_NULL_OBJECT, '-dm_section_view', ierr))
287 PetscCallA(PetscSectionDestroy(section, ierr))
291 PetscCallA(DMPlexGetMigrationSF(pdm, natPointSF, ierr))
292 PetscCallA(DMGetLocalSection(pdm, rootSection, ierr))
293 PetscCallA(PetscSFCreateInverseSF(natPointSF, natPointSFInv, ierr))
294 PetscCallA(PetscSectionCreate(PETSC_COMM_WORLD, leafSection, ierr))
295 … PetscCallA(PetscSFDistributeSection(natPointSFInv, rootSection, remoteOffsets, leafSection, ierr))
296 PetscCallA(PetscSFDestroyRemoteOffsets(remoteOffsets, ierr))
297 PetscCallA(DMSetLocalSection(dm, leafSection, ierr))
298 PetscCallA(DMPlexCreateGlobalToNaturalSF(pdm, leafSection, natPointSF, natSF, ierr))
299 PetscCallA(PetscSFDestroy(natPointSFInv, ierr))
300 PetscCallA(PetscSectionDestroy(leafSection, ierr))
301 PetscCallA(DMSetNaturalSF(pdm, natSF, ierr))
302 PetscCallA(PetscObjectDereference(natSF, ierr))
306 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldU], isU, dmU, ierr))
307 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldA], isA, dmA, ierr))
308 PetscCallA(DMCreateSubDM(pdm, 1_kPI, [fieldS], isS, dmS, ierr))
309 PetscCallA(DMCreateSubDM(pdm, 2_kPI, fieldUA, isUA, dmUA, ierr))
315 PetscCallA(DMCreateSuperDM(dmList, 2_kPI, PETSC_NULL_IS_POINTER, dmUA2, ierr))
318 PetscCallA(DMGetGlobalVector(pdm, X, ierr))
319 PetscCallA(DMGetGlobalVector(dmU, U, ierr))
320 PetscCallA(DMGetGlobalVector(dmA, A, ierr))
321 PetscCallA(DMGetGlobalVector(dmS, S, ierr))
322 PetscCallA(DMGetGlobalVector(dmUA, UA, ierr))
323 PetscCallA(DMGetGlobalVector(dmUA2, UA2, ierr))
325 PetscCallA(PetscObjectSetName(U, 'U', ierr))
326 PetscCallA(PetscObjectSetName(A, 'Alpha', ierr))
327 PetscCallA(PetscObjectSetName(S, 'Sigma', ierr))
328 PetscCallA(PetscObjectSetName(UA, 'UAlpha', ierr))
329 PetscCallA(PetscObjectSetName(UA2, 'UAlpha2', ierr))
330 PetscCallA(VecSet(X, -111.0_kPR, ierr))
333 PetscCallA(DMGetLocalSection(dmUA, sectionUA, ierr))
334 PetscCallA(DMGetLocalVector(dmUA, UALoc, ierr))
335 PetscCallA(VecGetArray(UALoc, cval, ierr))
336 PetscCallA(DMGetCoordinateSection(dmUA, coordSection, ierr))
337 PetscCallA(DMGetCoordinatesLocal(dmUA, coord, ierr))
338 PetscCallA(DMPlexGetChart(dmUA, pStart, pEnd, ierr))
341 PetscCallA(PetscSectionGetDof(sectionUA, p, dofUA, ierr))
343 PetscCallA(PetscSectionGetOffset(sectionUA, p, offUA, ierr))
344 PetscCallA(DMPlexVecGetClosure(dmUA, coordSection, coord, p, PETSC_NULL_INTEGER, xyz, ierr))
353 … PetscCallA(DMPlexVecRestoreClosure(dmUA, coordSection, coord, p, PETSC_NULL_INTEGER, xyz, ierr))
357 PetscCallA(VecRestoreArray(UALoc, cval, ierr))
358 PetscCallA(DMLocalToGlobalBegin(dmUA, UALoc, INSERT_VALUES, UA, ierr))
359 PetscCallA(DMLocalToGlobalEnd(dmUA, UALoc, INSERT_VALUES, UA, ierr))
360 PetscCallA(DMRestoreLocalVector(dmUA, UALoc, ierr))
363 PetscCallA(VecISCopy(X, isUA, SCATTER_FORWARD, UA, ierr))
365 PetscCallA(VecISCopy(X, isU, SCATTER_REVERSE, U, ierr))
366 PetscCallA(VecISCopy(X, isA, SCATTER_REVERSE, A, ierr))
367 PetscCallA(VecViewFromOptions(UA, PETSC_NULL_OBJECT, '-ua_vec_view', ierr))
368 PetscCallA(VecViewFromOptions(U, PETSC_NULL_OBJECT, '-u_vec_view', ierr))
369 PetscCallA(VecViewFromOptions(A, PETSC_NULL_OBJECT, '-a_vec_view', ierr))
371 PetscCallA(VecISCopy(X, isUA, SCATTER_REVERSE, UA2, ierr))
372 PetscCallA(VecViewFromOptions(UA2, PETSC_NULL_OBJECT, '-ua2_vec_view', ierr))
375 PetscCallA(DMSetOutputSequenceNumber(dmU, 0_kPI, time, ierr))
376 PetscCallA(DMSetOutputSequenceNumber(dmA, 0_kPI, time, ierr))
378 PetscCallA(VecView(U, viewer, ierr))
379 PetscCallA(VecView(A, viewer, ierr))
385 PetscCallA(DMSetOutputSequenceNumber(dmUA, 1_kPI, time, ierr))
386 PetscCallA(DMGetGlobalVector(dmUA, tmpVec, ierr))
387 PetscCallA(VecCopy(UA, tmpVec, ierr))
388 PetscCallA(PetscObjectSetName(tmpVec, 'U', ierr))
389 PetscCallA(VecView(tmpVec, viewer, ierr))
392 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
393 PetscCallA(VecLoad(tmpVec, viewer, ierr))
394 PetscCallA(VecAXPY(UA, -1.0_kPR, tmpVec, ierr))
395 PetscCallA(VecNorm(UA, NORM_INFINITY, norm, ierr))
399 PetscCallA(DMRestoreGlobalVector(dmUA, tmpVec, ierr))
402 PetscCallA(DMGetGlobalVector(dmUA2, tmpVec, ierr))
403 PetscCallA(VecCopy(UA2, tmpVec, ierr))
404 PetscCallA(PetscObjectSetName(tmpVec, 'U', ierr))
405 PetscCallA(DMSetOutputSequenceNumber(dmUA2, 2_kPI, time, ierr))
406 PetscCallA(VecView(tmpVec, viewer, ierr))
409 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
410 PetscCallA(VecLoad(tmpVec, viewer, ierr))
411 PetscCallA(VecAXPY(UA2, -1.0_kPR, tmpVec, ierr))
412 PetscCallA(VecNorm(UA2, NORM_INFINITY, norm, ierr))
416 PetscCallA(DMRestoreGlobalVector(dmUA2, tmpVec, ierr))
422 PetscCallA(DMGetCoordinateSection(dmS, coordSection, ierr))
423 PetscCallA(DMGetCoordinatesLocal(dmS, coord, ierr))
424 PetscCallA(DMGetLabelIdIS(dmS, 'Cell Sets', csIS, ierr))
425 PetscCallA(DMGetLabelSize(dmS, 'Cell Sets', numCS, ierr))
426 PetscCallA(ISGetIndices(csIS, csID, ierr))
429 PetscCallA(DMGetStratumIS(dmS, 'Cell Sets', csID(set), cellIS, ierr))
430 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
431 PetscCallA(ISGetSize(cellIS, numCells, ierr))
433 …allA(DMPlexVecGetClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), PETSC_NULL_INTEGER, cval, ierr))
434 …scCallA(DMPlexVecGetClosure(dmS, coordSection, coord, cellID(cell), PETSC_NULL_INTEGER, xyz, ierr))
442 …CallA(DMPlexVecSetClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), cval, INSERT_ALL_VALUES, ierr))
443 …(DMPlexVecRestoreClosure(dmS, PETSC_NULL_SECTION, S, cellID(cell), PETSC_NULL_INTEGER, cval, ierr))
444 …llA(DMPlexVecRestoreClosure(dmS, coordSection, coord, cellID(cell), PETSC_NULL_INTEGER, xyz, ierr))
446 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
447 PetscCallA(ISDestroy(cellIS, ierr))
449 PetscCallA(ISRestoreIndices(csIS, csID, ierr))
450 PetscCallA(ISDestroy(csIS, ierr))
451 PetscCallA(VecViewFromOptions(S, PETSC_NULL_OBJECT, '-s_vec_view', ierr))
454 PetscCallA(DMSetOutputSequenceNumber(dmS, 0_kPI, time, ierr))
455 PetscCallA(VecView(S, viewer, ierr))
458 PetscCallA(DMGetGlobalVector(dmS, tmpVec, ierr))
459 PetscCallA(VecSet(tmpVec, -1000.0_kPR, ierr))
460 PetscCallA(PetscObjectSetName(tmpVec, 'Sigma', ierr))
461 PetscCallA(VecLoad(tmpVec, viewer, ierr))
462 PetscCallA(VecAXPY(S, -1.0_kPR, tmpVec, ierr))
463 PetscCallA(VecNorm(S, NORM_INFINITY, norm, ierr))
467 PetscCallA(DMRestoreGlobalVector(dmS, tmpVec, ierr))
469 PetscCallA(DMRestoreGlobalVector(dmUA2, UA2, ierr))
470 PetscCallA(DMRestoreGlobalVector(dmUA, UA, ierr))
471 PetscCallA(DMRestoreGlobalVector(dmS, S, ierr))
472 PetscCallA(DMRestoreGlobalVector(dmA, A, ierr))
473 PetscCallA(DMRestoreGlobalVector(dmU, U, ierr))
474 PetscCallA(DMRestoreGlobalVector(pdm, X, ierr))
475 PetscCallA(DMDestroy(dmU, ierr))
476 PetscCallA(ISDestroy(isU, ierr))
477 PetscCallA(DMDestroy(dmA, ierr))
478 PetscCallA(ISDestroy(isA, ierr))
479 PetscCallA(DMDestroy(dmS, ierr))
480 PetscCallA(ISDestroy(isS, ierr))
481 PetscCallA(DMDestroy(dmUA, ierr))
482 PetscCallA(ISDestroy(isUA, ierr))
483 PetscCallA(DMDestroy(dmUA2, ierr))
484 PetscCallA(DMDestroy(pdm, ierr))
486 PetscCallA(DMDestroy(dm, ierr))
492 PetscCallA(PetscViewerDestroy(viewer, ierr))
493 PetscCallA(PetscFinalize(ierr))