[325] | 1 | /* ************************************************************************** * |
---|
[335] | 2 | * Copyright © IPSL/LSCE, xios, Avril 2010 - Octobre 2011 * |
---|
[325] | 3 | * ************************************************************************** */ |
---|
| 4 | |
---|
| 5 | #include <boost/multi_array.hpp> |
---|
| 6 | #include <boost/shared_ptr.hpp> |
---|
| 7 | #include <string> |
---|
| 8 | #include <iostream> |
---|
| 9 | |
---|
| 10 | |
---|
| 11 | #include "xmlioserver.hpp" |
---|
| 12 | #include "oasis_cinterface.hpp" |
---|
| 13 | |
---|
[352] | 14 | #include "attribute_template.hpp" |
---|
| 15 | #include "object_template.hpp" |
---|
| 16 | #include "group_template.hpp" |
---|
[325] | 17 | |
---|
| 18 | #include "icutil.hpp" |
---|
| 19 | #include "cxios.hpp" |
---|
[342] | 20 | #include "client.hpp" |
---|
[325] | 21 | #include "field.hpp" |
---|
[352] | 22 | #include "context.hpp" |
---|
[403] | 23 | #include "context_client.hpp" |
---|
[382] | 24 | #include "mpi.hpp" |
---|
[347] | 25 | #include "timer.hpp" |
---|
[369] | 26 | #include "array_new.hpp" |
---|
[325] | 27 | |
---|
[403] | 28 | |
---|
[325] | 29 | extern "C" |
---|
| 30 | { |
---|
| 31 | // /////////////////////////////// Définitions ////////////////////////////// // |
---|
| 32 | |
---|
| 33 | // ----------------------- Redéfinition de types ---------------------------- |
---|
| 34 | |
---|
| 35 | typedef enum { NETCDF4 = 0 } XFileType; |
---|
| 36 | |
---|
[345] | 37 | typedef xios::CContext * XContextPtr; |
---|
[325] | 38 | |
---|
| 39 | // -------------------- Traitement des données ------------------------------ |
---|
| 40 | void cxios_init_server(void) |
---|
| 41 | { |
---|
| 42 | CXios::initServerSide(); |
---|
| 43 | } |
---|
| 44 | |
---|
| 45 | void cxios_init_client(const char * client_id , int len_client_id, MPI_Fint* f_local_comm, MPI_Fint* f_return_comm ) |
---|
| 46 | { |
---|
| 47 | std::string str; |
---|
| 48 | MPI_Comm local_comm ; |
---|
| 49 | MPI_Comm return_comm ; |
---|
| 50 | |
---|
| 51 | if (!cstr2string(client_id, len_client_id, str)) return; |
---|
[347] | 52 | |
---|
[325] | 53 | int initialized ; |
---|
| 54 | MPI_Initialized(&initialized) ; |
---|
| 55 | if (initialized) local_comm=MPI_Comm_f2c(*f_local_comm) ; |
---|
| 56 | else local_comm=MPI_COMM_NULL ; |
---|
| 57 | CXios::initClientSide(str,local_comm,return_comm); |
---|
| 58 | *f_return_comm=MPI_Comm_c2f(return_comm) ; |
---|
[361] | 59 | CTimer::get("XIOS init").suspend() ; |
---|
| 60 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 61 | } |
---|
| 62 | |
---|
| 63 | void cxios_context_initialize(const char * context_id , int len_context_id, MPI_Fint* f_comm) |
---|
| 64 | { |
---|
| 65 | std::string str; |
---|
| 66 | MPI_Comm comm ; |
---|
| 67 | |
---|
| 68 | if (!cstr2string(context_id, len_context_id, str)) return; |
---|
[347] | 69 | CTimer::get("XIOS").resume() ; |
---|
| 70 | CTimer::get("XIOS init context").resume() ; |
---|
[325] | 71 | comm=MPI_Comm_f2c(*f_comm) ; |
---|
[342] | 72 | CClient::registerContext(str,comm) ; |
---|
[347] | 73 | CTimer::get("XIOS init context").suspend() ; |
---|
| 74 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 75 | } |
---|
| 76 | |
---|
| 77 | void cxios_context_close_definition() |
---|
| 78 | { |
---|
[347] | 79 | CTimer::get("XIOS").resume() ; |
---|
| 80 | CTimer::get("XIOS close definition").resume() ; |
---|
| 81 | CContext* context = CContext::getCurrent() ; |
---|
[325] | 82 | context->closeDefinition() ; |
---|
[347] | 83 | CTimer::get("XIOS close definition").suspend() ; |
---|
| 84 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 85 | } |
---|
| 86 | |
---|
| 87 | void cxios_context_finalize() |
---|
| 88 | { |
---|
[347] | 89 | CTimer::get("XIOS").resume() ; |
---|
| 90 | CTimer::get("XIOS context finalize").resume() ; |
---|
| 91 | CContext* context = CContext::getCurrent() ; |
---|
[325] | 92 | context->finalize() ; |
---|
[347] | 93 | CTimer::get("XIOS context finalize").suspend() ; |
---|
| 94 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 95 | } |
---|
| 96 | |
---|
| 97 | void cxios_finalize() |
---|
| 98 | { |
---|
[347] | 99 | CTimer::get("XIOS").resume() ; |
---|
| 100 | CTimer::get("XIOS finalize").resume() ; |
---|
[325] | 101 | CXios::clientFinalize() ; |
---|
| 102 | } |
---|
| 103 | |
---|
[445] | 104 | void cxios_solve_inheritance() |
---|
| 105 | { |
---|
| 106 | CTimer::get("XIOS").resume() ; |
---|
| 107 | CContext* context = CContext::getCurrent() ; |
---|
| 108 | context->solveAllInheritance(false) ; |
---|
| 109 | CTimer::get("XIOS").suspend() ; |
---|
| 110 | } |
---|
[325] | 111 | |
---|
| 112 | // ---------------------- Ecriture des données ------------------------------ |
---|
| 113 | |
---|
| 114 | void cxios_write_data_k81(const char * fieldid, int fieldid_size, double * data_k8, int data_Xsize) |
---|
| 115 | { |
---|
| 116 | std::string fieldid_str; |
---|
| 117 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
| 118 | |
---|
[347] | 119 | CTimer::get("XIOS").resume() ; |
---|
| 120 | CTimer::get("XIOS send field").resume() ; |
---|
| 121 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 122 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
[369] | 123 | CArray<double,(StdSize)1> data(data_k8,shape(data_Xsize),neverDeleteData) ; |
---|
[325] | 124 | CField::get(fieldid_str)->setData(data) ; |
---|
[369] | 125 | CField toto ; |
---|
| 126 | toto.setData(data) ; |
---|
[347] | 127 | CTimer::get("XIOS send field").suspend() ; |
---|
| 128 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 129 | } |
---|
| 130 | |
---|
| 131 | void cxios_write_data_k82(const char * fieldid, int fieldid_size, double * data_k8, int data_Xsize, int data_Ysize) |
---|
| 132 | { |
---|
| 133 | std::string fieldid_str; |
---|
| 134 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
| 135 | |
---|
[347] | 136 | CTimer::get("XIOS").resume() ; |
---|
| 137 | CTimer::get("XIOS send field").resume() ; |
---|
| 138 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 139 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
[347] | 140 | |
---|
[369] | 141 | CArray<double,2>data(data_k8,shape(data_Xsize,data_Ysize),neverDeleteData) ; |
---|
[325] | 142 | CField::get(fieldid_str)->setData(data) ; |
---|
[347] | 143 | CTimer::get("XIOS send field").suspend() ; |
---|
| 144 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 145 | } |
---|
| 146 | |
---|
| 147 | void cxios_write_data_k83(const char * fieldid, int fieldid_size, double * data_k8, int data_Xsize, int data_Ysize, int data_Zsize) |
---|
| 148 | { |
---|
| 149 | std::string fieldid_str; |
---|
| 150 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
[347] | 151 | |
---|
| 152 | CTimer::get("XIOS").resume() ; |
---|
| 153 | CTimer::get("XIOS send field").resume() ; |
---|
| 154 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 155 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
[369] | 156 | |
---|
| 157 | CArray<double,3>data(data_k8,shape(data_Xsize,data_Ysize,data_Zsize),neverDeleteData) ; |
---|
[325] | 158 | CField::get(fieldid_str)->setData(data) ; |
---|
[347] | 159 | CTimer::get("XIOS send field").suspend() ; |
---|
| 160 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 161 | } |
---|
| 162 | |
---|
| 163 | void cxios_write_data_k41(const char * fieldid, int fieldid_size, float * data_k4, int data_Xsize) |
---|
| 164 | { |
---|
| 165 | std::string fieldid_str; |
---|
[347] | 166 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
| 167 | |
---|
| 168 | CTimer::get("XIOS").resume() ; |
---|
| 169 | CTimer::get("XIOS send field").resume() ; |
---|
| 170 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 171 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
| 172 | |
---|
[369] | 173 | CArray<float,1> data_tmp(data_k4,shape(data_Xsize),neverDeleteData) ; |
---|
| 174 | CArray<double,1> data(data_Xsize) ; |
---|
| 175 | data=data_tmp ; |
---|
[325] | 176 | CField::get(fieldid_str)->setData(data) ; |
---|
[347] | 177 | CTimer::get("XIOS send field").suspend() ; |
---|
| 178 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 179 | } |
---|
| 180 | |
---|
| 181 | void cxios_write_data_k42(const char * fieldid, int fieldid_size, float * data_k4, int data_Xsize, int data_Ysize) |
---|
| 182 | { |
---|
| 183 | std::string fieldid_str; |
---|
| 184 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
[347] | 185 | |
---|
| 186 | CTimer::get("XIOS").resume() ; |
---|
| 187 | CTimer::get("XIOS send field").resume() ; |
---|
| 188 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 189 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
| 190 | |
---|
[369] | 191 | CArray<float,2> data_tmp(data_k4,shape(data_Xsize,data_Ysize),neverDeleteData) ; |
---|
| 192 | CArray<double,2> data(data_Xsize,data_Ysize) ; |
---|
| 193 | data=data_tmp ; |
---|
[325] | 194 | CField::get(fieldid_str)->setData(data) ; |
---|
[347] | 195 | CTimer::get("XIOS send field").suspend() ; |
---|
| 196 | CTimer::get("XIOS").suspend() ; |
---|
[325] | 197 | } |
---|
| 198 | |
---|
| 199 | void cxios_write_data_k43(const char * fieldid, int fieldid_size, float * data_k4, int data_Xsize, int data_Ysize, int data_Zsize) |
---|
| 200 | { |
---|
| 201 | std::string fieldid_str; |
---|
[347] | 202 | |
---|
[325] | 203 | if (!cstr2string(fieldid, fieldid_size, fieldid_str)) return; |
---|
[347] | 204 | |
---|
| 205 | CTimer::get("XIOS").resume() ; |
---|
| 206 | CTimer::get("XIOS send field").resume() ; |
---|
| 207 | CContext* context = CContext::getCurrent() ; |
---|
[403] | 208 | if (!context->hasServer) context->client->checkBuffers() ; |
---|
[369] | 209 | |
---|
| 210 | CArray<float,3> data_tmp(data_k4,shape(data_Xsize,data_Ysize,data_Zsize),neverDeleteData) ; |
---|
| 211 | CArray<double,3> data(data_Xsize,data_Ysize,data_Zsize) ; |
---|
| 212 | data=data_tmp ; |
---|
| 213 | |
---|
[325] | 214 | CField::get(fieldid_str)->setData(data) ; |
---|
[347] | 215 | CTimer::get("XIOS send field").suspend() ; |
---|
| 216 | CTimer::get("XIOS").suspend() ; |
---|
| 217 | |
---|
[325] | 218 | } |
---|
| 219 | |
---|
| 220 | } // extern "C" |
---|