int
UmfpackGenLinSOE::sendSelf(int cTag, Channel &theChannel)
{
LinearSOESolver *theSoeSolver = this->getSolver();
if (theSoeSolver != 0) {
if (theSoeSolver->sendSelf(cTag, theChannel) < 0) {
opserr <<"WARNING MumpsParallelSOE::sendSelf() - failed to send solver\n";
return -1;
}
} else {
opserr <<"WARNING MumpsParallelSOE::sendSelf() - no solver to send!\n";
return -1;
}
return 0;
}
int
MumpsParallelSOE::sendSelf(int commitTag, Channel &theChannel)
{
int sendID =0;
// if P0 check if already sent. If already sent use old processID; if not allocate a new process
// id for remote part of object, enlarge channel * to hold a channel * for this remote object.
// if not P0, send current processID
if (processID == 0) {
// check if already using this object
bool found = false;
for (int i=0; i<numChannels; i++)
if (theChannels[i] == &theChannel) {
sendID = i+1;
found = true;
}
// if new object, enlarge Channel pointers to hold new channel * & allocate new ID
if (found == false) {
int nextNumChannels = numChannels + 1;
Channel **nextChannels = new Channel *[nextNumChannels];
if (nextNumChannels == 0) {
opserr << "MumpsParallelSOE::sendSelf() - failed to allocate channel array of size: " <<
nextNumChannels << endln;
return -1;
}
for (int i=0; i<numChannels; i++)
nextChannels[i] = theChannels[i];
nextChannels[numChannels] = &theChannel;
numChannels = nextNumChannels;
if (theChannels != 0)
delete [] theChannels;
theChannels = nextChannels;
if (localCol != 0)
delete [] localCol;
localCol = new ID *[numChannels];
if (localCol == 0) {
opserr << "MumpsParallelSOE::sendSelf() - failed to allocate id array of size: " <<
nextNumChannels << endln;
return -1;
}
for (int i=0; i<numChannels; i++)
localCol[i] = 0;
// allocate new processID for remote object
sendID = numChannels;
}
} else
sendID = processID;
// send remotes processID
ID idData(2);
idData(0) = sendID;
idData(1) = matType;
int res = theChannel.sendID(0, commitTag, idData);
if (res < 0) {
opserr <<"WARNING MumpsParallelSOE::sendSelf() - failed to send data\n";
return -1;
}
LinearSOESolver *theSoeSolver = this->getSolver();
if (theSoeSolver != 0) {
if (theSoeSolver->sendSelf(commitTag, theChannel) < 0) {
opserr <<"WARNING MumpsParallelSOE::sendSelf() - failed to send solver\n";
return -1;
}
} else {
opserr <<"WARNING MumpsParallelSOE::sendSelf() - no solver to send!\n";
return -1;
}
return 0;
}
