Git Product home page Git Product logo

Comments (22)

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024
  • is any architecture being favored during compilation? (has given me trouble before)

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

I'm guessing it's gonna be a case of NaN values being returned from the SDF.

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

Found it

IF DEBUG_MODE:
print(f"Volume[0] does not have cached intersections...")
print(state.pos)
self.sdf = -self.SDF(state.pos)

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

Found it

IF DEBUG_MODE:
print(f"Volume[0] does not have cached intersections...")
print(state.pos)
self.sdf = -self.SDF(state.pos)

line 499 prints the correct value: {'x': 0.0, 'y': 0.0, 'z': 10.0}

when printed on the first line of the SDF however: {'x': 4.64496183116417e-310, 'y': 4.6449602498213e-310, 'z': 6.9156690045375e-310}

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

cdef double SDF(self, double3& pos)

from declaration file of the abstract class Volume. Declaration differs from almost all child classes

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

cdef double SDF(self, double3& pos)

from declaration file of the abstract class Volume. Declaration differs from almost all child classes

yup, that was the problem

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

<Electron energy = 1736163663.0886524
position: -1.1889099799771166, 0.7227456741028393, 0.0590154505296173
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736163663.0886524
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961731
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736163663.0886524
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961731
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164134.7492294
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961731
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164134.7492294
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961731
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 7.068983324437212e-18. s_max = 4.73998109816514e-12 Truncate?

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164363.4308345 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.4425162366595506e+19

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419032.1423253

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -243.35091576525065

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 6.251944704487478e-09

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -243.35105114074227 eV

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442518392984275e+19 eV/cm

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 1.8701641667453615e-18cm

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164363.4308345
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164427.81158
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164427.81158
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961732
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 1.284761120278323e-18. s_max = 3.851031170734386e-12 Truncate?

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164606.7818856 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442516719183089e+19

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419156.1063347

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -44.228116930492384

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 1.136267389213581e-09

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -44.22808045995023 eV

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.4425139243292885e+19 eV/cm

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 2.817301370675297e-19cm

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164606.7818856
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164616.4804847
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164616.4804847
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 4.719132462611445e-19. s_max = 3.515137022820144e-12 Truncate?

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164651.009966 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442516806879812e+19

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419178.6362976

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -16.245692892439383

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 4.173691315039781e-10

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -16.245678903713024 eV

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442513858727985e+19 eV/cm

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 3.2269981507154396e-19cm

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164651.009966
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164662.1189518
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164662.1189518
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 1.380944328260702e-18. s_max = 4.247347438615072e-12 Truncate?

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164667.2556448 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.4425168390922166e+19

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419186.9119121

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -47.53924168971775

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 1.2213337101341633e-09

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -47.53928807314909 eV

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.4425202450430984e+19 eV/cm

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 9.665391219945707e-19cm

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164667.2556448
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164700.5289497
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164700.5289497
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 7.217491997436703e-19. s_max = 3.089193670541057e-12 Truncate?

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164714.7949328 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442516933354499e+19

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419211.1286176

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -24.846338595315068

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 6.383288715726928e-10

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -24.846306231842487 eV

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.44251247395448e+19 eV/cm

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 4.362291804934386e-19cm

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164714.7949328
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164729.8121767
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164729.8121767
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 1.841414159580243e-19. s_max = 4.0071761829951894e-12 Truncate?

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164739.641239 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442516982620477e+19

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419223.7854265

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -6.33909952796549

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 1.6285820906028594e-10

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -6.339083956498245 eV

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.4425085326503965e+19 eV/cm

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 7.119549968427503e-20cm

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164739.641239
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164742.09215
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164742.09215
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 2.7311354277578016e-18. s_max = 3.1597766288199433e-12 Truncate?

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164745.9803228 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442516995189797e+19

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419227.0145814

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -94.01980380797414

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 2.4154688800989652e-09

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -94.01981706809644 eV

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.4425175739192295e+19 eV/cm

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 5.798433495801803e-19cm

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164745.9803228
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164765.9415321
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164765.9415321
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 3.675753816607478e-18. s_max = 4.3948767027075e-12 Truncate?

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164840.00014 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.4425171816150233e+19

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419274.9086565

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -126.53846130189149

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 3.2509078806336588e-09

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -126.53850010207555 eV

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442518362637891e+19 eV/cm

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 2.418345234202557e-19cm

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164840.00014
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736164848.325338
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736164848.325338
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 5.123814400561861e-18. s_max = 4.186683393717236e-12 Truncate?

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736164966.53864 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.4425174325192397e+19

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419339.36788

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -176.38821290948368

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 4.531601237777922e-09

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -176.38822456427587 eV

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442517834856269e+19 eV/cm

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 4.4411233010536264e-18cm

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.059015450529617326
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736164966.53864
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736165119.4251018
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736165119.4251018
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 1.1675952089380023e-18. s_max = 3.9846420848120486e-12 Truncate?

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736165142.9268646 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.442517782266939e+19

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419429.2207549

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -40.194673157869374

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 1.0326439241104534e-09

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -40.194681380787706 eV

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442518526377577e+19 eV/cm

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 3.793186821743406e-19cm

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736165142.9268646
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736165155.9849806
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736165155.9849806
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

{delta} :: Proposing displacement.

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Proposed displacement: 1.5112378525393498e-18. s_max = 2.8374204871755587e-12 Truncate?

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: i = 6155 E[6155] = 0.0 < 1736165183.121546 < 3.29455253542804e+180 = E[6156]

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = -3.4425178619661218e+19

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: SP = 884419449.6960881

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: avgW = -52.02463378992628

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: varW = 1.3365682099784455e-09

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

sample_w :: CASE 1

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is -52.024598337102 eV

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of -3.442515567598078e+19 eV/cm

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: The size of the first segment is 9.497262472823337e-19cm

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

MOVE IS DONE...

<Electron energy = 1736165183.121546
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Reducing energy...

<Electron energy = 1736165215.81602
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = 1736165215.81602
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

DOING HINGE...

<Electron energy = 1736165235.1461444
position: -1.1889099799771166, 0.7227456741028393, 0.05901545052961733
direction: 0.13997174426814477, -0.0645060960254009, 0.9880520605626482, norm = 1.0

False :: Temporary variance reduction.

<Electron energy = 1736165235.1461444
position: -1

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

The energy is going up ._.

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024
sample_w :: i = 6155 E[6155] = 0.0 < 1736164363.4308345 < 3.29455253542804e+180 = E[6156]

this kind of stuff should just give off a segmentation fault

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

1736164363.4308345

1736164363.4308345 eV = 1.7361643634308345 GeV

ye this is above the established energy axis

a negative soft stopping power is, somehow, being sampled, causing the electron to gain energy

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

False :: FIRST SEGMENT :: Moving...

<Electron energy = 2147372073.7788603
position: -0.0637114979755362, 1.7193390397608797, 0.10374359959689745
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

MOVE IS DONE...

<Electron energy = 2147372073.7788603
position: -0.0637114979755362, 1.7193390397608797, 0.10374359959689745
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

Reducing energy...

<Electron energy = nan
position: -0.0637114979755362, 1.7193390397608797, 0.10374359959689745
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

Checking if energy is bellow cut off 100000.0...

<Electron energy = nan
position: -0.0637114979755362, 1.7193390397608797, 0.10374359959689745
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

DOING HINGE...

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: Temporary variance reduction.

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

{delta} :: Proposing displacement.

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: Proposed displacement: nan. s_max = nan Truncate?

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: FIRST SEGMENT

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: i = 0 E[0] = 50.0 < nan < 53.78 = E[1]

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: SP = nan

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: SP = nan

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: avgW = nan

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: varW = nan

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

sample_w :: CASE 4

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: FIRST SEGMENT :: Sampled energy loss is nan eV

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: FIRST SEGMENT :: Resulting in a soft SP of nan eV/cm

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: FIRST SEGMENT :: The size of the first segment is nancm

<Electron energy = nan
position: nan, nan, nan
direction: nan, nan, nan, norm = nan

False :: FIRST SEGMENT :: Moving...

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

False :: Proposed displacement: 6.301456148313662e-13. s_max = 4.298253188622311e-12 Truncate?

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

False :: FIRST SEGMENT

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: i = 4999 E[4999] = 0.0 < 2147372073.7788603 < 8.704941672803754e+150 = E[5000]

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: SP = 8.33097868955957e+22

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: SP = 21220464259011.816

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: avgW = -589209423124.9545

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: varW = -476.9898949798073

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

sample_w :: CASE 1

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

False :: FIRST SEGMENT :: Sampled energy loss is nan eV

<Electron energy = 2147372073.7788603
position: -0.0637114979755093, 1.7193390397607833, 0.10374359959680093
direction: -0.1934953022990754, 0.6930826476831923, 0.694403349263701, norm = 1.0

False :: FIRST SEGMENT :: Resulting in a soft SP of nan eV/cm

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

<Electron energy = 150536663.46064615
position: 0.2448921759948367, 0.08463875357937813, 9.959750919041296
direction: 0.2546802863965955, 0.21173508370307542, 0.9435603881311465, norm = 1.0

True :: FIRST SEGMENT :: Moving...

<Electron energy = 150536663.46064615
position: -2.3019106879711186, -2.032712083451376, 0.524147037729831
direction: 0.2546802863965955, 0.21173508370307542, 0.9435603881311465, norm = 1.0

True :: FIRST SEGMENT :: Detected a surface!!

<Electron energy = 2813328456.3082485
position: -2.3019106879711186, -2.032712083451376, 0.524147037729831
direction: 0.2546802863965955, 0.21173508370307542, 0.9435603881311465, norm = 1.0

True :: Temporary variance reduction.

<Electron energy = 2813328456.3082485
position: -2.3019106879711186, -2.032712083451376, 0.524147037729831
direction: 0.2546802863965955, 0.21173508370307542, 0.9435603881311465, norm = 1.0

{delta} :: Proposing displacement.

right here

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

self.state.L = tau
# first segment, pos0
IF DEBUG_MODE: input(str(self) + " MOVING...")
if (<V> self.state.current_region).move(self.state, S_soft):
IF DEBUG_MODE: input(str(self) + " HAS CROSSED SOMETHING...")
self.state.E -= (tau - self.state.L)*S_soft

tau and self.state.L were resulting in numerical error, sometimes it was negative, and likely large enough in absolute value to cause the energy of the particle to go above the max(EAX)

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

gotta redo a lot of CSG.pyx in order to get it to work using only one declaration of SDF

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024
Universe(
    Sphere(radius = 1, material = gold)
    Translate(  # translates **everything** down tree
        Sphere(
            Cube(
                dx = 1, 
                dy = 1, 
                dz = 1, 
                material = air,
            ),
            radius = 1,
            material = water,
            translate = (0, 1, 2) # only translates this surface
        )
        
    )
)

this gonna be the new notation until I finish the transcompiler

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024
Universe(
    Sphere(radius = 1, material = gold)
    *Translate(  # translates **everything** down tree
        Sphere(
            Cube(
                dx = 1, 
                dy = 1, 
                dz = 1, 
                material = air,
            ),
            radius = 1,
            material = water,
            translate = (0, 1, 2) # only translates this surface
        )
        delta_x = 1,
        delta_y = 1,
        delta_z = 1,
    )
    opaque = True # particles that enter this volume get eliminated
)

need the unpack thingy on transforms

also the opaque keyword is important

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

How this thing be lookin:

class Volume               # template/abstract class 
class BVH(Volume)     # segregates the bounding volume hierarchy logic

# this where it gets funky, instances of these two guys are aliases of each other:
class CSGvolume(BVH) # represents a solid volume
class CSGproxy(BVH)    # represents the intersections made on the surface of a solid volume

# CSGvolume is further subdivided 
class Primitive(CSGvolume) # the atom of the csg world
    class Sphere(Primitive)
    class Cube(Primitive)
    class Cylinder(Primitive)
    ...
    etc...

# the stuff that changes the primitives
class Transform(CSGvolume)
class Operation(CSGvolume)

Rn, each primitive holds an instance of a Transform subtype. Which changes in accordance with the transformations that the user asked. This is important. It allows me to choose how to apply the transformations at run time, thus saving a great deal of computational time.

It does beg the question, why am I making Transform subclass CSGvol. Polymorphic behavior from Transform is not expected by the CSG hierarchy. That class is just a tool used by Primitive types.

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

How this thing be lookin:

class Volume               # template/abstract class 
class BVH(Volume)     # segregates the bounding volume hierarchy logic

# this where it gets funky, instances of these two guys are aliases of each other:
class CSGvolume(BVH) # represents a solid volume
class CSGproxy(BVH)    # represents the intersections made on the surface of a solid volume

# CSGvolume is further subdivided 
class Primitive(CSGvolume) # the atom of the csg world
    class Sphere(Primitive)
    class Cube(Primitive)
    class Cylinder(Primitive)
    ...
    etc...

# the stuff that changes the primitives
class Transform(CSGvolume)
class Operation(CSGvolume)

Rn, each primitive holds an instance of a Transform subtype. Which changes in accordance with the transformations that the user asked. This is important. It allows me to choose how to apply the transformations at run time, thus saving a great deal of computational time.

It does beg the question, why am I making Transform subclass CSGvol. Polymorphic behavior from Transform is not expected by the CSG hierarchy. That class is just a tool used by Primitive types.

Can't I hold the transformation matrices as attributes of the primitives, and just use function pointers to decide how transformations are applied?

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024
Universe(
    Sphere(radius = 1, material = gold)
    *Translate(  # translates **everything** down tree
        Sphere(
            Cube(dx = 1, dy = 1, dz = 1, material = air, ),
            Union(
                Sphere(
                    radius = 1, 
                    transforms = ('translate', (0, 0, .5) )
                    ),
                Sphere(
                    radius = 1,
                    transforms = ('translate', (0, 0, -.5) )
                ),
                material = air,
                ),
            radius = 1,
            material = water,
            transforms = [
                ('translate', (1, 2, 3) )
                ('rotate', (0, 0, 1), pi/2)
                ('translate', (1, -2, 3))
                ],
            ),
        delta_x = 1,
        delta_y = 1,
        delta_z = 1,
        ),
    opaque = True # particles that enter this volume get eliminated
    )

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

The "html" equivalent

<Universe opaque=True>
    <Sphere radius = 1, material = gold>
    <Translate delta_x = 1 delta_y = 1 delta_z = 1>
        <Sphere radius = 1 material = water
            transformations:
                translate by 1,2,3
                rotate pi/2 by axis 0,0,1
                translate by 1,-2,3
        >
            <Cube dx=1, dy=1, dz=1 material = air />
            <Union>
                <Sphere radius = 1 transforms: translate by 0,0, .5/>
                <Sphere radius = 1 transforms: translate by 0,0,-.5/>
            </Union>
        </Sphere>
    </Translate>
</Universe>

from montycarlo.

RuiFilipeCampos avatar RuiFilipeCampos commented on June 9, 2024

This issue has been solved once. Now using dockers.

from montycarlo.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.