OpenFPM_pdata  1.1.0
Project that contain the implementation of distributed structures
 All Data Structures Namespaces Functions Variables Typedefs Enumerations Friends Pages
vector_dist_operator_assign.hpp
1 /*
2  * vector_dist_operator_assign.ipp
3  *
4  * Created on: Jul 21, 2016
5  * Author: i-bird
6  */
7 
8 #ifndef OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATOR_ASSIGN_HPP_
9 #define OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATOR_ASSIGN_HPP_
10 
13 template<typename T>
14 struct construct_expression
15 {
17  static inline const T & construct(const T & e)
18  {
19  return e;
20  }
21 };
22 
24 template<>
25 struct construct_expression<double>
26 {
28  static inline vector_dist_expression<0,double> construct(double e)
29  {
30  return vector_dist_expression<0,double>(e);
31  }
32 };
33 
35 template<>
36 struct construct_expression<float>
37 {
39  static inline vector_dist_expression<0,float> construct(const float e)
40  {
41  return vector_dist_expression<0,float>(e);
42  }
43 };
44 
45 template<typename prp1, typename expr1, typename prp2, typename expr2> void assign(prp1 & p1, const expr1 & v_e1, prp2 & p2, const expr2 & v_e2)
46 {
47  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
48  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
49 
50  v_exp1.init();
51  v_exp2.init();
52 
53  auto it = p1.getVector().getDomainIterator();
54 
55  while (it.isNext())
56  {
57  auto key = it.get();
58 
59  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
60  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
61 
62  ++it;
63  }
64 }
65 
66 template<typename prp1, typename expr1, typename prp2, typename expr2, typename prp3, typename expr3> void assign(prp1 & p1, const expr1 & v_e1, prp2 & p2, const expr2 & v_e2, prp3 & p3, const expr3 & v_e3)
67 {
68  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
69  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
70  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
71 
72  v_exp1.init();
73  v_exp2.init();
74  v_exp3.init();
75 
76  auto it = p1.getVector().getDomainIterator();
77 
78  while (it.isNext())
79  {
80  auto key = it.get();
81 
82  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
83  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
84  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
85 
86  ++it;
87  }
88 }
89 
90 template<typename prp1, typename expr1,
91  typename prp2, typename expr2,
92  typename prp3, typename expr3,
93  typename prp4, typename expr4>
94 void assign(prp1 & p1, const expr1 & v_e1,
95  prp2 & p2, const expr2 & v_e2,
96  prp3 & p3, const expr3 & v_e3,
97  prp4 & p4, const expr4 & v_e4)
98 {
99  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
100  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
101  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
102  auto v_exp4 = construct_expression<expr4>::construct(v_e4);
103 
104  v_exp1.init();
105  v_exp2.init();
106  v_exp3.init();
107  v_exp4.init();
108 
109  auto it = p1.getVector().getDomainIterator();
110 
111  while (it.isNext())
112  {
113  auto key = it.get();
114 
115  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
116  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
117  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
118  pos_or_propL<typename prp4::vtype,prp4::prop>::value(p4.getVector(),key) = v_exp4.value(key);
119 
120  ++it;
121  }
122 }
123 
124 
125 template<typename prp1, typename expr1,
126  typename prp2, typename expr2,
127  typename prp3, typename expr3,
128  typename prp4, typename expr4,
129  typename prp5, typename expr5>
130 void assign(prp1 & p1, const expr1 & v_e1,
131  prp2 & p2, const expr2 & v_e2,
132  prp3 & p3, const expr3 & v_e3,
133  prp4 & p4, const expr4 & v_e4,
134  prp5 & p5, const expr5 & v_e5)
135 {
136  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
137  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
138  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
139  auto v_exp4 = construct_expression<expr4>::construct(v_e4);
140  auto v_exp5 = construct_expression<expr5>::construct(v_e5);
141 
142  v_exp1.init();
143  v_exp2.init();
144  v_exp3.init();
145  v_exp4.init();
146  v_exp5.init();
147 
148  auto it = p1.getVector().getDomainIterator();
149 
150  while (it.isNext())
151  {
152  auto key = it.get();
153 
154  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
155  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
156  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
157  pos_or_propL<typename prp4::vtype,prp4::prop>::value(p4.getVector(),key) = v_exp4.value(key);
158  pos_or_propL<typename prp5::vtype,prp5::prop>::value(p5.getVector(),key) = v_exp5.value(key);
159 
160  ++it;
161  }
162 }
163 
164 template<typename prp1, typename expr1,
165  typename prp2, typename expr2,
166  typename prp3, typename expr3,
167  typename prp4, typename expr4,
168  typename prp5, typename expr5,
169  typename prp6, typename expr6>
170 void assign(prp1 & p1, const expr1 & v_e1,
171  prp2 & p2, const expr2 & v_e2,
172  prp3 & p3, const expr3 & v_e3,
173  prp4 & p4, const expr4 & v_e4,
174  prp5 & p5, const expr5 & v_e5,
175  prp6 & p6, const expr6 & v_e6)
176 {
177  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
178  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
179  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
180  auto v_exp4 = construct_expression<expr4>::construct(v_e4);
181  auto v_exp5 = construct_expression<expr5>::construct(v_e5);
182  auto v_exp6 = construct_expression<expr6>::construct(v_e6);
183 
184  v_exp1.init();
185  v_exp2.init();
186  v_exp3.init();
187  v_exp4.init();
188  v_exp5.init();
189  v_exp6.init();
190 
191  auto it = p1.getVector().getDomainIterator();
192 
193  while (it.isNext())
194  {
195  auto key = it.get();
196 
197  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
198  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
199  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
200  pos_or_propL<typename prp4::vtype,prp4::prop>::value(p4.getVector(),key) = v_exp4.value(key);
201  pos_or_propL<typename prp5::vtype,prp5::prop>::value(p5.getVector(),key) = v_exp5.value(key);
202  pos_or_propL<typename prp6::vtype,prp6::prop>::value(p6.getVector(),key) = v_exp6.value(key);
203 
204  ++it;
205  }
206 }
207 
208 
209 
210 template<typename prp1, typename expr1,
211  typename prp2, typename expr2,
212  typename prp3, typename expr3,
213  typename prp4, typename expr4,
214  typename prp5, typename expr5,
215  typename prp6, typename expr6,
216  typename prp7, typename expr7>
217 void assign(prp1 & p1, const expr1 & v_e1,
218  prp2 & p2, const expr2 & v_e2,
219  prp3 & p3, const expr3 & v_e3,
220  prp4 & p4, const expr4 & v_e4,
221  prp5 & p5, const expr5 & v_e5,
222  prp6 & p6, const expr6 & v_e6,
223  prp7 & p7, const expr7 & v_e7)
224 {
225  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
226  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
227  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
228  auto v_exp4 = construct_expression<expr4>::construct(v_e4);
229  auto v_exp5 = construct_expression<expr5>::construct(v_e5);
230  auto v_exp6 = construct_expression<expr6>::construct(v_e6);
231  auto v_exp7 = construct_expression<expr7>::construct(v_e7);
232 
233  v_exp1.init();
234  v_exp2.init();
235  v_exp3.init();
236  v_exp4.init();
237  v_exp5.init();
238  v_exp6.init();
239  v_exp7.init();
240 
241  auto it = p1.getVector().getDomainIterator();
242 
243  while (it.isNext())
244  {
245  auto key = it.get();
246 
247  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
248  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
249  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
250  pos_or_propL<typename prp4::vtype,prp4::prop>::value(p4.getVector(),key) = v_exp4.value(key);
251  pos_or_propL<typename prp5::vtype,prp5::prop>::value(p5.getVector(),key) = v_exp5.value(key);
252  pos_or_propL<typename prp6::vtype,prp6::prop>::value(p6.getVector(),key) = v_exp6.value(key);
253  pos_or_propL<typename prp7::vtype,prp7::prop>::value(p7.getVector(),key) = v_exp7.value(key);
254 
255  ++it;
256  }
257 }
258 
259 
260 template<typename prp1, typename expr1,
261  typename prp2, typename expr2,
262  typename prp3, typename expr3,
263  typename prp4, typename expr4,
264  typename prp5, typename expr5,
265  typename prp6, typename expr6,
266  typename prp7, typename expr7,
267  typename prp8, typename expr8>
268 void assign(prp1 & p1, const expr1 & v_e1,
269  prp2 & p2, const expr2 & v_e2,
270  prp3 & p3, const expr3 & v_e3,
271  prp4 & p4, const expr4 & v_e4,
272  prp5 & p5, const expr5 & v_e5,
273  prp6 & p6, const expr6 & v_e6,
274  prp7 & p7, const expr7 & v_e7,
275  prp8 & p8, const expr8 & v_e8)
276 {
277  auto v_exp1 = construct_expression<expr1>::construct(v_e1);
278  auto v_exp2 = construct_expression<expr2>::construct(v_e2);
279  auto v_exp3 = construct_expression<expr3>::construct(v_e3);
280  auto v_exp4 = construct_expression<expr4>::construct(v_e4);
281  auto v_exp5 = construct_expression<expr5>::construct(v_e5);
282  auto v_exp6 = construct_expression<expr6>::construct(v_e6);
283  auto v_exp7 = construct_expression<expr7>::construct(v_e7);
284  auto v_exp8 = construct_expression<expr8>::construct(v_e8);
285 
286  v_exp1.init();
287  v_exp2.init();
288  v_exp3.init();
289  v_exp4.init();
290  v_exp5.init();
291  v_exp6.init();
292  v_exp7.init();
293  v_exp8.init();
294 
295  auto it = p1.getVector().getDomainIterator();
296 
297  while (it.isNext())
298  {
299  auto key = it.get();
300 
301  pos_or_propL<typename prp1::vtype,prp1::prop>::value(p1.getVector(),key) = v_exp1.value(key);
302  pos_or_propL<typename prp2::vtype,prp2::prop>::value(p2.getVector(),key) = v_exp2.value(key);
303  pos_or_propL<typename prp3::vtype,prp3::prop>::value(p3.getVector(),key) = v_exp3.value(key);
304  pos_or_propL<typename prp4::vtype,prp4::prop>::value(p4.getVector(),key) = v_exp4.value(key);
305  pos_or_propL<typename prp5::vtype,prp5::prop>::value(p5.getVector(),key) = v_exp5.value(key);
306  pos_or_propL<typename prp6::vtype,prp6::prop>::value(p6.getVector(),key) = v_exp6.value(key);
307  pos_or_propL<typename prp7::vtype,prp7::prop>::value(p7.getVector(),key) = v_exp7.value(key);
308  pos_or_propL<typename prp8::vtype,prp8::prop>::value(p8.getVector(),key) = v_exp8.value(key);
309 
310  ++it;
311  }
312 }
313 
314 #endif /* OPENFPM_NUMERICS_SRC_OPERATORS_VECTOR_VECTOR_DIST_OPERATOR_ASSIGN_HPP_ */
pos_or_propL::value
static auto value(vector &v, const vect_dist_key_dx &k) -> decltype(v.template getProp< prp >(k))
return the value (position or property) of the particle k in the vector v
Definition: vector_dist_operators.hpp:327
construct_expression
Definition: vector_dist_operator_assign.hpp:14
vector_dist_expression
Main class that encapsulate a vector properties operand to be used for expressions construction...
Definition: vector_dist_operators.hpp:433
construct_expression< double >::construct
static vector_dist_expression< 0, double > construct(double e)
It return an expression from a double.
Definition: vector_dist_operator_assign.hpp:28
key
This class is a trick to indicate the compiler a specific specialization pattern. ...
Definition: memory_c.hpp:201
construct_expression::construct
static const T & construct(const T &e)
It return the expression itself.
Definition: vector_dist_operator_assign.hpp:17
construct_expression< float >::construct
static vector_dist_expression< 0, float > construct(const float e)
It return an expression from a float.
Definition: vector_dist_operator_assign.hpp:39