ESyS-Particle  4.0.1
VectorParticleFieldSlave.hpp
00001 
00002 //                                                         //
00003 // Copyright (c) 2003-2011 by The University of Queensland //
00004 // Earth Systems Science Computational Centre (ESSCC)      //
00005 // http://www.uq.edu.au/esscc                              //
00006 //                                                         //
00007 // Primary Business: Brisbane, Queensland, Australia       //
00008 // Licensed under the Open Software License version 3.0    //
00009 // http://www.opensource.org/licenses/osl-3.0.php          //
00010 //                                                         //
00012 
00013 //-- STL includes --
00014 #include <vector>
00015 #include <utility>
00016 
00017 using std::vector;
00018 using std::pair;
00019 
00020 // -- IO includes --
00021 #include <iostream>
00022 
00023 using std::cout;
00024 using std::endl;
00025 
00033 template <typename T>
00034 VectorParticleFieldSlave<T>::VectorParticleFieldSlave(TML_Comm* comm,ParallelParticleArray<T>* ppa,typename T::VectorFieldFunction rdf):AFieldSlave(comm)
00035 {
00036   m_ppa=ppa;
00037   m_rdf=rdf;
00038 }
00039 
00043 template <typename T>
00044 void VectorParticleFieldSlave<T>::sendData()
00045 { 
00046   vector<pair<int,Vec3> > data_vec;
00047   vector<pair<int,Vec3> > pos_vec;
00048 
00049   data_vec=m_ppa->forAllInnerParticlesGetIndexed(m_rdf);
00050   pos_vec=m_ppa->forAllInnerParticlesGetIndexed(typename T::VectorFieldFunction(&T::getPos));
00051  
00052   // send data to master
00053   m_comm->send_gather(data_vec,0);
00054   m_comm->send_gather(pos_vec,0);
00055 }