19 #ifdef HAVE_MUELU_EPETRAEXT
29 #include <arpa/inet.h>
41 const std::string& levelName = inListEntry->first;
44 if (inList.
isSublist(levelName) && ((levelName.find(
"level ") == 0 && levelName.size() > 6) || levelName.find(
"user data") == 0)) {
45 int levelID = strtol(levelName.substr(6).c_str(), 0, 0);
47 if (levelName.find(
"user data") == std::string::npos) {
49 levelID = strtol(levelName.substr(6).c_str(), 0, 0);
50 if (maxLevel < levelID)
57 const std::string& name = levelListEntry->first;
59 bool isNumberedNullspace = (name.rfind(
"Nullspace", 0) == 0 && name.length() > 9 && std::isdigit(name.back(), std::locale::classic()));
61 if (name ==
"A" || name ==
"P" || name ==
"R" || name ==
"M" || name ==
"Mdiag" || name ==
"K" || name ==
"Nullspace" || name ==
"Material" || name ==
"Coordinates" || name ==
"D0" || name ==
"Dk_1" || name ==
"Dk_2" || name ==
"Mk_one" || name ==
"Mk_1_one" || name ==
"M1_beta" || name ==
"M1_alpha" || name ==
"invMk_1_invBeta" || name ==
"invMk_2_invAlpha" || name ==
"M1" || name ==
"Ms" || name ==
"M0inv" || name ==
"Pnodal" || name ==
"NodeMatrix" || name ==
"NodeAggMatrix" || name ==
"Node Comm" || name ==
"DualNodeID2PrimalNodeID"
62 #ifdef HAVE_MUELU_INTREPID2
63 || name ==
"pcoarsen: element to node map"
65 || name ==
"output stream" || isNumberedNullspace) {
66 nonSerialList.
sublist(levelName).
setEntry(name, levelListEntry->second);
68 #ifdef HAVE_MUELU_MATLAB
70 nonSerialList.
sublist(levelName).
setEntry(name, levelListEntry->second);
74 nonSerialList.
sublist(levelName).
setEntry(name, levelListEntry->second);
81 serialList.
setEntry(inListEntry->first, inListEntry->second);
91 char* buf = (
char*)malloc(stream.size() + 1);
92 strcpy(buf, stream.c_str());
93 char* token = strtok(buf, delimChars);
101 char* tokStart = token;
102 char* tokEnd = token + strlen(token) - 1;
103 while (*tokStart ==
' ' && tokStart < tokEnd)
105 while (*tokEnd ==
' ' && tokStart < tokEnd)
108 if (tokStart < tokEnd) {
109 std::string finishedToken(tokStart, tokEnd - tokStart);
110 tokenList.push_back(finishedToken);
112 token = strtok(NULL, delimChars);
119 char* str = (
char*)malloc(name.length() + 1);
120 strcpy(str, name.c_str());
122 char* firstWord = strtok(str,
" ");
127 char* secondWord = strtok(NULL,
" ");
132 char* thirdWord = strtok(NULL,
" ");
138 char* tolowerIt = firstWord;
140 *tolowerIt = (char)
tolower(*tolowerIt);
144 if (strstr(firstWord,
"matrix") ||
145 strstr(firstWord,
"multivector") ||
146 strstr(firstWord,
"map") ||
147 strstr(firstWord,
"ordinalvector") ||
148 strstr(firstWord,
"int") ||
149 strstr(firstWord,
"scalar") ||
150 strstr(firstWord,
"double") ||
151 strstr(firstWord,
"complex") ||
152 strstr(firstWord,
"string"))
165 char* str = (
char*)malloc(name.length() + 1);
166 strcpy(str, name.c_str());
168 char* firstWord = strtok(str,
" ");
173 char* secondWord = strtok(NULL,
" ");
178 char* thirdWord = strtok(NULL,
" ");
184 char* tolowerIt = firstWord;
186 *tolowerIt = (char)
tolower(*tolowerIt);
190 if (strstr(firstWord,
"matrix") ||
191 strstr(firstWord,
"multivector") ||
192 strstr(firstWord,
"map") ||
193 strstr(firstWord,
"ordinalvector") ||
194 strstr(firstWord,
"int") ||
195 strstr(firstWord,
"scalar") ||
196 strstr(firstWord,
"double") ||
197 strstr(firstWord,
"complex") ||
198 strstr(firstWord,
"string") ||
199 strstr(firstWord,
"array<go>") ||
200 strstr(firstWord,
"array<lo>") ||
201 strstr(firstWord,
"arrayrcp<lo>") ||
202 strstr(firstWord,
"arrayrcp<go>"))
216 int numRanks = baseComm->getSize();
218 NodeId = baseComm->getRank();
223 char hostname[MPI_MAX_PROCESSOR_NAME];
225 MPI_Get_processor_name(hostname, &len);
226 struct hostent* host = gethostbyname(hostname);
227 int myaddr = (int)inet_addr(inet_ntoa(*(
struct in_addr*)host->h_addr));
230 std::vector<int> addressList(numRanks);
231 Teuchos::gatherAll(*baseComm, 1, &myaddr, numRanks, &addressList[0]);
234 std::sort(addressList.begin(), addressList.end());
238 for (
int i = 0, prev = addressList[0]; i < numRanks && prev != myaddr; i++) {
239 if (prev != addressList[i]) {
240 prev = addressList[i];
250 if (reductionFactor != 1) {
253 int coresPerNode = 0;
254 for (
int i = 0, prev = addressList[0]; i < numRanks; i++) {
255 if (prev != addressList[i]) {
256 prev = addressList[i];
257 coresPerNode = std::max(i - lastI, coresPerNode);
261 coresPerNode = std::max(numRanks - lastI, coresPerNode);
264 if (coresPerNode % reductionFactor != 0)
265 throw std::runtime_error(
"Reduction factor does not evently divide # cores per node");
266 int reducedCPN = coresPerNode / reductionFactor;
267 int nodeDivision = newComm->getRank() / reducedCPN;
269 NodeId = numNodes * reductionFactor + nodeDivision;
270 newComm = baseComm->split(NodeId, baseComm->getRank());
275 NodeId = baseComm->getRank();
ParameterList & setEntry(const std::string &name, U &&entry)
ConstIterator end() const
std::string tolower(const std::string &str)
bool IsParamMuemexVariable(const std::string &name)
bool isSublist(const std::string &name) const
params_t::ConstIterator ConstIterator
ConstIterator begin() const
bool IsParamValidVariable(const std::string &name)
ParameterList & sublist(const std::string &name, bool mustAlreadyExist=false, const std::string &docString="")
void TokenizeStringAndStripWhiteSpace(const std::string &stream, std::vector< std::string > &tokenList, const char *delimChars)
Teuchos::RCP< const Teuchos::Comm< int > > GenerateNodeComm(RCP< const Teuchos::Comm< int > > &baseComm, int &NodeId, const int reductionFactor)
long ExtractNonSerializableData(const Teuchos::ParameterList &inList, Teuchos::ParameterList &serialList, Teuchos::ParameterList &nonSerialList)
Extract non-serializable data from level-specific sublists and move it to a separate parameter list...