static char help[] =
"...\n\n";
++cOunter;
return 0;
}
private:
int &cOunter;
};
--cOunter;
return 0;
}
private:
int &cOunter;
};
int main(
int argc,
char *argv[]) {
try {
PetscBool flg = PETSC_TRUE;
#if PETSC_VERSION_GE(3, 6, 4)
255, &flg);
#else
#endif
if (flg != PETSC_TRUE) {
SETERRQ(PETSC_COMM_SELF, 1, "*** ERROR -my_file (MESH FILE NEEDED)");
}
DMType dm_name = "DMMOFEM";
DM dm;
CHKERR DMCreate(PETSC_COMM_WORLD, &dm);
CHKERR DMSetType(dm, dm_name);
ParallelComm *pcomm = ParallelComm::get_pcomm(&moab,
MYPCOMM_INDEX);
auto moab_comm_wrap =
boost::make_shared<WrapMPIComm>(PETSC_COMM_WORLD, false);
if (pcomm == NULL)
pcomm = new ParallelComm(&moab, moab_comm_wrap->get_comm());
const char *option;
option = "PARALLEL=BCAST_DELETE;PARALLEL_RESOLVE_SHARED_ENTS;PARTITION="
"PARALLEL_PARTITION;";
bit_level0.set(0);
root_set, 3, bit_level0);
int field_rank = 3;
#if PETSC_VERSION_GE(3, 6, 4)
&flg);
#else
&field_rank, &flg);
#endif
"FIELD1");
"FIELD2");
"FE_ONLY_DATA");
dm, PETSC_FALSE);
->checkMPIAIJWithArraysMatrixFillIn<PetscGlobalIdx_mi_tag>(
"TEST_PROBLEM", -1, -1, 1);
->checkMPIAIJMatrixFillIn<PetscGlobalIdx_mi_tag>("TEST_PROBLEM", -1, -1,
1);
std::vector<std::string> fields_list;
fields_list.push_back("FIELD1");
PetscSection section;
§ion);
CHKERR PetscSectionView(section, PETSC_VIEWER_STDOUT_WORLD);
CHKERR DMSetSection(dm, section);
CHKERR PetscSectionDestroy(§ion);
PetscBool save_file = PETSC_TRUE;
#if PETSC_VERSION_GE(3, 6, 4)
&flg);
#else
&save_file, &flg);
#endif
if (save_file) {
PetscViewer viewer;
CHKERR PetscViewerASCIIOpen(PETSC_COMM_WORLD,
"dm_build_partitioned_mesh.txt", &viewer);
CHKERR PetscViewerPushFormat(viewer, PETSC_VIEWER_ASCII_INFO);
CHKERR PetscViewerDestroy(&viewer);
}
int count = 0;
boost::make_shared<CountUp>(count));
boost::make_shared<CountDown>(count));
if (count)
count);
auto check_consistency_of_uids = [&]() {
for (auto &fe : *field_ents_ptr) {
if (fe->getOwnerProc() !=
<< "UId owner proc "
<< "Entity owner proc " << fe->getOwnerProc();
"UId and entity handle inconsistency");
}
if (fe->getEnt() !=
<< "UId handle "
MOFEM_LOG(
"SELF", Sev::error) <<
"Entity handle " << fe->getEnt();
"UId and entity handle inconsistency");
}
const auto field_bit_number_from_uid =
if (fe->getBitNumber() != field_bit_number_from_uid) {
<< "UId field bit " << field_bit_number_from_uid;
<< "Entity owner proc " << fe->getBitNumber();
"UId and entity handle inconsistency");
}
field_bit_number_from_uid - 1))) {
<< "UId field bit "
<< "Entity owner proc " << fe->getName();
"UId and entity handle inconsistency");
}
}
};
CHKERR check_consistency_of_uids();
}
return 0;
}