load_tax_table: Optional deductible

[inverse-tax] / tax.js

diff --git a/tax.js b/tax.js
index 2319d21e0e44d4ed29b5aeb56eace3bf3790d819..10c44f3035e3c81fd21dcfe9554f754e0f93efa0 100644 (file)
--- a/tax.js
+++ b/tax.js
@@ -1,5 +1,13 @@
 "use strict";
 
 "use strict";
 

+function near(a, b, epsilon = 1e-6) {
+  return Math.abs(a - b) < epsilon;
+}
+
+function less_than_or_near(a, b, epsilon = 1e-6) {
+  return a < b || near(a, b, epsilon);
+}
+

 function parse_table(as_text) {
   function parse_line(line) {
     return line.trim().split(' ').filter(x => x !== '').map(x => parseFloat(x));
 function parse_table(as_text) {
   function parse_line(line) {
     return line.trim().split(' ').filter(x => x !== '').map(x => parseFloat(x));


@@ -7,11 +15,15 @@ function parse_table(as_text) {
   return as_text.trim().split('\n').map(parse_line);
 }
 
   return as_text.trim().split('\n').map(parse_line);
 }
 

-function parse_tax_table(as_text) {
-  return parse_table(as_text).map(([start, rate], i, table) =>
+function tax_table_from_table(table) {
+  return table.map(([start, rate], i, table) =>

     [start, i < table.length - 1 ? table[i+1][0] : Infinity, rate / 100.0]);
 }
 
     [start, i < table.length - 1 ? table[i+1][0] : Infinity, rate / 100.0]);
 }
 

+function parse_tax_table(as_text) {
+  return tax_table_from_table(parse_table(as_text));
+}
+

 function sum(nums) {
   return nums.reduce((total, num) => total + num, 0);
 }
 function sum(nums) {
   return nums.reduce((total, num) => total + num, 0);
 }


@@ -47,3 +59,42 @@ function merge_tax_tables(t1, t2) {
   }
   return merge_tax_tables(t2, t1);
 }
   }
   return merge_tax_tables(t2, t1);
 }

+
+function invert(table) {
+  // Here we solve
+  //   net = m * gross + b
+  // for gross:
+  //   net - b = m * gross
+  //   (net - b) / m = gross
+
+  const ms = table.map(([start, end, rate]) => 1 - rate);
+  const full_brackets = [0].concat(table.map(([start, end, rate]) => (end - start) * rate)).slice(0, table.length);
+  function sum_lower_brackets(remaining_brackets, acc = 0) {
+    if (remaining_brackets.length == 0) return [];
+    return [acc + remaining_brackets[0]].concat(sum_lower_brackets(remaining_brackets.slice(1), acc + remaining_brackets[0]));
+  }
+  const bs = sum_lower_brackets(full_brackets).map((lower_brackets, i) => {
+    const [start, end, rate] = table[i];
+    // Finding b:
+    //   net = gross - lower_brackets - rate * (gross - start)
+    //   net = gross - lower_brackets - rate * gross + rate * start
+    //   net = gross - rate * gross - lower_brackets + rate * start
+    //   net = (1 - rate) * gross - lower_brackets + rate * start
+    //   net = m * gross - lower_brackets + rate * start
+    //                   \_____________________________/  - here is b
+    return rate * start - lower_brackets;
+  });
+  const inverse_table = table.map(([start, end, rate], i) => {
+    const m = ms[i];
+    const b = bs[i];
+    return [m * start + b, m * end + b, m, b];
+  });
+  return function(net) {
+    for (const [start, end, m, b] of inverse_table) {
+      if (less_than_or_near(start, net) && less_than_or_near(net, end)) {
+        return (net - b) / m;
+      }
+    }
+    return net;
+  };
+}