diff --git a/mcstas-comps/examples/Tests_samples/Samples_Phonon/Samples_Phonon.instr b/mcstas-comps/examples/Tests_samples/Samples_Phonon/Samples_Phonon.instr index 537023ac7c..a225e23f60 100644 --- a/mcstas-comps/examples/Tests_samples/Samples_Phonon/Samples_Phonon.instr +++ b/mcstas-comps/examples/Tests_samples/Samples_Phonon/Samples_Phonon.instr @@ -17,7 +17,7 @@ * Simple test instrument for the Phonon_simple component. * Refer to the component documentation for further instructions. * -* %Example: E=10 -n 1e5 Detector: mon1_I=4.9e-25 +* %Example: E=10 -n 1e5 Detector: mon1_I=2.86265e-25 * * %Parameters * E: [meV] Mean energy at source diff --git a/mcstas-comps/samples/Phonon_simple.comp b/mcstas-comps/samples/Phonon_simple.comp index 557205eb10..d11400a7a7 100644 --- a/mcstas-comps/samples/Phonon_simple.comp +++ b/mcstas-comps/samples/Phonon_simple.comp @@ -319,9 +319,9 @@ double zridd_gpu(double x1, double x2, struct neutron_params *neutron, struct ph // Then in energy gain for the neutron for (int i=0; ie_steps_high_; i++){ - root = zridd(f, brack_low+range_high*i/ phonon->e_steps_high_, - brack_low+range_high*(i+1)/ phonon->e_steps_high_, - neutron, phonon, ROOTACC); + root = zridd(f, brack_mid+range_high*i/ phonon->e_steps_high_, + brack_mid+range_high*(i+1)/ phonon->e_steps_high_, + neutron, phonon, ROOTACC); if (root != UNUSED){ list[(*index)++]=root; } @@ -348,9 +348,9 @@ double zridd_gpu(double x1, double x2, struct neutron_params *neutron, struct ph // Then in energy gain for the neutron for (int i=0; ie_steps_high_; i++){ - root = zridd_gpu(brack_low+range_high*i/ phonon->e_steps_high_, - brack_low+range_high*(i+1)/ phonon->e_steps_high_, - neutron, phonon, ROOTACC); + root = zridd_gpu(brack_mid+range_high*i/ phonon->e_steps_high_, + brack_mid+range_high*(i+1)/ phonon->e_steps_high_, + neutron, phonon, ROOTACC); if (root != UNUSED){ list[(*index)++]=root; }