traccc_simulate qop assertion failure

[beomki-yeo]

Reported by Yushi in the mattermost (https://mattermost.web.cern.ch/acts/pl/fei6midmijbzmpa73rjy9s4oew)
traccc_simulate_toy_detector is OK but traccc_simulate fails with toy geometry

[beomki@pc-38 traccc_build]$  ./bin/traccc_simulate_toy_detector --output-directory=detray_simulation/toy_detector/n_particles_2000/ --gen-events=10 --gen-nparticles=10 --gen-mom-gev=10:100 --gen-eta=-3:3 --constraint-step-size-mm=1

Running Toy-Detector Simulation

>>> Particle Generation Options <<<
  Number of events to generate   : 10
  Number of particles to generate: 10
  Vertex                         : 0:0:0 mm
  Vertex standard deviation      : 0:0:0 mm
  Momentum range                 : 10:100 GeV
  Phi range                      : -180:180 deg
  Eta range                      : -3:3
  Theta range                    : 5.70047:174.3 deg
  Charge                         : -1
>>> Output Data Options <<<
  Output data format: csv
  Output directory  : detray_simulation/toy_detector/n_particles_2000/
>>> Track Propagation Options <<<
Navigation
----------------------------
  Min. mask tolerance   : 1e-05 [mm]
  Max. mask tolerance   : 1 [mm]
  Mask tolerance scalor : 0.01
  Path tolerance        : 1 [um]
  Overstep tolerance    : -100 [um]
  Search window         : 0 x 0

Parameter Transport
----------------------------
  Min. Stepsize         : 0.0001 [mm]
  Runge-Kutta tolerance : 0.0001 [mm]
  Max. step updates     : 10000
  Stepsize  constraint  : 1 [mm]
  Path limit            : 5 [m]
  Use Bethe energy loss : true
  Do cov. transport     : true
  Use eloss gradient    : false
  Use B-field gradient  : false



WARNING: No entries in volume finder

[beomki@pc-38 traccc_build]$ ./bin/traccc_simulate --detector-file=detray_json/toy_detector_geometry.json --output-directory=detray_simulation/toy_detector/n_particles_2000/ --gen-events=10 --gen-nparticles=10 --gen-mom-gev=10:100 --gen-eta=-3:3 --constraint-step-size-mm=1 

Running Detector Simulation

>>> Detector Options <<<
  Detector file       : detray_json/toy_detector_geometry.json
  Material file       : 
  Surface grid file   : 
  Use detray::detector: no
  Digitization file   : tml_detector/default-geometric-config-generic.json
>>> Particle Generation Options <<<
  Number of events to generate   : 10
  Number of particles to generate: 10
  Vertex                         : 0:0:0 mm
  Vertex standard deviation      : 0:0:0 mm
  Momentum range                 : 10:100 GeV
  Phi range                      : -180:180 deg
  Eta range                      : -3:3
  Theta range                    : 5.70047:174.3 deg
  Charge                         : -1
>>> Output Data Options <<<
  Output data format: csv
  Output directory  : detray_simulation/toy_detector/n_particles_2000/
>>> Track Propagation Options <<<
Navigation
----------------------------
  Min. mask tolerance   : 1e-05 [mm]
  Max. mask tolerance   : 1 [mm]
  Mask tolerance scalor : 0.01
  Path tolerance        : 1 [um]
  Overstep tolerance    : -100 [um]
  Search window         : 0 x 0

Parameter Transport
----------------------------
  Min. Stepsize         : 0.0001 [mm]
  Runge-Kutta tolerance : 0.0001 [mm]
  Max. step updates     : 10000
  Stepsize  constraint  : 1 [mm]
  Path limit            : 5 [m]
  Use Bethe energy loss : true
  Do cov. transport     : true
  Use eloss gradient    : false
  Use B-field gradient  : false



WARNING: No material in detector

WARNING: No entries in volume finder

Detector check: OK
traccc_simulate: /home/beomki/projects/traccc/traccc_build/_deps/detray-src/core/include/detray/tracks/detail/track_helper.hpp:99: detray::detail::track_helper<matrix_operator_t>::scalar_type detray::detail::track_helper<matrix_operator_t>::p(const bound_vector&) const [with matrix_operator_t = algebra::cmath::matrix::actor<long unsigned int, std::array, algebra::array::matrix_type, float, algebra::cmath::matrix::determinant::actor<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::matrix::determinant::partial_pivot_lud<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::element_getter<long unsigned int, std::array, float> >, algebra::cmath::matrix::determinant::hard_coded<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::element_getter<long unsigned int, std::array, float>, 2, 4> >, algebra::cmath::matrix::inverse::actor<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::matrix::inverse::partial_pivot_lud<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::element_getter<long unsigned int, std::array, float> >, algebra::cmath::matrix::inverse::hard_coded<long unsigned int, algebra::array::matrix_type, float, algebra::cmath::element_getter<long unsigned int, std::array, float>, 2, 4> >, algebra::cmath::element_getter<long unsigned int, std::array, float>, algebra::cmath::block_getter<long unsigned int, std::array, float> >; detray::detail::track_helper<matrix_operator_t>::scalar_type = float; detray::detail::track_helper<matrix_operator_t>::bound_vector = std::array<std::array<float, 6>, 1>]: Assertion `qop(bound_vec) != 0.f' failed.
Aborted (core dumped)

I don't have a good guess on why. The geometry built by json input can cause a such problem? @niermann999