@trebco/treb 37.0.0 → 37.1.0

package/build/treb-calculator/src/functions/fp.js

@@ -0,0 +1,463 @@
+/*
+ * This file is part of TREB.
+ *
+ * TREB is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any
+ * later version.
+ *
+ * TREB is distributed in the hope that it will be useful, but WITHOUT ANY
+ * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+ * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
+ * details.
+ *
+ * You should have received a copy of the GNU General Public License along
+ * with TREB. If not, see <https://www.gnu.org/licenses/>.
+ *
+ * Copyright 2022-2026 trebco, llc.
+ * info@treb.app
+ *
+ */
+import { Box, ValueType } from 'treb-base-types';
+import { ArgumentError, ValueError } from '../function-error';
+import { FlattenCellValues } from '../utilities';
+export const FPFunctionLibrary = {
+    MakeArray: {
+        description: 'Create an array using a function',
+        arguments: [
+            { name: 'rows' },
+            { name: 'columns' },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (rows, columns, lambda) {
+            if (rows > 0 && columns > 0 && lambda?.type === ValueType.function) {
+                const apply = this?.apply;
+                if (!apply) {
+                    return ValueError();
+                }
+                const value = [];
+                for (let c = 0; c < columns; c++) {
+                    const values_col = [];
+                    for (let r = 0; r < rows; r++) {
+                        values_col.push(apply(lambda, [
+                            { type: ValueType.number, value: r + 1 },
+                            { type: ValueType.number, value: c + 1 },
+                        ]));
+                    }
+                    value.push(values_col);
+                }
+                return {
+                    type: ValueType.array,
+                    value,
+                };
+            }
+            return ArgumentError();
+        },
+    },
+    ChooseCols: {
+        description: 'Returns one or more columns from an array, by column index',
+        arguments: [{
+                name: 'array',
+            }, {
+                name: 'column',
+                repeat: true,
+            }],
+        fn: function (data, ...args) {
+            if (!Array.isArray(data)) {
+                data = [[data]];
+            }
+            const flat = FlattenCellValues(args);
+            // console.info({flat});
+            const columns = data.length;
+            const result = [];
+            for (let arg of flat) {
+                if (typeof arg === 'number') {
+                    if (arg < 0) {
+                        arg = columns + arg;
+                    }
+                    else {
+                        arg--;
+                    }
+                    const col = data[arg];
+                    if (col) {
+                        const column = [];
+                        for (const value of col) {
+                            column.push(Box(value));
+                        }
+                        result.push(column);
+                    }
+                    else {
+                        return ArgumentError(); // invalid column
+                    }
+                }
+            }
+            if (!result.length) {
+                result.push([]);
+            }
+            return {
+                type: ValueType.array,
+                value: result,
+            };
+        },
+    },
+    ChooseRows: {
+        description: 'Returns one or more rows from an array, by row index',
+        arguments: [{
+                name: 'array',
+            }, {
+                name: 'row',
+                repeat: true,
+            }],
+        fn: function (data, ...args) {
+            if (!Array.isArray(data)) {
+                data = [[data]];
+            }
+            const flat = FlattenCellValues(args);
+            // console.info({flat});
+            const rows = data[0]?.length || 0;
+            const result = [];
+            for (let arg of flat) {
+                if (typeof arg === 'number') {
+                    if (arg < 0) {
+                        arg = rows + arg;
+                    }
+                    else {
+                        arg--;
+                    }
+                    for (const [column_index, col] of data.entries()) {
+                        const column = result[column_index] || [];
+                        if (arg < 0 || arg >= col.length) {
+                            return ArgumentError();
+                        }
+                        column.push(Box(col[arg]));
+                        result[column_index] = column;
+                    }
+                }
+            }
+            if (!result.length) {
+                result.push([]);
+            }
+            return {
+                type: ValueType.array,
+                value: result,
+            };
+        },
+    },
+    Take: {
+        description: 'Returns some number of rows/columns from the start or end of an array',
+        arguments: [{
+                name: 'array',
+            }, {
+                name: 'rows',
+            }, {
+                name: 'columns',
+            }],
+        fn: function (data, rows, columns) {
+            // we need one of rows, columns to be defined. neither can === 0.
+            if ((!rows && !columns) || rows === 0 || columns === 0 || (typeof rows !== 'number' && typeof rows !== 'undefined') || (typeof columns !== 'number' && typeof columns !== 'undefined')) {
+                return ArgumentError();
+            }
+            if (!Array.isArray(data)) {
+                data = [[data]];
+            }
+            const data_columns = data.length;
+            const data_rows = data[0].length;
+            const result = [];
+            // I guess we can compose?
+            // data is column-first
+            let start_column = 0;
+            let end_column = data_columns - 1;
+            let start_row = 0;
+            let end_row = data_rows - 1;
+            if (typeof columns === 'number') {
+                // clip data so it has the first (or last) X columns
+                if (columns > 0) {
+                    end_column = Math.min(columns, data_columns) - 1;
+                }
+                else if (columns < 0) {
+                    end_column = data_columns - 1;
+                    start_column = Math.max(0, end_column + columns + 1);
+                }
+            }
+            if (typeof rows === 'number') {
+                // clip data so it has the first (or last) X columns
+                if (rows > 0) {
+                    end_row = Math.min(rows, data_rows) - 1;
+                }
+                else if (rows < 0) {
+                    end_row = data_rows - 1;
+                    start_row = Math.max(0, end_row + rows + 1);
+                }
+            }
+            for (let c = start_column; c <= end_column; c++) {
+                const column = [];
+                for (let r = start_row; r <= end_row; r++) {
+                    column.push(Box(data[c][r]));
+                }
+                result.push(column);
+            }
+            return {
+                type: ValueType.array,
+                value: result,
+            };
+        },
+    },
+    Reduce: {
+        description: 'Accumulates a value by applying a function to a set of values',
+        arguments: [
+            {
+                name: 'initial value',
+                boxed: true,
+            },
+            {
+                name: 'data',
+                description: 'Input data',
+                boxed: true,
+            },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (initial, data, lambda) {
+            if (!this?.apply) {
+                return ValueError();
+            }
+            if (lambda.type !== ValueType.function) {
+                return ArgumentError();
+            }
+            if (data.type !== ValueType.array) {
+                data = {
+                    type: ValueType.array,
+                    value: [[data]],
+                };
+            }
+            const cols = data.value.length;
+            const rows = data.value[0].length;
+            for (let r = 0; r < rows; r++) {
+                for (let c = 0; c < cols; c++) {
+                    const apply_args = [initial, data.value[c][r]];
+                    const result = this.apply(lambda, apply_args);
+                    initial = result;
+                }
+            }
+            return initial;
+        },
+    },
+    Scan: {
+        description: 'Applies a function to a set of values, iteratively',
+        arguments: [
+            {
+                name: 'initial value',
+                boxed: true,
+            },
+            {
+                name: 'data',
+                description: 'Input data',
+                boxed: true,
+            },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (initial, data, lambda) {
+            if (!this?.apply) {
+                return ValueError();
+            }
+            if (lambda.type !== ValueType.function) {
+                return ArgumentError();
+            }
+            if (data.type !== ValueType.array) {
+                data = {
+                    type: ValueType.array,
+                    value: [[data]],
+                };
+            }
+            const results = [];
+            const cols = data.value.length;
+            const rows = data.value[0].length;
+            for (let i = 0; i < cols; i++) {
+                results.push([]);
+            }
+            for (let r = 0; r < rows; r++) {
+                for (let c = 0; c < cols; c++) {
+                    const apply_args = [initial, data.value[c][r]];
+                    const result = this.apply(lambda, apply_args);
+                    results[c][r] = result;
+                    initial = result;
+                }
+            }
+            return {
+                type: ValueType.array,
+                value: results,
+            };
+        },
+    },
+    ByRow: {
+        description: 'Apply a function to each row in an array',
+        arguments: [
+            {
+                name: 'data',
+                description: 'Input data',
+                repeat: true,
+                boxed: true,
+            },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (data, lambda) {
+            if (!this?.apply) {
+                return ValueError();
+            }
+            if (lambda.type !== ValueType.function) {
+                return ArgumentError();
+            }
+            if (data.type !== ValueType.array) {
+                data = {
+                    type: ValueType.array,
+                    value: [[data]],
+                };
+            }
+            const cols = data.value.length;
+            const rows = data.value[0].length;
+            const value = [[]];
+            for (let r = 0; r < rows; r++) {
+                const args = [];
+                for (let c = 0; c < cols; c++) {
+                    args.push(data.value[c][r]);
+                }
+                value[0].push(this.apply(lambda, [{
+                        type: ValueType.array, value: [args],
+                    }]));
+            }
+            return {
+                type: ValueType.array,
+                value,
+            };
+        }
+    },
+    ByCol: {
+        description: 'Apply a function to each column in an array',
+        arguments: [
+            {
+                name: 'data',
+                description: 'Input data',
+                repeat: true,
+                boxed: true,
+            },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (data, lambda) {
+            if (!this?.apply) {
+                return ValueError();
+            }
+            if (lambda.type !== ValueType.function) {
+                return ArgumentError();
+            }
+            if (data.type !== ValueType.array) {
+                data = {
+                    type: ValueType.array,
+                    value: [[data]],
+                };
+            }
+            const cols = data.value.length;
+            const rows = data.value[0].length;
+            const value = [];
+            for (let c = 0; c < cols; c++) {
+                const args = [];
+                for (let r = 0; r < rows; r++) {
+                    args.push(data.value[c][r]);
+                }
+                value.push([this.apply(lambda, [{
+                            type: ValueType.array, value: [args],
+                        }])]);
+            }
+            return {
+                type: ValueType.array,
+                value,
+            };
+        }
+    },
+    Map: {
+        description: 'Apply a function to a set of values',
+        arguments: [
+            {
+                name: 'data',
+                description: 'Input data',
+                repeat: true,
+                boxed: true,
+            },
+            {
+                name: 'lambda',
+                description: 'Function to apply',
+                boxed: true,
+            },
+        ],
+        fp: true,
+        fn: function (...args) {
+            if (args.length < 2) {
+                return ArgumentError();
+            }
+            const lambda = args[args.length - 1];
+            if (lambda.type !== ValueType.function) {
+                return ArgumentError();
+            }
+            const apply = this?.apply;
+            if (!apply) {
+                return ValueError();
+            }
+            for (let i = 0; i < args.length - 1; i++) {
+                const arg = args[i];
+                if (arg.type !== ValueType.array) {
+                    args[i] = {
+                        type: ValueType.array,
+                        value: [[arg]],
+                    };
+                }
+            }
+            const key = args[0];
+            const results = [];
+            if (key.type === ValueType.array) {
+                for (let r = 0; r < key.value.length; r++) {
+                    const row = key.value[r];
+                    const results_row = [];
+                    for (let c = 0; c < row.length; c++) {
+                        const apply_args = [row[c]];
+                        for (let i = 1; i < args.length - 1; i++) {
+                            const arg = args[i];
+                            if (arg.type === ValueType.array) {
+                                apply_args.push(arg.value[r][c]);
+                            }
+                        }
+                        results_row.push(apply(lambda, apply_args));
+                    }
+                    results.push(results_row);
+                }
+            }
+            return {
+                type: ValueType.array,
+                value: results,
+            };
+        },
+    },
+};
+//# sourceMappingURL=fp.js.map

package/build/treb-calculator/src/functions/fp.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"fp.js","sourceRoot":"","sources":["../../../../treb-calculator/src/functions/fp.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAIH,OAAO,EAAE,GAAG,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AACjD,OAAO,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AAC9D,OAAO,EAAE,iBAAiB,EAAE,MAAM,cAAc,CAAC;AAEjD,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAE5C,SAAS,EAAE;QACT,WAAW,EAAE,kCAAkC;QAC/C,SAAS,EAAE;YACT,EAAE,IAAI,EAAE,MAAM,EAAE;YAChB,EAAE,IAAI,EAAE,SAAS,EAAE;YACnB;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QACD,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,IAAY,EAAE,OAAe,EAAE,MAAqB;YAC/D,IAAI,IAAI,GAAG,CAAC,IAAI,OAAO,GAAG,CAAC,IAAI,MAAM,EAAE,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACnE,MAAM,KAAK,GAAG,IAAI,EAAE,KAAK,CAAC;gBAC1B,IAAI,CAAC,KAAK,EAAE,CAAC;oBACX,OAAO,UAAU,EAAE,CAAC;gBACtB,CAAC;gBAED,MAAM,KAAK,GAAmB,EAAE,CAAC;gBACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;oBACjC,MAAM,UAAU,GAAiB,EAAE,CAAC;oBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC9B,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE;4BAC5B,EAAE,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;4BACxC,EAAE,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE;yBACzC,CAAC,CAAC,CAAC;oBACN,CAAC;oBACD,KAAK,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;gBACzB,CAAC;gBAED,OAAO;oBACL,IAAI,EAAE,SAAS,CAAC,KAAK;oBACrB,KAAK;iBACN,CAAC;YAEJ,CAAC;YACD,OAAO,aAAa,EAAE,CAAC;QACzB,CAAC;KACF;IAED,UAAU,EAAE;QACV,WAAW,EAAE,4DAA4D;QACzE,SAAS,EAAE,CAAC;gBACV,IAAI,EAAE,OAAO;aACd,EAAE;gBACD,IAAI,EAAE,QAAQ;gBACd,MAAM,EAAE,IAAI;aACb,CAAC;QACF,EAAE,EAAE,UAAS,IAA6B,EAAE,GAAG,IAAc;YAE3D,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;gBACzB,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YAClB,CAAC;YAED,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;YACrC,wBAAwB;YAExB,MAAM,OAAO,GAAG,IAAI,CAAC,MAAM,CAAC;YAC5B,MAAM,MAAM,GAAmB,EAAE,CAAC;YAElC,KAAK,IAAI,GAAG,IAAI,IAAI,EAAE,CAAC;gBACrB,IAAI,OAAO,GAAG,KAAK,QAAQ,EAAE,CAAC;oBAE5B,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;wBACZ,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;oBACtB,CAAC;yBACI,CAAC;wBACJ,GAAG,EAAE,CAAC;oBACR,CAAC;oBAED,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;oBACtB,IAAI,GAAG,EAAE,CAAC;wBACR,MAAM,MAAM,GAAiB,EAAE,CAAC;wBAChC,KAAK,MAAM,KAAK,IAAI,GAAG,EAAE,CAAC;4BACxB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;wBAC1B,CAAC;wBACD,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;oBACtB,CAAC;yBACI,CAAC;wBACJ,OAAO,aAAa,EAAE,CAAC,CAAC,iBAAiB;oBAC3C,CAAC;gBACH,CAAC;YACH,CAAC;YAED,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC;gBACnB,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClB,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK,EAAE,MAAM;aACd,CAAC;QAEJ,CAAC;KACF;IAED,UAAU,EAAE;QACV,WAAW,EAAE,sDAAsD;QACnE,SAAS,EAAE,CAAC;gBACV,IAAI,EAAE,OAAO;aACd,EAAE;gBACD,IAAI,EAAE,KAAK;gBACX,MAAM,EAAE,IAAI;aACb,CAAC;QACF,EAAE,EAAE,UAAS,IAA6B,EAAE,GAAG,IAAc;YAE3D,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;gBACzB,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YAClB,CAAC;YAED,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;YACrC,wBAAwB;YAExB,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,IAAI,CAAC,CAAC;YAClC,MAAM,MAAM,GAAmB,EAAE,CAAC;YAElC,KAAK,IAAI,GAAG,IAAI,IAAI,EAAE,CAAC;gBACrB,IAAI,OAAO,GAAG,KAAK,QAAQ,EAAE,CAAC;oBAE5B,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;wBACZ,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC;oBACnB,CAAC;yBACI,CAAC;wBACJ,GAAG,EAAE,CAAC;oBACR,CAAC;oBAED,KAAK,MAAM,CAAC,YAAY,EAAE,GAAG,CAAC,IAAI,IAAI,CAAC,OAAO,EAAE,EAAE,CAAC;wBACjD,MAAM,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC,IAAI,EAAE,CAAC;wBAC1C,IAAI,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,CAAC,MAAM,EAAE,CAAC;4BACjC,OAAO,aAAa,EAAE,CAAC;wBACzB,CAAC;wBACD,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBAC3B,MAAM,CAAC,YAAY,CAAC,GAAG,MAAM,CAAC;oBAChC,CAAC;gBAEH,CAAC;YACH,CAAC;YAED,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC;gBACnB,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClB,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK,EAAE,MAAM;aACd,CAAC;QAEJ,CAAC;KACF;IAED,IAAI,EAAE;QACJ,WAAW,EAAE,uEAAuE;QACpF,SAAS,EAAE,CAAC;gBACV,IAAI,EAAE,OAAO;aACd,EAAE;gBACD,IAAI,EAAE,MAAM;aACb,EAAE;gBACD,IAAI,EAAE,SAAS;aAChB,CAAC;QACF,EAAE,EAAE,UAAS,IAA6B,EAAE,IAAa,EAAE,OAAgB;YAEzE,iEAAiE;YAEjE,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,OAAO,CAAC,IAAI,IAAI,KAAK,CAAC,IAAI,OAAO,KAAK,CAAC,IAAI,CAAC,OAAO,IAAI,KAAK,QAAQ,IAAI,OAAO,IAAI,KAAK,WAAW,CAAC,IAAI,CAAC,OAAO,OAAO,KAAK,QAAQ,IAAI,OAAO,OAAO,KAAK,WAAW,CAAC,EAAE,CAAC;gBACvL,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;gBACzB,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YAClB,CAAC;YAED,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC;YACjC,MAAM,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;YACjC,MAAM,MAAM,GAAmB,EAAE,CAAC;YAElC,0BAA0B;YAE1B,uBAAuB;YAEvB,IAAI,YAAY,GAAG,CAAC,CAAC;YACrB,IAAI,UAAU,GAAG,YAAY,GAAG,CAAC,CAAC;YAElC,IAAI,SAAS,GAAG,CAAC,CAAC;YAClB,IAAI,OAAO,GAAG,SAAS,GAAG,CAAC,CAAC;YAE5B,IAAI,OAAO,OAAO,KAAK,QAAQ,EAAE,CAAC;gBAEhC,oDAAoD;gBAEpD,IAAI,OAAO,GAAG,CAAC,EAAE,CAAC;oBAChB,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,YAAY,CAAC,GAAG,CAAC,CAAC;gBACnD,CAAC;qBACI,IAAI,OAAO,GAAG,CAAC,EAAE,CAAC;oBACrB,UAAU,GAAG,YAAY,GAAG,CAAC,CAAC;oBAC9B,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,UAAU,GAAG,OAAO,GAAG,CAAC,CAAC,CAAC;gBACvD,CAAC;YAEH,CAAC;YAED,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE,CAAC;gBAE7B,oDAAoD;gBAEpD,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;oBACb,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC;gBAC1C,CAAC;qBACI,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;oBAClB,OAAO,GAAG,SAAS,GAAG,CAAC,CAAC;oBACxB,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;gBAC9C,CAAC;YAEH,CAAC;YAED,KAAK,IAAI,CAAC,GAAG,YAAY,EAAE,CAAC,IAAI,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;gBAChD,MAAM,MAAM,GAAiB,EAAE,CAAC;gBAChC,KAAK,IAAI,CAAC,GAAG,SAAS,EAAE,CAAC,IAAI,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC1C,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/B,CAAC;gBACD,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YACtB,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK,EAAE,MAAM;aACd,CAAA;QAEH,CAAC;KACF;IAED,MAAM,EAAE;QACN,WAAW,EAAE,+DAA+D;QAC5E,SAAS,EAAE;YACT;gBACE,IAAI,EAAE,eAAe;gBACrB,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,MAAM;gBACZ,WAAW,EAAE,YAAY;gBACzB,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QACD,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,OAAmB,EAAE,IAAgB,EAAE,MAAqB;YAEvE,IAAI,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC;gBACjB,OAAO,UAAU,EAAE,CAAC;YACtB,CAAC;YAED,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACvC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gBAClC,IAAI,GAAG;oBACL,IAAI,EAAE,SAAS,CAAC,KAAK;oBACrB,KAAK,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;iBAChB,CAAA;YACH,CAAC;YAED,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;YAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;YAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,MAAM,UAAU,GAAiB,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC7D,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,UAAU,CAAC,CAAA;oBAC7C,OAAO,GAAG,MAAM,CAAC;gBACnB,CAAC;YACH,CAAC;YAED,OAAO,OAAO,CAAC;QAEjB,CAAC;KACF;IAED,IAAI,EAAE;QACJ,WAAW,EAAE,oDAAoD;QACjE,SAAS,EAAE;YACT;gBACE,IAAI,EAAE,eAAe;gBACrB,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,MAAM;gBACZ,WAAW,EAAE,YAAY;gBACzB,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QACD,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,OAAmB,EAAE,IAAgB,EAAE,MAAqB;YAEvE,IAAI,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC;gBACjB,OAAO,UAAU,EAAE,CAAC;YACtB,CAAC;YAED,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACvC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gBAClC,IAAI,GAAG;oBACL,IAAI,EAAE,SAAS,CAAC,KAAK;oBACrB,KAAK,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;iBAChB,CAAA;YACH,CAAC;YAED,MAAM,OAAO,GAAmB,EAAE,CAAC;YACnC,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;YAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;YAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAAC,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;YAAA,CAAC;YAElD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,MAAM,UAAU,GAAiB,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC7D,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,UAAU,CAAC,CAAA;oBAC7C,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;oBACvB,OAAO,GAAG,MAAM,CAAC;gBACnB,CAAC;YACH,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK,EAAE,OAAO;aACf,CAAA;QAEH,CAAC;KACF;IAED,KAAK,EAAE;QACL,WAAW,EAAE,0CAA0C;QACvD,SAAS,EAAE;YACT;gBACE,IAAI,EAAE,MAAM;gBACZ,WAAW,EAAE,YAAY;gBACzB,MAAM,EAAE,IAAI;gBACZ,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QACD,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,IAAgB,EAAE,MAAqB;YAClD,IAAI,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC;gBAAC,OAAO,UAAU,EAAE,CAAC;YAAC,CAAC;YAE1C,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACvC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gBAClC,IAAI,GAAG;oBACL,IAAI,EAAE,SAAS,CAAC,KAAK;oBACrB,KAAK,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;iBAChB,CAAC;YACJ,CAAC;YAED,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;YAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;YAElC,MAAM,KAAK,GAAmB,CAAC,EAAE,CAAC,CAAC;YACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,IAAI,GAAiB,EAAE,CAAC;gBAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9B,CAAC;gBACD,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;wBAChC,IAAI,EAAE,SAAS,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,IAAI,CAAC;qBACrC,CAAC,CAAC,CAAC,CAAC;YACP,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK;aACN,CAAA;QAEH,CAAC;KACF;IAGD,KAAK,EAAE;QACL,WAAW,EAAE,6CAA6C;QAC1D,SAAS,EAAE;YACT;gBACE,IAAI,EAAE,MAAM;gBACZ,WAAW,EAAE,YAAY;gBACzB,MAAM,EAAE,IAAI;gBACZ,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QACD,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,IAAgB,EAAE,MAAqB;YAClD,IAAI,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC;gBAAC,OAAO,UAAU,EAAE,CAAC;YAAC,CAAC;YAE1C,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACvC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gBAClC,IAAI,GAAG;oBACL,IAAI,EAAE,SAAS,CAAC,KAAK;oBACrB,KAAK,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;iBAChB,CAAC;YACJ,CAAC;YAED,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;YAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;YAElC,MAAM,KAAK,GAAmB,EAAE,CAAC;YACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,IAAI,GAAiB,EAAE,CAAC;gBAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9B,CAAC;gBACD,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;4BAC9B,IAAI,EAAE,SAAS,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,IAAI,CAAC;yBACrC,CAAC,CAAC,CAAC,CAAC,CAAC;YACR,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK;aACN,CAAA;QAEH,CAAC;KACF;IAED,GAAG,EAAE;QACH,WAAW,EAAE,qCAAqC;QAClD,SAAS,EAAE;YACT;gBACE,IAAI,EAAE,MAAM;gBACZ,WAAW,EAAE,YAAY;gBACzB,MAAM,EAAE,IAAI;gBACZ,KAAK,EAAE,IAAI;aACZ;YACD;gBACE,IAAI,EAAE,QAAQ;gBACd,WAAW,EAAE,mBAAmB;gBAChC,KAAK,EAAE,IAAI;aACZ;SACF;QAED,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,UAAS,GAAG,IAAkB;YAEhC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACpB,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;YACrC,IAAI,MAAM,CAAC,IAAI,KAAK,SAAS,CAAC,QAAQ,EAAE,CAAC;gBACvC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YAED,MAAM,KAAK,GAAG,IAAI,EAAE,KAAK,CAAC;YAC1B,IAAI,CAAC,KAAK,EAAE,CAAC;gBACX,OAAO,UAAU,EAAE,CAAC;YACtB,CAAC;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzC,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;gBACpB,IAAI,GAAG,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;oBACjC,IAAI,CAAC,CAAC,CAAC,GAAG;wBACR,IAAI,EAAE,SAAS,CAAC,KAAK;wBACrB,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;qBACf,CAAC;gBACJ,CAAC;YACH,CAAC;YAED,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACpB,MAAM,OAAO,GAAmB,EAAE,CAAC;YACnC,IAAI,GAAG,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gBACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC1C,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBACzB,MAAM,WAAW,GAAiB,EAAE,CAAC;oBACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;wBACpC,MAAM,UAAU,GAAiB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;wBAE1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BACzC,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;4BACpB,IAAI,GAAG,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;gCACjC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;4BAClC,CAAC;wBACH,CAAC;wBAED,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC,CAAC;oBAC9C,CAAC;oBACD,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;gBAC5B,CAAC;YACH,CAAC;YAED,OAAO;gBACL,IAAI,EAAE,SAAS,CAAC,KAAK;gBACrB,KAAK,EAAE,OAAO;aACf,CAAC;QAEJ,CAAC;KAEF;CAEF,CAAC"}

package/build/treb-calculator/src/functions/function-utilities.d.ts

@@ -0,0 +1,2 @@
+import type { UnionValue, Complex } from 'treb-base-types';
+export declare const CoerceComplex: (value: UnionValue) => Complex | false;

package/build/treb-calculator/src/functions/function-utilities.js

@@ -0,0 +1,36 @@
+/*
+ * This file is part of TREB.
+ *
+ * TREB is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any
+ * later version.
+ *
+ * TREB is distributed in the hope that it will be useful, but WITHOUT ANY
+ * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+ * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
+ * details.
+ *
+ * You should have received a copy of the GNU General Public License along
+ * with TREB. If not, see <https://www.gnu.org/licenses/>.
+ *
+ * Copyright 2022-2026 trebco, llc.
+ * info@treb.app
+ *
+ */
+import { ValueType } from 'treb-base-types';
+export const CoerceComplex = (value) => {
+    switch (value.type) {
+        case ValueType.complex:
+            return value.value;
+        case ValueType.number:
+            return { real: value.value, imaginary: 0 };
+        case ValueType.boolean:
+            return { real: value.value ? 1 : 0, imaginary: 0 };
+        // controversial? ...
+        case ValueType.undefined:
+            return { real: 0, imaginary: 0 };
+    }
+    return false;
+};
+//# sourceMappingURL=function-utilities.js.map

package/build/treb-calculator/src/functions/function-utilities.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"function-utilities.js","sourceRoot":"","sources":["../../../../treb-calculator/src/functions/function-utilities.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAIH,OAAO,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AAE5C,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,KAAiB,EAAiB,EAAE;IAEhE,QAAQ,KAAK,CAAC,IAAI,EAAE,CAAC;QACnB,KAAK,SAAS,CAAC,OAAO;YACpB,OAAO,KAAK,CAAC,KAAK,CAAC;QAErB,KAAK,SAAS,CAAC,MAAM;YACnB,OAAO,EAAE,IAAI,EAAE,KAAK,CAAC,KAAK,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC;QAE7C,KAAK,SAAS,CAAC,OAAO;YACpB,OAAO,EAAE,IAAI,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC;QAErD,qBAAqB;QACrB,KAAK,SAAS,CAAC,SAAS;YACtB,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,SAAS,EAAE,CAAC,EAAE,CAAC;IAErC,CAAC;IAED,OAAO,KAAK,CAAC;AAEf,CAAC,CAAC"}

package/build/treb-calculator/src/functions/gamma.d.ts

@@ -0,0 +1,20 @@
+/**
+ * faster approximation for real numbers
+ */
+export declare const gamma_ln: (z: number) => number;
+/**
+ * series representation
+ */
+export declare const gamma_series: (a: number, x: number) => number | false;
+/**
+ * continued fraction
+ */
+export declare const gamma_cf: (a: number, x: number) => number | false;
+/**
+ * regularized gamma function; CDF of gamma distribution with scale 1
+ */
+export declare const gamma_p: (a: number, x: number) => number | false;
+/**
+ * regularized gamma function
+ */
+export declare const gamma_q: (a: number, x: number) => number | false;

package/build/treb-calculator/src/functions/gamma.js

@@ -0,0 +1,142 @@
+/*
+ * This file is part of TREB.
+ *
+ * TREB is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any
+ * later version.
+ *
+ * TREB is distributed in the hope that it will be useful, but WITHOUT ANY
+ * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+ * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
+ * details.
+ *
+ * You should have received a copy of the GNU General Public License along
+ * with TREB. If not, see <https://www.gnu.org/licenses/>.
+ *
+ * Copyright 2022-2026 trebco, llc.
+ * info@treb.app
+ *
+ */
+//
+// functions and constants for gamma distribution. we're
+// returning false instead of throwing so we can return
+// a spreadsheet-style eror. TODO: optional?
+//
+const max_iterations = 1000;
+const epsilon = 3.0e-9;
+const min_value = 1.0e-30;
+/**
+ * faster approximation for real numbers
+ */
+export const gamma_ln = (z) => {
+    let x = z - 1.0;
+    let y = x + 5.5;
+    y -= (x + 0.5) * Math.log(y);
+    let a = 1.0;
+    const coefficients = [
+        76.18009173,
+        -86.50532033,
+        24.01409822,
+        -1.231739516,
+        0.120858003e-2,
+        -0.536382e-5,
+    ];
+    for (const coeff of coefficients) {
+        a += coeff / (++x);
+    }
+    return (-y + Math.log(2.50662827465 * a));
+};
+/**
+ * series representation
+ */
+export const gamma_series = (a, x) => {
+    const gamma_ln_a = gamma_ln(a);
+    if (x === 0) {
+        return 0;
+    }
+    let ap = a;
+    let sum = 1.0 / a;
+    let del = sum;
+    for (let n = 1; n < max_iterations; n++) {
+        ++ap;
+        del *= x / ap;
+        sum += del;
+        if (Math.abs(del) < Math.abs(sum) * epsilon) {
+            return sum * Math.exp(-x + a * Math.log(x) - (gamma_ln_a));
+        }
+    }
+    // throw new Error('too many iterations');
+    console.error('too many iterations');
+    return false;
+};
+/**
+ * continued fraction
+ */
+export const gamma_cf = (a, x) => {
+    const gamma_ln_a = gamma_ln(a);
+    let b = x + 1.0 - a;
+    let c = 1.0 / min_value;
+    let d = 1.0 / b;
+    let h = d;
+    for (let i = 1; i <= max_iterations; i++) {
+        const an = -i * (i - a);
+        b += 2.0;
+        d = an * d + b;
+        if (Math.abs(d) < min_value)
+            d = min_value;
+        c = b + an / c;
+        if (Math.abs(c) < min_value)
+            c = min_value;
+        d = 1.0 / d;
+        const del = d * c;
+        h *= del;
+        if (Math.abs(del - 1.0) < epsilon) {
+            return Math.exp(-x + a * Math.log(x) - (gamma_ln_a)) * h;
+        }
+    }
+    // throw new Error('too many iterations');
+    console.error('too many iterations');
+    return false;
+};
+/**
+ * regularized gamma function; CDF of gamma distribution with scale 1
+ */
+export const gamma_p = (a, x) => {
+    if (x < 0.0 || a <= 0.0) {
+        // throw new Error('invalid parameter');
+        console.error('invalid parameter');
+        return false;
+    }
+    if (x < (a + 1.0)) {
+        return gamma_series(a, x);
+    }
+    else {
+        const cf = gamma_cf(a, x);
+        if (cf === false) {
+            return false;
+        }
+        return 1 - cf;
+    }
+};
+/**
+ * regularized gamma function
+ */
+export const gamma_q = (a, x) => {
+    if (x < 0.0 || a <= 0.0) {
+        // throw new Error('invalid parameter');
+        console.error('invalid parameter');
+        return false;
+    }
+    if (x < (a + 1.0)) {
+        const series = gamma_series(a, x);
+        if (series === false) {
+            return false;
+        }
+        return 1 - series;
+    }
+    else {
+        return gamma_cf(a, x);
+    }
+};
+//# sourceMappingURL=gamma.js.map