simple-squiggle

A restricted subset of Squiggle
Log | Files | Refs | README

matrixFromColumns.js (4631B)

      1 "use strict";
      2 
      3 Object.defineProperty(exports, "__esModule", {
      4   value: true
      5 });
      6 exports.createMatrixFromColumns = void 0;
      7 
      8 var _factory = require("../../utils/factory.js");
      9 
     10 function _createForOfIteratorHelper(o, allowArrayLike) { var it = typeof Symbol !== "undefined" && o[Symbol.iterator] || o["@@iterator"]; if (!it) { if (Array.isArray(o) || (it = _unsupportedIterableToArray(o)) || allowArrayLike && o && typeof o.length === "number") { if (it) o = it; var i = 0; var F = function F() {}; return { s: F, n: function n() { if (i >= o.length) return { done: true }; return { done: false, value: o[i++] }; }, e: function e(_e) { throw _e; }, f: F }; } throw new TypeError("Invalid attempt to iterate non-iterable instance.\nIn order to be iterable, non-array objects must have a [Symbol.iterator]() method."); } var normalCompletion = true, didErr = false, err; return { s: function s() { it = it.call(o); }, n: function n() { var step = it.next(); normalCompletion = step.done; return step; }, e: function e(_e2) { didErr = true; err = _e2; }, f: function f() { try { if (!normalCompletion && it.return != null) it.return(); } finally { if (didErr) throw err; } } }; }
     11 
     12 function _unsupportedIterableToArray(o, minLen) { if (!o) return; if (typeof o === "string") return _arrayLikeToArray(o, minLen); var n = Object.prototype.toString.call(o).slice(8, -1); if (n === "Object" && o.constructor) n = o.constructor.name; if (n === "Map" || n === "Set") return Array.from(o); if (n === "Arguments" || /^(?:Ui|I)nt(?:8|16|32)(?:Clamped)?Array$/.test(n)) return _arrayLikeToArray(o, minLen); }
     13 
     14 function _arrayLikeToArray(arr, len) { if (len == null || len > arr.length) len = arr.length; for (var i = 0, arr2 = new Array(len); i < len; i++) { arr2[i] = arr[i]; } return arr2; }
     15 
     16 var name = 'matrixFromColumns';
     17 var dependencies = ['typed', 'matrix', 'flatten', 'size'];
     18 var createMatrixFromColumns = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
     19   var typed = _ref.typed,
     20       matrix = _ref.matrix,
     21       flatten = _ref.flatten,
     22       size = _ref.size;
     23 
     24   /**
     25    * Create a dense matrix from vectors as individual columns.
     26    * If you pass row vectors, they will be transposed (but not conjugated!)
     27    *
     28    * Syntax:
     29    *
     30    *    math.matrixFromColumns(...arr)
     31    *    math.matrixFromColumns(col1, col2)
     32    *    math.matrixFromColumns(col1, col2, col3)
     33    *
     34    * Examples:
     35    *
     36    *    math.matrixFromColumns([1, 2, 3], [[4],[5],[6]])
     37    *    math.matrixFromColumns(...vectors)
     38    *
     39    * See also:
     40    *
     41    *    matrix, matrixFromRows, matrixFromFunction, zeros
     42    *
     43    * @param {... Array | Matrix} cols Multiple columns
     44    * @return { number[][] | Matrix } if at least one of the arguments is an array, an array will be returned
     45    */
     46   return typed(name, {
     47     '...Array': function Array(arr) {
     48       return _createArray(arr);
     49     },
     50     '...Matrix': function Matrix(arr) {
     51       return matrix(_createArray(arr.map(function (m) {
     52         return m.toArray();
     53       })));
     54     } // TODO implement this properly for SparseMatrix
     55 
     56   });
     57 
     58   function _createArray(arr) {
     59     if (arr.length === 0) throw new TypeError('At least one column is needed to construct a matrix.');
     60     var N = checkVectorTypeAndReturnLength(arr[0]); // create an array with empty rows
     61 
     62     var result = [];
     63 
     64     for (var i = 0; i < N; i++) {
     65       result[i] = [];
     66     } // loop columns
     67 
     68 
     69     var _iterator = _createForOfIteratorHelper(arr),
     70         _step;
     71 
     72     try {
     73       for (_iterator.s(); !(_step = _iterator.n()).done;) {
     74         var col = _step.value;
     75         var colLength = checkVectorTypeAndReturnLength(col);
     76 
     77         if (colLength !== N) {
     78           throw new TypeError('The vectors had different length: ' + (N | 0) + ' ≠ ' + (colLength | 0));
     79         }
     80 
     81         var f = flatten(col); // push a value to each row
     82 
     83         for (var _i = 0; _i < N; _i++) {
     84           result[_i].push(f[_i]);
     85         }
     86       }
     87     } catch (err) {
     88       _iterator.e(err);
     89     } finally {
     90       _iterator.f();
     91     }
     92 
     93     return result;
     94   }
     95 
     96   function checkVectorTypeAndReturnLength(vec) {
     97     var s = size(vec);
     98 
     99     if (s.length === 1) {
    100       // 1D vector
    101       return s[0];
    102     } else if (s.length === 2) {
    103       // 2D vector
    104       if (s[0] === 1) {
    105         // row vector
    106         return s[1];
    107       } else if (s[1] === 1) {
    108         // col vector
    109         return s[0];
    110       } else {
    111         throw new TypeError('At least one of the arguments is not a vector.');
    112       }
    113     } else {
    114       throw new TypeError('Only one- or two-dimensional vectors are supported.');
    115     }
    116   }
    117 });
    118 exports.createMatrixFromColumns = createMatrixFromColumns;