#include "xios_spl.hpp" #include "cxios.hpp" #include "client.hpp" #include "server.hpp" #include "xml_parser.hpp" #include #include "mpi.hpp" #include "memory.hpp" #include #include "memtrack.hpp" #include "registry.hpp" using namespace ep_lib; namespace xios { const string CXios::rootFile="./iodef.xml" ; const string CXios::xiosCodeId="xios.x" ; const string CXios::clientFile="./xios_client"; const string CXios::serverFile="./xios_server"; const string CXios::serverPrmFile="./xios_server1"; const string CXios::serverSndFile="./xios_server2"; bool CXios::isClient ; bool CXios::isServer ; MPI_Comm CXios::globalComm ; bool CXios::usingOasis ; bool CXios::usingServer = false; bool CXios::usingServer2 = false; int CXios::ratioServer2 = 50; int CXios::nbPoolsServer2 = 1; double CXios::bufferSizeFactor = 1.0; const double CXios::defaultBufferSizeFactor = 1.0; StdSize CXios::minBufferSize = 1024 * sizeof(double); StdSize CXios::maxBufferSize = std::numeric_limits::max() ; bool CXios::printLogs2Files; bool CXios::isOptPerformance = true; CRegistry* CXios::globalRegistry = 0; double CXios::recvFieldTimeout = 300.0; bool CXios::checkEventSync=false ; //! Parse configuration file and create some objects from it void CXios::initialize() { set_new_handler(noMemory); int tmp_rank; MPI_Comm_rank(MPI_COMM_WORLD, &tmp_rank); #pragma omp critical { std::cout<<"thread "<("using_oasis",false) ; usingServer=getin("using_server",false) ; usingServer2=getin("using_server2",false) ; ratioServer2=getin("ratio_server2",50); nbPoolsServer2=getin("number_pools_server2",1); info.setLevel(getin("info_level",0)) ; report.setLevel(getin("info_level",50)); printLogs2Files=getin("print_file",false); StdString bufMemory("memory"); StdString bufPerformance("performance"); StdString bufOpt = getin("optimal_buffer_size", bufPerformance); std::transform(bufOpt.begin(), bufOpt.end(), bufOpt.begin(), ::tolower); if (0 == bufOpt.compare(bufMemory)) isOptPerformance = false; else if (0 != bufOpt.compare(bufPerformance)) { ERROR("CXios::parseXiosConfig()", << "optimal_buffer_size must be memory or performance "<< endl ); } bufferSizeFactor = getin("buffer_size_factor", defaultBufferSizeFactor); minBufferSize = getin("min_buffer_size", 1024 * sizeof(double)); maxBufferSize = getin("max_buffer_size", std::numeric_limits::max()); recvFieldTimeout = getin("recv_field_timeout", recvFieldTimeout); if (recvFieldTimeout < 0.0) ERROR("CXios::parseXiosConfig()", "recv_field_timeout cannot be negative."); checkEventSync = getin("check_event_sync", checkEventSync); //globalComm=MPI_COMM_WORLD ; int num_ep; if(isClient) { num_ep = omp_get_num_threads(); } if(isServer) { num_ep = 1; } MPI_Info info; #pragma omp master { MPI_Comm *ep_comm; MPI_Comm_create_endpoints(MPI_COMM_WORLD->mpi_comm, num_ep, info, ep_comm); // servers should reach here too. passage = ep_comm; //::MPI_Comm_group(to_mpi_comm(MPI_COMM_WORLD->mpi_comm), &MPI_GROUP_WORLD); } #pragma omp barrier CXios::globalComm = passage[omp_get_thread_num()]; } /*! Initialize client \param [in] codeId identity of context \param [in] localComm local communicator \param [in/out] returnComm communicator corresponding to group of client with same codeId */ void CXios::initClientSide(const string& codeId, MPI_Comm& localComm, MPI_Comm& returnComm) { isClient = true; initialize() ; CClient::initialize(codeId,localComm,returnComm) ; if (CClient::getRank()==0) globalRegistry = new CRegistry(returnComm) ; // If there are no server processes then we are in attached mode // and the clients are also servers isServer = !usingServer; if (printLogs2Files) { #pragma omp critical CClient::openInfoStream(clientFile); CClient::openErrorStream(clientFile); } else { CClient::openInfoStream(); CClient::openErrorStream(); } } void CXios::clientFinalize(void) { CClient::finalize() ; if (CClient::getRank()==0) { #pragma omp critical (_output) info(80)<<"Write data base Registry"<toString()<toFile("xios_registry.bin") ; delete globalRegistry ; } #ifdef XIOS_MEMTRACK #ifdef XIOS_MEMTRACK_LIGHT report(10) << " Memory report : current memory used by XIOS : "<< MemTrack::getCurrentMemorySize()*1.0/(1024*1024)<<" Mbyte" << endl ; report(10) << " Memory report : maximum memory used by XIOS : "<< MemTrack::getMaxMemorySize()*1.0/(1024*1024)<<" Mbyte" << endl ; #endif #ifdef XIOS_MEMTRACK_FULL MemTrack::TrackListMemoryUsage() ; MemTrack::TrackDumpBlocks(); #endif CClient::closeInfoStream(); #endif } //! Init server by parsing only xios part of config file void CXios::initServer() { set_new_handler(noMemory); std::set parseList; parseList.insert("xios"); xml::CXMLParser::ParseFile(rootFile, parseList); parseXiosConfig(); } //! Initialize server then put it into listening state void CXios::initServerSide(void) { isClient = false; isServer = true; initServer(); // Initialize all aspects MPI CServer::initialize(); if (CServer::getRank()==0 && CServer::serverLevel != 1) globalRegistry = new CRegistry(CServer::intraComm) ; if (printLogs2Files) { if (CServer::serverLevel == 0) { CServer::openInfoStream(serverFile); CServer::openErrorStream(serverFile); } else if (CServer::serverLevel == 1) { CServer::openInfoStream(serverPrmFile); CServer::openErrorStream(serverPrmFile); } else { CServer::openInfoStream(serverSndFile); CServer::openErrorStream(serverSndFile); } } else { CServer::openInfoStream(); CServer::openErrorStream(); } // Enter the loop to listen message from Client CServer::eventLoop(); // Finalize if (CServer::serverLevel == 0) { if (CServer::getRank()==0) { info(80)<<"Write data base Registry"<toString()<toFile("xios_registry.bin") ; delete globalRegistry ; } } else { // If using two server levels: // (1) merge registries on each pool // (2) send merged registries to the first pool // (3) merge received registries on the first pool if (CServer::serverLevel == 2) { vector& secondaryServerGlobalRanks = CServer::getSecondaryServerGlobalRanks(); int firstPoolGlobalRank = secondaryServerGlobalRanks[0]; int rankGlobal; MPI_Comm_rank(globalComm, &rankGlobal); // Merge registries defined on each pools CRegistry globalRegistrySndServers (CServer::intraComm); // All pools (except the first): send globalRegistry to the first pool for (int i=1; i