2424 * Created on: 12/11/2018
2525 * Author: Andrea Contu
2626 */
27- #include < Evolve.h>
28- #include < simulation.h>
27+ #include < simulation/ Evolve.h>
28+ #include < simulation/simulation .h>
2929
3030void signal_simulation (cfg::Setting const & root, int icfg=-1 ){
3131
@@ -53,14 +53,12 @@ void signal_simulation(cfg::Setting const& root, int icfg=-1){
5353
5454 // efield vecs
5555 VecDev_t<double > _xvec_ef, _yvec_ef, _zvec_ef;
56-
5756 hydra::multiarray<double ,3 ,hydra::host::sys_t > _vec_ef;
58- VecDev_t<double > efieldvecx, efieldvecy, efieldvecz;
59- // wfield vecs
6057
61- hydra::multiarray< double , 3 ,hydra::host:: sys_t > _vec_wf;
58+ // wfield vecs
6259 VecDev_t<double > _xvec_wf, _yvec_wf, _zvec_wf;
63- VecDev_t<double > wfieldvecx, wfieldvecy, wfieldvecz;
60+ hydra::multiarray<double ,3 ,hydra::host::sys_t > _vec_wf;
61+
6462
6563 // emob vecs
6664 VecDev_t<double > _xvec_emob, _yvec_emob, _zvec_emob;
@@ -85,7 +83,6 @@ void signal_simulation(cfg::Setting const& root, int icfg=-1){
8583 pm.prepare ();
8684
8785 pm.fillspacepoints <VecDev_t<double > >(_xvec,_yvec,_zvec);
88- pm.fillallquantities <VecDev_t<double > >(efieldvecx,efieldvecy,efieldvecz,emobvec,hmobvec,wfieldvecx,wfieldvecy,wfieldvecz);
8986 pm.fillallquantities <hydra::multiarray<double ,8 ,hydra::host::sys_t > >(_vec_map);
9087
9188 _xmin = *(pm.getx ().begin ());
@@ -115,7 +112,6 @@ void signal_simulation(cfg::Setting const& root, int icfg=-1){
115112 efieldm.prepare ();
116113
117114 efieldm.fillspacepoints <VecDev_t<double > >(_xvec_ef,_yvec_ef,_zvec_ef);
118- efieldm.fillvector <VecDev_t<double > >(efieldvecx,efieldvecy,efieldvecz);
119115 efieldm.fillvector < hydra::multiarray<double ,3 ,hydra::host::sys_t > >(_vec_ef);
120116
121117 for (auto cf:configlist){
@@ -132,7 +128,6 @@ void signal_simulation(cfg::Setting const& root, int icfg=-1){
132128 maps::physmap<double ,6 > wfieldm (wfieldmap_path); // Weighting field is a vector
133129 wfieldm.prepare ();
134130 wfieldm.fillspacepoints <VecDev_t<double > >(_xvec_wf,_yvec_wf,_zvec_wf);
135- wfieldm.fillvector <VecDev_t<double > >(wfieldvecx,wfieldvecy,wfieldvecz);
136131 wfieldm.fillvector < hydra::multiarray<double ,3 ,hydra::host::sys_t > >(_vec_wf);
137132
138133 std::string emobmap_path=root[" PhysicsMaps" ][" emob" ];
0 commit comments