Handle arbitrary systems

1 files changed, 0 insertions, 22 deletions

diff --git a/packages/sundials/src/helpers.c b/packages/sundials/src/helpers.c
index 420d3be..f0ca592 100644
--- a/packages/sundials/src/helpers.c
+++ b/packages/sundials/src/helpers.c
@@ -42,25 +42,3 @@ int check_flag(void *flagvalue, const char *funcname, int opt)
 
   return 0;
 }
-
-/* check the computed solution */
-int check_ans(N_Vector y, realtype t, realtype rtol, realtype atol)
-{
-  int      passfail=0;     /* answer pass (0) or fail (1) flag     */  
-  realtype ans, err, ewt;  /* answer data, error, and error weight */
-  realtype ONE=RCONST(1.0);
-
-  /* compute solution error */
-  ans  = atan(t);
-  ewt = ONE / (rtol * SUNRabs(ans) + atol);
-  err = ewt * SUNRabs(NV_Ith_S(y,0) - ans);
-
-  /* is the solution within the tolerances? */
-  passfail = (err < ONE) ? 0 : 1; 
-
-  if (passfail) {
-    fprintf(stdout, "\nSUNDIALS_WARNING: check_ans error=%g \n\n", err);
-  }
-
-  return(passfail);
-}