tadd residual function - granular - granular dynamics simulation
 (HTM) git clone git://src.adamsgaard.dk/granular
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
 (DIR) commit ad2bc6049831e53de63e8d2d20402a9cdb8e99a9
 (DIR) parent 6e12b2e9d44b7138f7e8fe80010783ef2bc6f4d5
 (HTM) 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