From ce4be21ca210c61f1e73f8d8432e5d3f8d9e61a5 Mon Sep 17 00:00:00 2001
From: "Douglas B. Rumbaugh" <doug@douglasrumbaugh.com>
Date: Mon, 6 Jan 2025 12:47:54 -0500
Subject: cdf: Added "reverse" cdf functionality

---
 src/cdf.c | 43 ++++++++++++++++++++++++++++++++++---------
 1 file changed, 34 insertions(+), 9 deletions(-)

(limited to 'src/cdf.c')

diff --git a/src/cdf.c b/src/cdf.c
index 8fe1ea6..2ab1d5c 100644
--- a/src/cdf.c
+++ b/src/cdf.c
@@ -122,25 +122,50 @@ static int read_data_uint(DistRecord **records, size_t capacity, FILE *file) {
 }
 
 static int print_data_fp(DistRecord *records, long double *freqs, size_t cnt) {
-  for (size_t i=0; i<cnt; i++) {
-    fprintf(stdout, "%Lf\t%lf\n", freqs[i], records[i].data.d);
-  }  
+  if (ARG_REVERSE_CDF) {
+    long double total_freq = 1.0;
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%lf\n", total_freq, records[i].data.d);
+      total_freq -= freqs[i];
+    }
+    
+  } else {
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%lf\n", freqs[i], records[i].data.d);
+    }  
+  }
 
   return 1;
 }
 
 static int print_data_int(DistRecord *records, long double *freqs, size_t cnt) {
-  for (size_t i=0; i<cnt; i++) {
-    fprintf(stdout, "%Lf\t%ld\n", freqs[i], records[i].data.i);
-  }  
+  if (ARG_REVERSE_CDF) {
+    long double total_freq = 1.0;
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%ld\n", total_freq, records[i].data.i);
+      total_freq -= freqs[i];
+    }
+  } else {
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%ld\n", freqs[i], records[i].data.i);
+    }  
+  }
 
   return 1;
 }
 
 static int print_data_uint(DistRecord *records, long double *freqs, size_t cnt) {
-  for (size_t i=0; i<cnt; i++) {
-    fprintf(stdout, "%Lf\t%ld\n", freqs[i], records[i].data.u);
-  }  
+  if (ARG_REVERSE_CDF) {
+    long double total_freq = 1.0;
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%ld\n", total_freq, records[i].data.u);
+      total_freq -= freqs[i];
+    }
+  } else {
+    for (size_t i=0; i<cnt; i++) {
+      fprintf(stdout, "%Lf\t%ld\n", freqs[i], records[i].data.u);
+    }  
+  }
 
   return 1;
 }
-- 
cgit v1.2.3