Update Evolve.h

acontu · web-flow · commit b8c60cfa0990 · 2020-06-17T15:06:11.000+02:00
Avoid dot product bug in Hydra
diff --git a/include/simulation/Evolve.h b/include/simulation/Evolve.h
@@ -241,7 +241,7 @@ namespace evolve{
             auto Vdrift = k1 +ref;
             
             //calculate current with Ramo's theorem
-            hydra::get<5>(p)=hydra::get<_tc_charge>(p)*HCHARGE*(Vdrift.dot(wfvec));
+            hydra::get<5>(p)=hydra::get<_tc_charge>(p)*HCHARGE*(Vdrift.get(0)*wfvec.get(0)+Vdrift.get(1)*wfvec.get(1)+Vdrift.get(2)*wfvec.get(2));
         }
         
         size_t _fN;
@@ -373,7 +373,7 @@ namespace evolve{
             auto Vdrift = k1 +ref;
             
             //calculate current with Ramo's theorem
-            hydra::get<5>(p)=hydra::get<_tc_charge>(p)*HCHARGE*(Vdrift.dot(wfvec));
+            hydra::get<5>(p)=hydra::get<_tc_charge>(p)*HCHARGE*(Vdrift.get(0)*wfvec.get(0)+Vdrift.get(1)*wfvec.get(1)+Vdrift.get(2)*wfvec.get(2));
             
         }
         

Original file line number	Diff line number	Diff line change
`@@ -241,7 +241,7 @@ namespace evolve{`
`241`	`241`	`auto Vdrift = k1 +ref;`
`242`	`242`
`243`	`243`	`//calculate current with Ramo's theorem`
`244`		`- hydra::get<5>(p)=hydra::get<_tc_charge>(p)HCHARGE(Vdrift.dot(wfvec));`
	`244`	`+ hydra::get<5>(p)=hydra::get<_tc_charge>(p)HCHARGE(Vdrift.get(0)wfvec.get(0)+Vdrift.get(1)wfvec.get(1)+Vdrift.get(2)*wfvec.get(2));`
`245`	`245`	`}`
`246`	`246`
`247`	`247`	`size_t _fN;`
`@@ -373,7 +373,7 @@ namespace evolve{`
`373`	`373`	`auto Vdrift = k1 +ref;`
`374`	`374`
`375`	`375`	`//calculate current with Ramo's theorem`
`376`		`- hydra::get<5>(p)=hydra::get<_tc_charge>(p)HCHARGE(Vdrift.dot(wfvec));`
	`376`	`+ hydra::get<5>(p)=hydra::get<_tc_charge>(p)HCHARGE(Vdrift.get(0)wfvec.get(0)+Vdrift.get(1)wfvec.get(1)+Vdrift.get(2)*wfvec.get(2));`
`377`	`377`
`378`	`378`	`}`
`379`	`379`