@@ -221,10 +221,11 @@ mod tests {
221221 mod planar_inversion {
222222 use super :: * ;
223223
224- const C_HALF : f64 = 20.0 ; // C/2 = 20 kcal/mol
224+ const C : f64 = 40.0 ;
225+ const C_HALF : f64 = C / 2.0 ;
225226
226227 fn params ( ) -> f64 {
227- C_HALF
228+ PlanarInversion :: precompute ( C )
228229 }
229230
230231 // --------------------------------------------------------------------
@@ -248,7 +249,7 @@ mod tests {
248249 fn sanity_f32_f64_consistency ( ) {
249250 let cos_psi = 0.4 ;
250251 let p64 = params ( ) ;
251- let p32 = C_HALF as f32 ;
252+ let p32 = PlanarInversion :: precompute ( C as f32 ) ;
252253
253254 let e64 = PlanarInversion :: energy ( cos_psi, p64) ;
254255 let e32 = PlanarInversion :: energy ( cos_psi as f32 , p32) ;
@@ -363,13 +364,13 @@ mod tests {
363364
364365 #[ test]
365366 fn precompute_values ( ) {
366- let c_half = PlanarInversion :: precompute ( C_HALF * 2.0 ) ;
367+ let c_half = PlanarInversion :: precompute ( C ) ;
367368 assert_relative_eq ! ( c_half, C_HALF , epsilon = 1e-14 ) ;
368369 }
369370
370371 #[ test]
371372 fn precompute_round_trip ( ) {
372- let p = PlanarInversion :: precompute ( 40.0 ) ;
373+ let p = PlanarInversion :: precompute ( C ) ;
373374 let e = PlanarInversion :: energy ( 0.0 , p) ;
374375 assert_relative_eq ! ( e, 0.0 , epsilon = 1e-14 ) ;
375376 }
0 commit comments