29 #include "Foundation/triplet.h"
30 #include "Foundation/console.h"
55 vector<pair<Vec3,pair<bool,double> > > raw_data;
56 vector<pair<Vec3,double> > data;
59 raw_data=this->m_pis->forAllInnerInteractionsGetWithPos(this->m_rdf);
62 for(vector<pair<
Vec3,pair<bool,double> > >::iterator iter=raw_data.begin();
65 if(iter->second.first){
66 data.push_back(make_pair(iter->first,iter->second.second));
71 this->m_comm->send_gather(data,0);
82 vector<pair<esys::lsm::quintuple<Vec3,double,Vec3,double,Vec3>,pair<bool,double> > > raw_data;
83 vector<pair<esys::lsm::quintuple<Vec3,double,Vec3,double,Vec3>,
double> > data;
86 raw_data=this->m_pis->forAllInnerInteractionsGetRaw2(this->m_rdf);
92 if(iter->second.first){
93 data.push_back(make_pair(iter->first,iter->second.second));
98 this->m_comm->send_gather(data,0);
106 template <
typename T>
109 vector<pair<bool,double> > data_vec;
112 this->m_pis->forAllInnerInteractionsGet(data_vec,this->m_rdf);
116 for(vector<pair<bool,double> >::iterator iter=data_vec.begin();
117 iter!=data_vec.end();
119 if(iter->first) sum+=iter->second;
122 vector<double> sum_vec;
123 sum_vec.push_back(sum);
124 this->m_comm->send_gather(sum_vec,0);
130 template <
typename T>
133 vector<pair<bool,double> > data_vec;
136 this->m_pis->forAllInnerInteractionsGet(data_vec,this->m_rdf);
141 for(vector<pair<bool,double> >::iterator iter=data_vec.begin();
142 iter!=data_vec.end();
146 max=(iter->second > max) ? iter->second : max;
154 vector<double> max_vec;
155 max_vec.push_back(max);
156 this->m_comm->send_gather(max_vec,0);
162 template <
typename T>
165 vector<pair<triplet<int,int,Vec3>, pair<bool,double> > > raw_data;
166 vector<pair<triplet<int,int,Vec3>,
double> > data;
169 raw_data=this->m_pis->forAllInnerInteractionsGetDataWithID(this->m_rdf);
171 console.
XDebug() <<
"got " << raw_data.size() <<
" raw data\n";
175 iter!=raw_data.end();
177 if(iter->second.first){
178 data.push_back(make_pair(iter->first,iter->second.second));
182 console.
XDebug() <<
"got " << data.size() <<
" filtered data\n";
185 this->m_comm->send_gather(data,0);