]> Creatis software - bbtk.git/blob - packages/std/src/bbstdVectorFilterDouble.cxx
#3068 BBTK Feature New Normal - bbstdVectorFilterDouble addition and sustraction
[bbtk.git] / packages / std / src / bbstdVectorFilterDouble.cxx
1 //===== 
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)
3 //===== 
4 #include "bbstdVectorFilterDouble.h"
5 #include "bbstdPackage.h"
6 namespace bbstd
7 {
8
9 BBTK_ADD_BLACK_BOX_TO_PACKAGE(std,VectorFilterDouble)
10 BBTK_BLACK_BOX_IMPLEMENTATION(VectorFilterDouble,bbtk::AtomicBlackBox);
11 //===== 
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)
13 //===== 
14 void VectorFilterDouble::Process()
15 {
16
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 
23 //    Where :
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')
28
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();
41
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 ); }
52
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;
63
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 );
74
75
76         if (bbGetInputType()==0)  // Erase duplicate lines
77         {
78                 bool okSizeVec=true;
79                 int ipLstvec; 
80                 for (ipLstvec=1;ipLstvec<pLstVec.size();ipLstvec++)
81                 {
82                         if ( (*pLstVec[ipLstvec]).size()!=(*pLstVec[0]).size() )  { okSizeVec=false; }
83                 }
84                 if ( okSizeVec==false)
85                 {
86                         printf("EED VectorFilterDouble::Process WARNING! vectors are not of the same size.\n");
87                 } else {
88                         int iLine;
89                         int ipLstvec2; 
90                         bool okLine;
91                         for (iLine=0 ; iLine < (*pLstVec[0]).size() ; iLine++ ) 
92                         {
93                                 okLine=false;
94                                 for ( ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
95                                 {
96                                         if ( (*pLstVec[ipLstvec2])[iLine]!=(*pLstVec[ipLstvec2])[iLine-1] )  { okLine=true; }
97                                 } // for ipLstVec2
98                                                 
99                                 if ( (okLine==true) || (iLine==0) ) 
100                                 {
101                                         for (ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
102                                         {
103                                                 (*pLstVecOut[ipLstvec2]).push_back( (*pLstVec[ipLstvec2])[iLine] );
104                                         } // for
105                                 } // if okLine
106                         } // for iLine
107                 } // if okSizeVec
108         } // Type==0
109
110
111         if (bbGetInputType()==1) // Resize vector
112         {
113                 int             ipLstvec2;
114                 double  iLine;
115                 int     sizeLines       = (*pLstVec[0]).size();
116                 double  step            = (double)sizeLines/(double)bbGetInputk1();
117                 for (iLine=0 ; iLine<sizeLines ; iLine=iLine+step)
118                 {
119                         for (ipLstvec2=0 ; ipLstvec2<pLstVec.size() ; ipLstvec2++)
120                         {
121                                 (*pLstVecOut[ipLstvec2]).push_back( (*pLstVec[ipLstvec2])[(int)iLine] );
122                         } // for                        
123                 } // for ipLstVec
124         } // Type==1
125
126         if (bbGetInputType()==2) // Insert intermediate points
127         {
128                 int             ipLstvec;
129                 int     i;
130                 int     size;
131                 double  result;
132                 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
133                 {
134                         size=(*pLstVec[ipLstvec]).size();
135                         if (size>=2)
136                         {
137                                 for (i=0;i<size-1;i++)
138                                 {
139                                         (*pLstVecOut[ipLstvec]).push_back( (*pLstVec[ipLstvec])[i] );   // First item                   
140                                         result=((*pLstVec[ipLstvec])[i] + (*pLstVec[ipLstvec])[i+1]  ) / 2;
141                                         (*pLstVecOut[ipLstvec]).push_back( result );                            
142                                 }// for size
143                                 (*pLstVecOut[ipLstvec]).push_back( (*pLstVec[ipLstvec])[size-1] );      // Last item            
144                         } // if size>=2
145                 } // for pLstVec                
146         } // Type==2
147
148
149         if (bbGetInputType()==3) // Addition  k1
150         {
151                 int             ipLstvec;
152                 int     i;
153                 int     size;
154                 double  result;
155                 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
156                 {
157                         size=(*pLstVec[ipLstvec]).size();
158                         for (i=0;i<size;i++)
159                         {
160                                 result= (*pLstVec[ipLstvec])[i] + bbGetInputk1() ;
161                                 (*pLstVecOut[ipLstvec]).push_back( result );                            
162                         }// for size
163                 } // for pLstVec                
164         } // Type==3
165
166         if (bbGetInputType()==4) // Substraction  k1
167         {
168                 int             ipLstvec;
169                 int     i;
170                 int     size;
171                 double  result;
172                 for (ipLstvec=0 ; ipLstvec<pLstVec.size() ; ipLstvec++)
173                 {
174                         size=(*pLstVec[ipLstvec]).size();
175                         for (i=0;i<size;i++)
176                         {
177                                 result= (*pLstVec[ipLstvec])[i] - bbGetInputk1() ;
178                                 (*pLstVecOut[ipLstvec]).push_back( result );                            
179                         }// for size
180                 } // for pLstVec                
181         } // Type==4
182
183         bbSetOutputOut0( Out0 );
184         bbSetOutputOut1( Out1 );
185         bbSetOutputOut2( Out2 );
186         bbSetOutputOut3( Out3 );
187         bbSetOutputOut4( Out4 );
188         bbSetOutputOut5( Out5 );
189         bbSetOutputOut6( Out6 );
190         bbSetOutputOut7( Out7 );
191         bbSetOutputOut8( Out8 );
192         bbSetOutputOut9( Out9 );
193 }
194 //===== 
195 // 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)
196 //===== 
197 void VectorFilterDouble::bbUserSetDefaultValues()
198 {
199
200 //  SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX 
201 //    Here we initialize the input 'In' to 0
202    bbSetInputType(0);
203    bbSetInputk1(0);
204   
205 }
206 //===== 
207 // 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)
208 //===== 
209 void VectorFilterDouble::bbUserInitializeProcessing()
210 {
211
212 //  THE INITIALIZATION METHOD BODY :
213 //    Here does nothing 
214 //    but this is where you should allocate the internal/output pointers 
215 //    if any 
216
217   
218 }
219 //===== 
220 // 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)
221 //===== 
222 void VectorFilterDouble::bbUserFinalizeProcessing()
223 {
224
225 //  THE FINALIZATION METHOD BODY :
226 //    Here does nothing 
227 //    but this is where you should desallocate the internal/output pointers 
228 //    if any
229   
230 }
231 }
232 // EO namespace bbstd
233
234