tadd residual function - granular - granular dynamics simulation
  HTML git clone git://src.adamsgaard.dk/granular
   DIR Log
   DIR Files
   DIR Refs
   DIR README
   DIR LICENSE
       ---
   DIR commit ad2bc6049831e53de63e8d2d20402a9cdb8e99a9
   DIR parent 6e12b2e9d44b7138f7e8fe80010783ef2bc6f4d5
  HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Thu, 18 Mar 2021 08:16:04 +0100
       
       add residual function
       
       Diffstat:
         M arrays.c                            |       6 ++++++
         M arrays.h                            |       2 ++
       
       2 files changed, 8 insertions(+), 0 deletions(-)
       ---
   DIR diff --git a/arrays.c b/arrays.c
       t@@ -314,3 +314,9 @@ euclidean_distance(const double *a, const double *b, const int n)
        
                return sqrt(out);
        }
       +
       +double
       +residual(const double new_val, const double old_val)
       +{
       +        return (new_val - old_val) / (old_val + 1e-16);
       +}
   DIR diff --git a/arrays.h b/arrays.h
       t@@ -54,4 +54,6 @@ double * normalize(const double *in, const int n);
        double euclidean_norm(const double *a, const int n);
        double euclidean_distance(const double *a, const double *b, const int n);
        
       +double residual(const double new_val, const double old_val);
       +
        #endif