2 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
4 #include "bbstdVectorFilterDouble.h"
5 #include "bbstdPackage.h"
9 BBTK_ADD_BLACK_BOX_TO_PACKAGE(std,VectorFilterDouble)
10 BBTK_BLACK_BOX_IMPLEMENTATION(VectorFilterDouble,bbtk::AtomicBlackBox);
12 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
14 void VectorFilterDouble::Process()
17 // THE MAIN PROCESSING METHOD BODY
18 // Here we simply set the input 'In' value to the output 'Out'
19 // And print out the output value
20 // INPUT/OUTPUT ACCESSORS ARE OF THE FORM :
21 // void bbSet{Input|Output}NAME(const TYPE&)
22 // const TYPE& bbGet{Input|Output}NAME() const
24 // * NAME is the name of the input/output
25 // (the one provided in the attribute 'name' of the tag 'input')
26 // * TYPE is the C++ type of the input/output
27 // (the one provided in the attribute 'type' of the tag 'input')
29 std::vector< std::vector<double> * > pLstVec;
30 std::vector< std::vector<double> * > pLstVecOut;
31 std::vector<double> In0 = bbGetInputIn0();
32 std::vector<double> In1 = bbGetInputIn1();
33 std::vector<double> In2 = bbGetInputIn2();
34 std::vector<double> In3 = bbGetInputIn3();
35 std::vector<double> In4 = bbGetInputIn4();
36 std::vector<double> In5 = bbGetInputIn5();
37 std::vector<double> In6 = bbGetInputIn6();
38 std::vector<double> In7 = bbGetInputIn7();
39 std::vector<double> In8 = bbGetInputIn8();
40 std::vector<double> In9 = bbGetInputIn9();
42 if (bbGetInputIn0().size()!=0) { pLstVec.push_back( &In0 ); }
43 if (bbGetInputIn1().size()!=0) { pLstVec.push_back( &In1 ); }
44 if (bbGetInputIn2().size()!=0) { pLstVec.push_back( &In2 ); }
45 if (bbGetInputIn3().size()!=0) { pLstVec.push_back( &In3 ); }
46 if (bbGetInputIn4().size()!=0) { pLstVec.push_back( &In4 ); }
47 if (bbGetInputIn5().size()!=0) { pLstVec.push_back( &In5 ); }
48 if (bbGetInputIn6().size()!=0) { pLstVec.push_back( &In6 ); }
49 if (bbGetInputIn7().size()!=0) { pLstVec.push_back( &In7 ); }
50 if (bbGetInputIn8().size()!=0) { pLstVec.push_back( &In8 ); }
51 if (bbGetInputIn9().size()!=0) { pLstVec.push_back( &In9 ); }
53 std::vector<double> Out0;
54 std::vector<double> Out1;
55 std::vector<double> Out2;
56 std::vector<double> Out3;
57 std::vector<double> Out4;
58 std::vector<double> Out5;
59 std::vector<double> Out6;
60 std::vector<double> Out7;
61 std::vector<double> Out8;
62 std::vector<double> Out9;
64 pLstVecOut.push_back( &Out0 );
65 pLstVecOut.push_back( &Out1 );
66 pLstVecOut.push_back( &Out2 );
67 pLstVecOut.push_back( &Out3 );
68 pLstVecOut.push_back( &Out4 );
69 pLstVecOut.push_back( &Out5 );
70 pLstVecOut.push_back( &Out6 );
71 pLstVecOut.push_back( &Out7 );
72 pLstVecOut.push_back( &Out8 );
73 pLstVecOut.push_back( &Out9 );
76 if (bbGetInputType()==0) // Erase duplicate lines
80 for (ipLstvec=1;ipLstvec<pLstVec.size();ipLstvec++)
82 if ( (*pLstVec[ipLstvec]).size()!=(*pLstVec[0]).size() ) { okSizeVec=false; }
84 if ( okSizeVec==false)
86 printf("EED VectorFilterDouble::Process WARNING! vectors are not of the same size.\n");
91 for (iLine=0 ; iLine < (*pLstVec[0]).size() ; iLine++ )
94 for ( ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
96 if ( (*pLstVec[ipLstvec2])[iLine]!=(*pLstVec[ipLstvec2])[iLine-1] ) { okLine=true; }
99 if ( (okLine==true) || (iLine==0) )
101 for (ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
103 (*pLstVecOut[ipLstvec2]).push_back( (*pLstVec[ipLstvec2])[iLine] );
111 if (bbGetInputType()==1) // Resize vector
115 int sizeLines = (*pLstVec[0]).size();
116 double step = (double)sizeLines/((double)bbGetInputk1()-1);
117 for (iLine=0 ; iLine<sizeLines ; iLine=iLine+step)
119 for (ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
121 (*pLstVecOut[ipLstvec2]).push_back( (*pLstVec[ipLstvec2])[(int)iLine] );
125 // adding the lastone
126 for (ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
128 (*pLstVecOut[ipLstvec2]).push_back( (*pLstVec[ipLstvec2])[(int)(sizeLines-1)] );
133 if (bbGetInputType()==2) // Insert intermediate points
139 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
141 size=(*pLstVec[ipLstvec]).size();
144 for (i=0;i<size-1;i++)
146 (*pLstVecOut[ipLstvec]).push_back( (*pLstVec[ipLstvec])[i] ); // First item
147 result=((*pLstVec[ipLstvec])[i] + (*pLstVec[ipLstvec])[i+1] ) / 2;
148 (*pLstVecOut[ipLstvec]).push_back( result );
150 (*pLstVecOut[ipLstvec]).push_back( (*pLstVec[ipLstvec])[size-1] ); // Last item
156 if (bbGetInputType()==3) // Addition k1
162 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
164 size=(*pLstVec[ipLstvec]).size();
167 result= (*pLstVec[ipLstvec])[i] + bbGetInputk1() ;
168 (*pLstVecOut[ipLstvec]).push_back( result );
173 if (bbGetInputType()==4) // Substraction k1
179 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
181 size=(*pLstVec[ipLstvec]).size();
184 result= (*pLstVec[ipLstvec])[i] - bbGetInputk1() ;
185 (*pLstVecOut[ipLstvec]).push_back( result );
190 if (bbGetInputType()==5) // Multiplication k1
196 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
198 size=(*pLstVec[ipLstvec]).size();
201 result= (*pLstVec[ipLstvec])[i] * bbGetInputk1() ;
202 (*pLstVecOut[ipLstvec]).push_back( result );
207 if (bbGetInputType()==6) // Division k1ng gadget (node check, owner GtkCheckButton)
215 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
217 size=(*pLstVec[ipLstvec]).size();
220 if (bbGetInputk1()!=0)
222 result= (*pLstVec[ipLstvec])[i] / bbGetInputk1() ;
223 (*pLstVecOut[ipLstvec]).push_back( result );
225 printf("EED VectorFilterDouble::Process Warning!! Div by 0 in Type 6 \n");
226 (*pLstVecOut[ipLstvec]).push_back( -1 );
234 if (bbGetInputType()==7)
237 for (j=0;j<In4.size(); j++)
239 Out0.push_back( In4[j] );
240 Out1.push_back( In5[j] );
241 Out2.push_back( In6[j] );
245 int iGeneral=0,sizeI,sizeLst1Indexes=In3.size();
246 for (i=0; i<sizeLst1Indexes-1;i++) { iGeneral=iGeneral+In3[i]; }
247 sizeI=iGeneral+In3[ sizeLst1Indexes-1 ];
250 int sizeJ = In7[ 0 ];
252 double dist2,dist2Min,dx,dy,dz;
255 for (j=jGeneral;j<sizeJ;j++)
258 for (i=iGeneral;i<sizeI;i++)
263 dist2 = dx*dx + dy*dy + dz*dz;
280 if (bbGetInputType()==8) // order all vectors using as base the vector In0
282 printf("EED VectorFilterDouble::Process Type 8 Start\n");
285 sizeII=pLstVec.size();
289 for (i=0;i<size; i++)
291 for (ii=0;ii<sizeII;ii++)
293 (*pLstVecOut[ii]).push_back( (*pLstVec[ii])[i] );
300 if ((*pLstVecOut[0])[j] <= (*pLstVecOut[0])[i] )
303 for (ii=0;ii<sizeII;ii++)
305 tmp = (*pLstVecOut[ii])[i];
306 (*pLstVecOut[ii])[i] = (*pLstVecOut[ii])[j];
307 (*pLstVecOut[ii])[j] = tmp;
308 // tmp=Out0[i]; Out0[i]=Out0[j]; Out0[j]=tmp;
313 printf("EED VectorFilterDouble::Process Type 8 End\n");
320 bbSetOutputOut0( Out0 );
321 bbSetOutputOut1( Out1 );
322 bbSetOutputOut2( Out2 );
323 bbSetOutputOut3( Out3 );
324 bbSetOutputOut4( Out4 );
325 bbSetOutputOut5( Out5 );
326 bbSetOutputOut6( Out6 );
327 bbSetOutputOut7( Out7 );
328 bbSetOutputOut8( Out8 );
329 bbSetOutputOut9( Out9 );
332 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
334 void VectorFilterDouble::bbUserSetDefaultValues()
337 // SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX
338 // Here we initialize the input 'In' to 0
344 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
346 void VectorFilterDouble::bbUserInitializeProcessing()
349 // THE INITIALIZATION METHOD BODY :
351 // but this is where you should allocate the internal/output pointers
357 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
359 void VectorFilterDouble::bbUserFinalizeProcessing()
362 // THE FINALIZATION METHOD BODY :
364 // but this is where you should desallocate the internal/output pointers
369 // EO namespace bbstd