[624] | 1 | /*! |
---|
| 2 | \file axis_algorithm_inverse.hpp |
---|
| 3 | \author Ha NGUYEN |
---|
| 4 | \since 14 May 2015 |
---|
[630] | 5 | \date 29 June 2015 |
---|
[624] | 6 | |
---|
| 7 | \brief Algorithm for inversing an axis.. |
---|
| 8 | */ |
---|
[623] | 9 | #include "axis_algorithm_inverse.hpp" |
---|
[630] | 10 | #include "transformation_mapping.hpp" |
---|
| 11 | #include "context.hpp" |
---|
| 12 | #include "context_client.hpp" |
---|
[620] | 13 | |
---|
| 14 | namespace xios { |
---|
| 15 | |
---|
[623] | 16 | CAxisAlgorithmInverse::CAxisAlgorithmInverse(CAxis* axisDestination, CAxis* axisSource) |
---|
[622] | 17 | : CAxisAlgorithmTransformation(axisDestination, axisSource) |
---|
[620] | 18 | { |
---|
[666] | 19 | if (axisDestination->n_glo.getValue() != axisSource->n_glo.getValue()) |
---|
[620] | 20 | { |
---|
[623] | 21 | ERROR("CAxisAlgorithmInverse::CAxisAlgorithmInverse(CAxis* axisDestination, CAxis* axisSource)", |
---|
[666] | 22 | << "Two axis have different global size" |
---|
| 23 | << "Size of axis source " <<axisSource->getId() << " is " << axisSource->n_glo.getValue() << std::endl |
---|
| 24 | << "Size of axis destionation " <<axisDestination->getId() << " is " << axisDestination->n_glo.getValue()); |
---|
[620] | 25 | } |
---|
[827] | 26 | } |
---|
| 27 | |
---|
| 28 | void CAxisAlgorithmInverse::computeIndexSourceMapping_(const std::vector<CArray<double,1>* >& dataAuxInputs) |
---|
| 29 | { |
---|
| 30 | this->transformationMapping_.resize(1); |
---|
| 31 | this->transformationWeight_.resize(1); |
---|
| 32 | |
---|
[833] | 33 | TransformationIndexMap& transMap = this->transformationMapping_[0]; |
---|
| 34 | TransformationWeightMap& transWeight = this->transformationWeight_[0]; |
---|
[827] | 35 | |
---|
| 36 | int globalIndexSize = axisDestGlobalIndex_.size(); |
---|
| 37 | for (int idx = 0; idx < globalIndexSize; ++idx) |
---|
| 38 | { |
---|
| 39 | transMap[axisDestGlobalIndex_[idx]].push_back(axisDestGlobalSize_-axisDestGlobalIndex_[idx]-1); |
---|
| 40 | transWeight[axisDestGlobalIndex_[idx]].push_back(1.0); |
---|
| 41 | } |
---|
| 42 | |
---|
[666] | 43 | int niSrc = axisSrc_->n.getValue(); |
---|
| 44 | int sizeSrc = axisSrc_->n_glo.getValue(); |
---|
[630] | 45 | if (niSrc != sizeSrc) updateAxisValue(); |
---|
| 46 | else |
---|
| 47 | { |
---|
| 48 | for (int idx = 0; idx < sizeSrc; ++idx) |
---|
| 49 | { |
---|
| 50 | axisDest_->value(idx) = axisSrc_->value(sizeSrc-idx-1); |
---|
| 51 | } |
---|
| 52 | } |
---|
[620] | 53 | } |
---|
| 54 | |
---|
[630] | 55 | /*! |
---|
| 56 | Update value on axis after inversing |
---|
| 57 | After an axis is inversed, not only the data on it must be inversed but also the value |
---|
| 58 | */ |
---|
| 59 | void CAxisAlgorithmInverse::updateAxisValue() |
---|
| 60 | { |
---|
| 61 | CContext* context = CContext::getCurrent(); |
---|
| 62 | CContextClient* client=context->client; |
---|
| 63 | |
---|
[666] | 64 | int niSrc = axisSrc_->n.getValue(); |
---|
| 65 | int ibeginSrc = axisSrc_->begin.getValue(); |
---|
[630] | 66 | |
---|
| 67 | CTransformationMapping transformationMap(axisDest_, axisSrc_); |
---|
| 68 | |
---|
[833] | 69 | TransformationIndexMap& transMap = this->transformationMapping_[0]; |
---|
| 70 | TransformationWeightMap& transWeight = this->transformationWeight_[0]; |
---|
[630] | 71 | |
---|
[829] | 72 | CTransformationMapping::DestinationIndexMap globaIndexMapFromDestToSource; |
---|
[833] | 73 | TransformationIndexMap::const_iterator it = transMap.begin(), ite = transMap.end(); |
---|
[829] | 74 | int localIndex = 0; |
---|
[630] | 75 | for (; it != ite; ++it) |
---|
| 76 | { |
---|
[829] | 77 | globaIndexMapFromDestToSource[it->first].push_back(make_pair(localIndex,make_pair((it->second)[0], (transWeight[it->first])[0]))); |
---|
| 78 | ++localIndex; |
---|
[630] | 79 | } |
---|
| 80 | |
---|
| 81 | transformationMap.computeTransformationMapping(globaIndexMapFromDestToSource); |
---|
| 82 | |
---|
[829] | 83 | const CTransformationMapping::ReceivedIndexMap& globalIndexToReceive = transformationMap.getGlobalIndexReceivedOnGridDestMapping(); |
---|
| 84 | const CTransformationMapping::SentIndexMap& globalIndexToSend = transformationMap.getGlobalIndexSendToGridDestMapping(); |
---|
[630] | 85 | |
---|
| 86 | // Sending global index of original grid source |
---|
[829] | 87 | CTransformationMapping::SentIndexMap::const_iterator itbSend = globalIndexToSend.begin(), itSend, |
---|
| 88 | iteSend = globalIndexToSend.end(); |
---|
[630] | 89 | int sendBuffSize = 0; |
---|
| 90 | for (itSend = itbSend; itSend != iteSend; ++itSend) sendBuffSize += (itSend->second).size(); |
---|
| 91 | |
---|
| 92 | typedef double Scalar; |
---|
| 93 | Scalar* sendBuff, *currentSendBuff; |
---|
| 94 | if (0 != sendBuffSize) sendBuff = new Scalar[sendBuffSize]; |
---|
| 95 | for (StdSize idx = 0; idx < sendBuffSize; ++idx) sendBuff[idx] = NumTraits<Scalar>::sfmax(); |
---|
| 96 | |
---|
| 97 | int currentBuffPosition = 0; |
---|
| 98 | for (itSend = itbSend; itSend != iteSend; ++itSend) |
---|
| 99 | { |
---|
| 100 | int destRank = itSend->first; |
---|
[829] | 101 | const std::vector<std::pair<int,size_t> >& globalIndexOfCurrentGridSourceToSend = itSend->second; |
---|
[630] | 102 | int countSize = globalIndexOfCurrentGridSourceToSend.size(); |
---|
| 103 | for (int idx = 0; idx < (countSize); ++idx) |
---|
| 104 | { |
---|
[829] | 105 | int index = globalIndexOfCurrentGridSourceToSend[idx].first; |
---|
[630] | 106 | sendBuff[idx+currentBuffPosition] = (axisSrc_->value)(index); |
---|
| 107 | } |
---|
| 108 | currentSendBuff = sendBuff + currentBuffPosition; |
---|
| 109 | MPI_Send(currentSendBuff, countSize, MPI_DOUBLE, destRank, 14, client->intraComm); |
---|
| 110 | currentBuffPosition += countSize; |
---|
| 111 | } |
---|
| 112 | |
---|
| 113 | // Receiving global index of grid source sending from current grid source |
---|
[829] | 114 | CTransformationMapping::ReceivedIndexMap::const_iterator itbRecv = globalIndexToReceive.begin(), itRecv, |
---|
| 115 | iteRecv = globalIndexToReceive.end(); |
---|
[630] | 116 | int recvBuffSize = 0; |
---|
| 117 | for (itRecv = itbRecv; itRecv != iteRecv; ++itRecv) recvBuffSize += (itRecv->second).size(); |
---|
| 118 | |
---|
| 119 | Scalar* recvBuff, *currentRecvBuff; |
---|
| 120 | if (0 != recvBuffSize) recvBuff = new Scalar [recvBuffSize]; |
---|
| 121 | for (StdSize idx = 0; idx < recvBuffSize; ++idx) recvBuff[idx] = NumTraits<Scalar>::sfmax(); |
---|
| 122 | |
---|
| 123 | int currentRecvBuffPosition = 0; |
---|
| 124 | for (itRecv = itbRecv; itRecv != iteRecv; ++itRecv) |
---|
| 125 | { |
---|
| 126 | MPI_Status status; |
---|
| 127 | int srcRank = itRecv->first; |
---|
| 128 | int countSize = (itRecv->second).size(); |
---|
| 129 | currentRecvBuff = recvBuff + currentRecvBuffPosition; |
---|
| 130 | MPI_Recv(currentRecvBuff, countSize, MPI_DOUBLE, srcRank, 14, client->intraComm, &status); |
---|
| 131 | currentRecvBuffPosition += countSize; |
---|
| 132 | } |
---|
| 133 | |
---|
[666] | 134 | int ibeginDest = axisDest_->begin.getValue(); |
---|
[630] | 135 | currentRecvBuff = recvBuff; |
---|
| 136 | for (itRecv = itbRecv; itRecv != iteRecv; ++itRecv) |
---|
| 137 | { |
---|
| 138 | int countSize = (itRecv->second).size(); |
---|
| 139 | for (int idx = 0; idx < countSize; ++idx, ++currentRecvBuff) |
---|
| 140 | { |
---|
[842] | 141 | int index = ((itRecv->second)[idx]).localIndex - ibeginDest; |
---|
| 142 | (axisDest_->value)(index) = *currentRecvBuff; |
---|
| 143 | // int ssize = (itRecv->second)[idx].size(); |
---|
| 144 | // for (int i = 0; i < ssize; ++i) |
---|
| 145 | // { |
---|
| 146 | // int index = ((itRecv->second)[idx][i]).first - ibeginDest; |
---|
| 147 | // (axisDest_->value)(index) = *currentRecvBuff; |
---|
| 148 | // } |
---|
[630] | 149 | } |
---|
| 150 | } |
---|
| 151 | |
---|
| 152 | if (0 != sendBuffSize) delete [] sendBuff; |
---|
| 153 | if (0 != recvBuffSize) delete [] recvBuff; |
---|
[620] | 154 | } |
---|
[630] | 155 | |
---|
| 156 | } |
---|