mathjs 11.11.2 → 11.12.0
Sign up to get free protection for your applications and to get access to all the features.
- package/HISTORY.md +11 -0
- package/lib/browser/math.js +1 -1
- package/lib/browser/math.js.LICENSE.txt +2 -2
- package/lib/browser/math.js.map +1 -1
- package/lib/cjs/core/create.js +2 -2
- package/lib/cjs/defaultInstance.js +2 -2
- package/lib/cjs/entry/allFactoriesAny.js +2 -2
- package/lib/cjs/entry/allFactoriesNumber.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesDet.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesDistance.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesLsolve.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesLsolveAll.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesLup.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesQr.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesStirlingS2.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesSubtract.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesSubtractScalar.generated.js +17 -0
- package/lib/cjs/entry/dependenciesAny/dependenciesUnitClass.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesUsolve.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny/dependenciesUsolveAll.generated.js +2 -2
- package/lib/cjs/entry/dependenciesAny.generated.js +7 -0
- package/lib/cjs/entry/dependenciesNumber/dependenciesStirlingS2.generated.js +2 -2
- package/lib/cjs/entry/dependenciesNumber/dependenciesSubtractScalar.generated.js +17 -0
- package/lib/cjs/entry/dependenciesNumber.generated.js +7 -0
- package/lib/cjs/entry/impureFunctionsAny.generated.js +4 -3
- package/lib/cjs/entry/impureFunctionsNumber.generated.js +2 -1
- package/lib/cjs/entry/pureFunctionsAny.generated.js +41 -37
- package/lib/cjs/entry/pureFunctionsNumber.generated.js +16 -13
- package/lib/cjs/expression/parse.js +50 -69
- package/lib/cjs/factoriesAny.js +7 -0
- package/lib/cjs/factoriesNumber.js +2 -1
- package/lib/cjs/function/algebra/decomposition/lup.js +3 -3
- package/lib/cjs/function/algebra/decomposition/qr.js +5 -5
- package/lib/cjs/function/algebra/solver/lsolve.js +4 -4
- package/lib/cjs/function/algebra/solver/lsolveAll.js +6 -6
- package/lib/cjs/function/algebra/solver/usolve.js +4 -4
- package/lib/cjs/function/algebra/solver/usolveAll.js +6 -6
- package/lib/cjs/function/arithmetic/subtract.js +7 -37
- package/lib/cjs/function/arithmetic/subtractScalar.js +52 -0
- package/lib/cjs/function/combinatorics/stirlingS2.js +2 -2
- package/lib/cjs/function/geometry/distance.js +10 -10
- package/lib/cjs/function/matrix/det.js +4 -4
- package/lib/cjs/header.js +2 -2
- package/lib/cjs/type/unit/Unit.js +5 -5
- package/lib/cjs/utils/object.js +4 -2
- package/lib/cjs/utils/snapshot.js +3 -3
- package/lib/cjs/utils/string.js +12 -16
- package/lib/cjs/version.js +1 -1
- package/lib/esm/entry/dependenciesAny/dependenciesDet.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesDistance.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesLsolve.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesLsolveAll.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesLup.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesQr.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesStirlingS2.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesSubtract.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesSubtractScalar.generated.js +10 -0
- package/lib/esm/entry/dependenciesAny/dependenciesUnitClass.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesUsolve.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny/dependenciesUsolveAll.generated.js +2 -2
- package/lib/esm/entry/dependenciesAny.generated.js +1 -0
- package/lib/esm/entry/dependenciesNumber/dependenciesStirlingS2.generated.js +2 -2
- package/lib/esm/entry/dependenciesNumber/dependenciesSubtractScalar.generated.js +10 -0
- package/lib/esm/entry/dependenciesNumber.generated.js +1 -0
- package/lib/esm/entry/impureFunctionsAny.generated.js +5 -4
- package/lib/esm/entry/impureFunctionsNumber.generated.js +3 -2
- package/lib/esm/entry/pureFunctionsAny.generated.js +40 -37
- package/lib/esm/entry/pureFunctionsNumber.generated.js +16 -13
- package/lib/esm/expression/parse.js +50 -69
- package/lib/esm/factoriesAny.js +1 -0
- package/lib/esm/factoriesNumber.js +1 -0
- package/lib/esm/function/algebra/decomposition/lup.js +3 -3
- package/lib/esm/function/algebra/decomposition/qr.js +5 -5
- package/lib/esm/function/algebra/solver/lsolve.js +4 -4
- package/lib/esm/function/algebra/solver/lsolveAll.js +6 -6
- package/lib/esm/function/algebra/solver/usolve.js +4 -4
- package/lib/esm/function/algebra/solver/usolveAll.js +6 -6
- package/lib/esm/function/arithmetic/subtract.js +7 -27
- package/lib/esm/function/arithmetic/subtractScalar.js +46 -0
- package/lib/esm/function/combinatorics/stirlingS2.js +2 -2
- package/lib/esm/function/geometry/distance.js +10 -10
- package/lib/esm/function/matrix/det.js +4 -4
- package/lib/esm/type/unit/Unit.js +5 -5
- package/lib/esm/utils/object.js +5 -3
- package/lib/esm/utils/snapshot.js +1 -1
- package/lib/esm/utils/string.js +12 -16
- package/lib/esm/version.js +1 -1
- package/package.json +13 -13
- package/types/index.d.ts +31 -0
@@ -5,7 +5,8 @@ Object.defineProperty(exports, "__esModule", {
|
|
5
5
|
});
|
6
6
|
exports.efimovFactor = exports.e = exports.dotPow = exports.dotMultiply = exports.dotDivide = exports.dot = exports.divideScalar = exports.divide = exports.distance = exports.diff = exports.diag = exports.deuteronMass = exports.det = exports.deepEqual = exports.cumsum = exports.cube = exports.ctranspose = exports.csch = exports.csc = exports.cross = exports.createUnit = exports.count = exports.coulomb = exports.coth = exports.cot = exports.cosh = exports.cos = exports.corr = exports.conj = exports.conductanceQuantum = exports.concat = exports.composition = exports.complex = exports.compareText = exports.compareNatural = exports.compare = exports.combinationsWithRep = exports.combinations = exports.column = exports.clone = exports.classicalElectronRadius = exports.ceil = exports.cbrt = exports.catalan = exports["boolean"] = exports.boltzmann = exports.bohrRadius = exports.bohrMagneton = exports.bitXor = exports.bitOr = exports.bitNot = exports.bitAnd = exports.bin = exports.bignumber = exports.bellNumbers = exports.avogadro = exports.atomicMass = exports.atanh = exports.atan2 = exports.atan = exports.asinh = exports.asin = exports.asech = exports.asec = exports.arg = exports.apply = exports.and = exports.addScalar = exports.add = exports.acsch = exports.acsc = exports.acoth = exports.acot = exports.acosh = exports.acos = exports.abs = exports._true = exports._null = exports._false = exports._NaN = exports._Infinity = exports.Unit = exports.SparseMatrix = exports.Spa = exports.SQRT2 = exports.SQRT1_2 = exports.ResultSet = exports.Range = exports.Matrix = exports.LOG2E = exports.LOG10E = exports.LN2 = exports.LN10 = exports.Index = exports.ImmutableDenseMatrix = exports.Fraction = exports.FibonacciHeap = exports.DenseMatrix = exports.Complex = exports.BigNumber = void 0;
|
7
7
|
exports.ones = exports.oct = exports.numeric = exports.number = exports.nuclearMagneton = exports.nthRoots = exports.nthRoot = exports.not = exports.norm = exports.neutronMass = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.molarVolume = exports.molarPlanckConstant = exports.molarMassC12 = exports.molarMass = exports.mode = exports.mod = exports.min = exports.median = exports.mean = exports.max = exports.matrixFromRows = exports.matrixFromFunction = exports.matrixFromColumns = exports.matrix = exports.map = exports.magneticFluxQuantum = exports.magneticConstant = exports.mad = exports.lyap = exports.lusolve = exports.lup = exports.lsolveAll = exports.lsolve = exports.loschmidt = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.kron = exports.klitzing = exports.kldivergence = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.invmod = exports.inverseConductanceQuantum = exports.inv = exports.intersect = exports.index = exports.im = exports.ifft = exports.identity = exports.i = exports.hypot = exports.hex = exports.hasNumericValue = exports.hartreeEnergy = exports.gravity = exports.gravitationConstant = exports.getMatrixDataType = exports.gcd = exports.gasConstant = exports.gamma = exports.freqz = exports.fraction = exports.format = exports.forEach = exports.floor = exports.flatten = exports.fix = exports.firstRadiation = exports.fineStructure = exports.filter = exports.fft = exports.fermiCoupling = exports.faraday = exports.factorial = exports.expm1 = exports.expm = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.elementaryCharge = exports.electronMass = exports.electricConstant = exports.eigs = void 0;
|
8
|
-
exports.
|
8
|
+
exports.zeta = exports.zeros = exports.xor = exports.xgcd = exports.wienDisplacement = exports.weakMixingAngle = exports.version = exports.variance = exports.vacuumImpedance = exports.usolveAll = exports.usolve = exports.unit = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.transpose = exports.trace = exports.to = exports.thomsonCrossSection = exports.tau = exports.tanh = exports.tan = exports.sylvester = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.stefanBoltzmann = exports.std = exports.squeeze = exports.square = exports.sqrtm = exports.sqrt = exports.splitUnit = exports.speedOfLight = exports.sparse = exports.sort = exports.solveODE = exports.smallerEq = exports.smaller = exports.slu = exports.size = exports.sinh = exports.sin = exports.sign = exports.setUnion = exports.setSymDifference = exports.setSize = exports.setPowerset = exports.setMultiplicity = exports.setIsSubset = exports.setIntersect = exports.setDistinct = exports.setDifference = exports.setCartesian = exports.secondRadiation = exports.sech = exports.sec = exports.schur = exports.sackurTetrode = exports.rydberg = exports.row = exports.round = exports.rotationMatrix = exports.rotate = exports.rightLogShift = exports.rightArithShift = exports.resize = exports.reshape = exports.replacer = exports.reducedPlanckConstant = exports.re = exports.range = exports.randomInt = exports.random = exports.quantumOfCirculation = exports.quantileSeq = exports.qr = exports.protonMass = exports.prod = exports.print = exports.pow = exports.polynomialRoot = exports.planckTime = exports.planckTemperature = exports.planckMass = exports.planckLength = exports.planckConstant = exports.planckCharge = exports.pinv = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = void 0;
|
9
|
+
exports.zpk2tf = void 0;
|
9
10
|
var _configReadonly = require("./configReadonly.js");
|
10
11
|
var _factoriesAny = require("../factoriesAny.js");
|
11
12
|
/**
|
@@ -301,6 +302,9 @@ var square = exports.square = /* #__PURE__ */(0, _factoriesAny.createSquare)({
|
|
301
302
|
var string = exports.string = /* #__PURE__ */(0, _factoriesAny.createString)({
|
302
303
|
typed: typed
|
303
304
|
});
|
305
|
+
var subtractScalar = exports.subtractScalar = /* #__PURE__ */(0, _factoriesAny.createSubtractScalar)({
|
306
|
+
typed: typed
|
307
|
+
});
|
304
308
|
var tan = exports.tan = /* #__PURE__ */(0, _factoriesAny.createTan)({
|
305
309
|
typed: typed
|
306
310
|
});
|
@@ -537,6 +541,15 @@ var leftShift = exports.leftShift = /* #__PURE__ */(0, _factoriesAny.createLeftS
|
|
537
541
|
typed: typed,
|
538
542
|
zeros: zeros
|
539
543
|
});
|
544
|
+
var lsolve = exports.lsolve = /* #__PURE__ */(0, _factoriesAny.createLsolve)({
|
545
|
+
DenseMatrix: DenseMatrix,
|
546
|
+
divideScalar: divideScalar,
|
547
|
+
equalScalar: equalScalar,
|
548
|
+
matrix: matrix,
|
549
|
+
multiplyScalar: multiplyScalar,
|
550
|
+
subtractScalar: subtractScalar,
|
551
|
+
typed: typed
|
552
|
+
});
|
540
553
|
var matrixFromColumns = exports.matrixFromColumns = /* #__PURE__ */(0, _factoriesAny.createMatrixFromColumns)({
|
541
554
|
flatten: flatten,
|
542
555
|
matrix: matrix,
|
@@ -556,6 +569,23 @@ var ones = exports.ones = /* #__PURE__ */(0, _factoriesAny.createOnes)({
|
|
556
569
|
matrix: matrix,
|
557
570
|
typed: typed
|
558
571
|
});
|
572
|
+
var qr = exports.qr = /* #__PURE__ */(0, _factoriesAny.createQr)({
|
573
|
+
addScalar: addScalar,
|
574
|
+
complex: complex,
|
575
|
+
conj: conj,
|
576
|
+
divideScalar: divideScalar,
|
577
|
+
equal: equal,
|
578
|
+
identity: identity,
|
579
|
+
isZero: isZero,
|
580
|
+
matrix: matrix,
|
581
|
+
multiplyScalar: multiplyScalar,
|
582
|
+
sign: sign,
|
583
|
+
sqrt: sqrt,
|
584
|
+
subtractScalar: subtractScalar,
|
585
|
+
typed: typed,
|
586
|
+
unaryMinus: unaryMinus,
|
587
|
+
zeros: zeros
|
588
|
+
});
|
559
589
|
var resize = exports.resize = /* #__PURE__ */(0, _factoriesAny.createResize)({
|
560
590
|
config: _configReadonly.config,
|
561
591
|
matrix: matrix
|
@@ -585,10 +615,10 @@ var smaller = exports.smaller = /* #__PURE__ */(0, _factoriesAny.createSmaller)(
|
|
585
615
|
});
|
586
616
|
var subtract = exports.subtract = /* #__PURE__ */(0, _factoriesAny.createSubtract)({
|
587
617
|
DenseMatrix: DenseMatrix,
|
588
|
-
addScalar: addScalar,
|
589
618
|
concat: concat,
|
590
619
|
equalScalar: equalScalar,
|
591
620
|
matrix: matrix,
|
621
|
+
subtractScalar: subtractScalar,
|
592
622
|
typed: typed,
|
593
623
|
unaryMinus: unaryMinus
|
594
624
|
});
|
@@ -611,7 +641,7 @@ var usolve = exports.usolve = /* #__PURE__ */(0, _factoriesAny.createUsolve)({
|
|
611
641
|
equalScalar: equalScalar,
|
612
642
|
matrix: matrix,
|
613
643
|
multiplyScalar: multiplyScalar,
|
614
|
-
|
644
|
+
subtractScalar: subtractScalar,
|
615
645
|
typed: typed
|
616
646
|
});
|
617
647
|
var xor = exports.xor = /* #__PURE__ */(0, _factoriesAny.createXor)({
|
@@ -695,7 +725,7 @@ var distance = exports.distance = /* #__PURE__ */(0, _factoriesAny.createDistanc
|
|
695
725
|
divideScalar: divideScalar,
|
696
726
|
multiplyScalar: multiplyScalar,
|
697
727
|
sqrt: sqrt,
|
698
|
-
|
728
|
+
subtractScalar: subtractScalar,
|
699
729
|
typed: typed
|
700
730
|
});
|
701
731
|
var dot = exports.dot = /* #__PURE__ */(0, _factoriesAny.createDot)({
|
@@ -761,13 +791,13 @@ var log = exports.log = /* #__PURE__ */(0, _factoriesAny.createLog)({
|
|
761
791
|
divideScalar: divideScalar,
|
762
792
|
typed: typed
|
763
793
|
});
|
764
|
-
var
|
794
|
+
var lsolveAll = exports.lsolveAll = /* #__PURE__ */(0, _factoriesAny.createLsolveAll)({
|
765
795
|
DenseMatrix: DenseMatrix,
|
766
796
|
divideScalar: divideScalar,
|
767
797
|
equalScalar: equalScalar,
|
768
798
|
matrix: matrix,
|
769
799
|
multiplyScalar: multiplyScalar,
|
770
|
-
|
800
|
+
subtractScalar: subtractScalar,
|
771
801
|
typed: typed
|
772
802
|
});
|
773
803
|
var matrixFromRows = exports.matrixFromRows = /* #__PURE__ */(0, _factoriesAny.createMatrixFromRows)({
|
@@ -819,23 +849,6 @@ var partitionSelect = exports.partitionSelect = /* #__PURE__ */(0, _factoriesAny
|
|
819
849
|
isNumeric: isNumeric,
|
820
850
|
typed: typed
|
821
851
|
});
|
822
|
-
var qr = exports.qr = /* #__PURE__ */(0, _factoriesAny.createQr)({
|
823
|
-
addScalar: addScalar,
|
824
|
-
complex: complex,
|
825
|
-
conj: conj,
|
826
|
-
divideScalar: divideScalar,
|
827
|
-
equal: equal,
|
828
|
-
identity: identity,
|
829
|
-
isZero: isZero,
|
830
|
-
matrix: matrix,
|
831
|
-
multiplyScalar: multiplyScalar,
|
832
|
-
sign: sign,
|
833
|
-
sqrt: sqrt,
|
834
|
-
subtract: subtract,
|
835
|
-
typed: typed,
|
836
|
-
unaryMinus: unaryMinus,
|
837
|
-
zeros: zeros
|
838
|
-
});
|
839
852
|
var rightLogShift = exports.rightLogShift = /* #__PURE__ */(0, _factoriesAny.createRightLogShift)({
|
840
853
|
DenseMatrix: DenseMatrix,
|
841
854
|
concat: concat,
|
@@ -879,7 +892,7 @@ var usolveAll = exports.usolveAll = /* #__PURE__ */(0, _factoriesAny.createUsolv
|
|
879
892
|
equalScalar: equalScalar,
|
880
893
|
matrix: matrix,
|
881
894
|
multiplyScalar: multiplyScalar,
|
882
|
-
|
895
|
+
subtractScalar: subtractScalar,
|
883
896
|
typed: typed
|
884
897
|
});
|
885
898
|
var zpk2tf = exports.zpk2tf = /* #__PURE__ */(0, _factoriesAny.createZpk2tf)({
|
@@ -929,7 +942,7 @@ var det = exports.det = /* #__PURE__ */(0, _factoriesAny.createDet)({
|
|
929
942
|
isZero: isZero,
|
930
943
|
matrix: matrix,
|
931
944
|
multiply: multiply,
|
932
|
-
|
945
|
+
subtractScalar: subtractScalar,
|
933
946
|
typed: typed,
|
934
947
|
unaryMinus: unaryMinus
|
935
948
|
});
|
@@ -999,15 +1012,6 @@ var log1p = exports.log1p = /* #__PURE__ */(0, _factoriesAny.createLog1p)({
|
|
999
1012
|
log: log,
|
1000
1013
|
typed: typed
|
1001
1014
|
});
|
1002
|
-
var lsolveAll = exports.lsolveAll = /* #__PURE__ */(0, _factoriesAny.createLsolveAll)({
|
1003
|
-
DenseMatrix: DenseMatrix,
|
1004
|
-
divideScalar: divideScalar,
|
1005
|
-
equalScalar: equalScalar,
|
1006
|
-
matrix: matrix,
|
1007
|
-
multiplyScalar: multiplyScalar,
|
1008
|
-
subtract: subtract,
|
1009
|
-
typed: typed
|
1010
|
-
});
|
1011
1015
|
var max = exports.max = /* #__PURE__ */(0, _factoriesAny.createMax)({
|
1012
1016
|
config: _configReadonly.config,
|
1013
1017
|
larger: larger,
|
@@ -1139,7 +1143,7 @@ var lup = exports.lup = /* #__PURE__ */(0, _factoriesAny.createLup)({
|
|
1139
1143
|
larger: larger,
|
1140
1144
|
matrix: matrix,
|
1141
1145
|
multiplyScalar: multiplyScalar,
|
1142
|
-
|
1146
|
+
subtractScalar: subtractScalar,
|
1143
1147
|
typed: typed,
|
1144
1148
|
unaryMinus: unaryMinus
|
1145
1149
|
});
|
@@ -1214,7 +1218,7 @@ var Unit = exports.Unit = /* #__PURE__ */(0, _factoriesAny.createUnitClass)({
|
|
1214
1218
|
number: number,
|
1215
1219
|
pow: pow,
|
1216
1220
|
round: round,
|
1217
|
-
|
1221
|
+
subtractScalar: subtractScalar
|
1218
1222
|
});
|
1219
1223
|
var vacuumImpedance = exports.vacuumImpedance = /* #__PURE__ */(0, _factoriesAny.createVacuumImpedance)({
|
1220
1224
|
BigNumber: BigNumber,
|
@@ -1547,7 +1551,7 @@ var stirlingS2 = exports.stirlingS2 = /* #__PURE__ */(0, _factoriesAny.createSti
|
|
1547
1551
|
multiplyScalar: multiplyScalar,
|
1548
1552
|
number: number,
|
1549
1553
|
pow: pow,
|
1550
|
-
|
1554
|
+
subtractScalar: subtractScalar,
|
1551
1555
|
typed: typed
|
1552
1556
|
});
|
1553
1557
|
var unit = exports.unit = /* #__PURE__ */(0, _factoriesAny.createUnitFunction)({
|
@@ -4,7 +4,7 @@ Object.defineProperty(exports, "__esModule", {
|
|
4
4
|
value: true
|
5
5
|
});
|
6
6
|
exports.min = exports.median = exports.mean = exports.max = exports.matrix = exports.map = exports.mad = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.index = exports.hypot = exports.hasNumericValue = exports.gcd = exports.gamma = exports.format = exports.forEach = exports.floor = exports.fix = exports.filter = exports.factorial = exports.expm1 = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.e = exports.divideScalar = exports.divide = exports.deepEqual = exports.cumsum = exports.cube = exports.csch = exports.csc = exports.coth = exports.cot = exports.cosh = exports.cos = exports.corr = exports.composition = exports.compareText = exports.compareNatural = exports.compare = exports.combinationsWithRep = exports.combinations = exports.clone = exports.ceil = exports.cbrt = exports.catalan = exports["boolean"] = exports.bitXor = exports.bitOr = exports.bitNot = exports.bitAnd = exports.bellNumbers = exports.atanh = exports.atan2 = exports.atan = exports.asinh = exports.asin = exports.asech = exports.asec = exports.apply = exports.and = exports.addScalar = exports.add = exports.acsch = exports.acsc = exports.acoth = exports.acot = exports.acosh = exports.acos = exports.abs = exports._true = exports._null = exports._false = exports._NaN = exports._Infinity = exports.SQRT2 = exports.SQRT1_2 = exports.ResultSet = exports.Range = exports.LOG2E = exports.LOG10E = exports.LN2 = exports.LN10 = void 0;
|
7
|
-
exports.zeta = exports.xor = exports.xgcd = exports.version = exports.variance = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.tau = exports.tanh = exports.tan = exports.sum = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.std = exports.square = exports.sqrt = exports.smallerEq = exports.smaller = exports.size = exports.sinh = exports.sin = exports.sign = exports.sech = exports.sec = exports.round = exports.rightLogShift = exports.rightArithShift = exports.replacer = exports.range = exports.randomInt = exports.random = exports.quantileSeq = exports.prod = exports.print = exports.pow = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.numeric = exports.number = exports.nthRoot = exports.not = exports.norm = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.mode = exports.mod = void 0;
|
7
|
+
exports.zeta = exports.xor = exports.xgcd = exports.version = exports.variance = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.tau = exports.tanh = exports.tan = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.std = exports.square = exports.sqrt = exports.smallerEq = exports.smaller = exports.size = exports.sinh = exports.sin = exports.sign = exports.sech = exports.sec = exports.round = exports.rightLogShift = exports.rightArithShift = exports.replacer = exports.range = exports.randomInt = exports.random = exports.quantileSeq = exports.prod = exports.print = exports.pow = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.numeric = exports.number = exports.nthRoot = exports.not = exports.norm = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.mode = exports.mod = void 0;
|
8
8
|
var _configReadonly = require("./configReadonly.js");
|
9
9
|
var _factoriesNumber = require("../factoriesNumber.js");
|
10
10
|
/**
|
@@ -386,18 +386,7 @@ var smaller = exports.smaller = /* #__PURE__ */(0, _factoriesNumber.createSmalle
|
|
386
386
|
config: _configReadonly.config,
|
387
387
|
typed: typed
|
388
388
|
});
|
389
|
-
var
|
390
|
-
addScalar: addScalar,
|
391
|
-
combinations: combinations,
|
392
|
-
divideScalar: divideScalar,
|
393
|
-
factorial: factorial,
|
394
|
-
isInteger: isInteger,
|
395
|
-
isNegative: isNegative,
|
396
|
-
larger: larger,
|
397
|
-
multiplyScalar: multiplyScalar,
|
398
|
-
number: number,
|
399
|
-
pow: pow,
|
400
|
-
subtract: subtract,
|
389
|
+
var subtractScalar = exports.subtractScalar = /* #__PURE__ */(0, _factoriesNumber.createSubtractScalar)({
|
401
390
|
typed: typed
|
402
391
|
});
|
403
392
|
var _true = exports._true = /* #__PURE__ */(0, _factoriesNumber.createTrue)({});
|
@@ -519,6 +508,20 @@ var std = exports.std = /* #__PURE__ */(0, _factoriesNumber.createStd)({
|
|
519
508
|
typed: typed,
|
520
509
|
variance: variance
|
521
510
|
});
|
511
|
+
var stirlingS2 = exports.stirlingS2 = /* #__PURE__ */(0, _factoriesNumber.createStirlingS2)({
|
512
|
+
addScalar: addScalar,
|
513
|
+
combinations: combinations,
|
514
|
+
divideScalar: divideScalar,
|
515
|
+
factorial: factorial,
|
516
|
+
isInteger: isInteger,
|
517
|
+
isNegative: isNegative,
|
518
|
+
larger: larger,
|
519
|
+
multiplyScalar: multiplyScalar,
|
520
|
+
number: number,
|
521
|
+
pow: pow,
|
522
|
+
subtractScalar: subtractScalar,
|
523
|
+
typed: typed
|
524
|
+
});
|
522
525
|
var ceil = exports.ceil = /* #__PURE__ */(0, _factoriesNumber.createCeil)({
|
523
526
|
config: _configReadonly.config,
|
524
527
|
round: round,
|
@@ -160,6 +160,19 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
|
|
160
160
|
undefined: undefined
|
161
161
|
};
|
162
162
|
var NUMERIC_CONSTANTS = ['NaN', 'Infinity'];
|
163
|
+
var ESCAPE_CHARACTERS = {
|
164
|
+
'"': '"',
|
165
|
+
"'": "'",
|
166
|
+
'\\': '\\',
|
167
|
+
'/': '/',
|
168
|
+
b: '\b',
|
169
|
+
f: '\f',
|
170
|
+
n: '\n',
|
171
|
+
r: '\r',
|
172
|
+
t: '\t'
|
173
|
+
// note that \u is handled separately in parseStringToken()
|
174
|
+
};
|
175
|
+
|
163
176
|
function initialState() {
|
164
177
|
return {
|
165
178
|
extraNodes: {},
|
@@ -1234,7 +1247,7 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
|
|
1234
1247
|
node = parseAccessors(state, node);
|
1235
1248
|
return node;
|
1236
1249
|
}
|
1237
|
-
return
|
1250
|
+
return parseString(state);
|
1238
1251
|
}
|
1239
1252
|
|
1240
1253
|
/**
|
@@ -1320,61 +1333,14 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
|
|
1320
1333
|
}
|
1321
1334
|
|
1322
1335
|
/**
|
1323
|
-
* Parse a double
|
1336
|
+
* Parse a single or double quoted string.
|
1324
1337
|
* @return {Node} node
|
1325
1338
|
* @private
|
1326
1339
|
*/
|
1327
|
-
function
|
1340
|
+
function parseString(state) {
|
1328
1341
|
var node, str;
|
1329
|
-
if (state.token === '"') {
|
1330
|
-
str =
|
1331
|
-
|
1332
|
-
// create constant
|
1333
|
-
node = new ConstantNode(str);
|
1334
|
-
|
1335
|
-
// parse index parameters
|
1336
|
-
node = parseAccessors(state, node);
|
1337
|
-
return node;
|
1338
|
-
}
|
1339
|
-
return parseSingleQuotesString(state);
|
1340
|
-
}
|
1341
|
-
|
1342
|
-
/**
|
1343
|
-
* Parse a string surrounded by double quotes "..."
|
1344
|
-
* @return {string}
|
1345
|
-
*/
|
1346
|
-
function parseDoubleQuotesStringToken(state) {
|
1347
|
-
var str = '';
|
1348
|
-
while (currentCharacter(state) !== '' && currentCharacter(state) !== '"') {
|
1349
|
-
if (currentCharacter(state) === '\\') {
|
1350
|
-
// escape character, immediately process the next
|
1351
|
-
// character to prevent stopping at a next '\"'
|
1352
|
-
var cNext = nextCharacter(state);
|
1353
|
-
if (cNext !== "'") {
|
1354
|
-
str += currentCharacter(state);
|
1355
|
-
}
|
1356
|
-
next(state);
|
1357
|
-
}
|
1358
|
-
str += currentCharacter(state);
|
1359
|
-
next(state);
|
1360
|
-
}
|
1361
|
-
getToken(state);
|
1362
|
-
if (state.token !== '"') {
|
1363
|
-
throw createSyntaxError(state, 'End of string " expected');
|
1364
|
-
}
|
1365
|
-
getToken(state);
|
1366
|
-
return JSON.parse('"' + str + '"'); // unescape escaped characters
|
1367
|
-
}
|
1368
|
-
|
1369
|
-
/**
|
1370
|
-
* Parse a single quotes string.
|
1371
|
-
* @return {Node} node
|
1372
|
-
* @private
|
1373
|
-
*/
|
1374
|
-
function parseSingleQuotesString(state) {
|
1375
|
-
var node, str;
|
1376
|
-
if (state.token === '\'') {
|
1377
|
-
str = parseSingleQuotesStringToken(state);
|
1342
|
+
if (state.token === '"' || state.token === "'") {
|
1343
|
+
str = parseStringToken(state, state.token);
|
1378
1344
|
|
1379
1345
|
// create constant
|
1380
1346
|
node = new ConstantNode(str);
|
@@ -1387,30 +1353,47 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
|
|
1387
1353
|
}
|
1388
1354
|
|
1389
1355
|
/**
|
1390
|
-
* Parse a string surrounded by single quotes
|
1356
|
+
* Parse a string surrounded by single or double quotes
|
1357
|
+
* @param {Object} state
|
1358
|
+
* @param {"'" | "\""} quote
|
1391
1359
|
* @return {string}
|
1392
1360
|
*/
|
1393
|
-
function
|
1361
|
+
function parseStringToken(state, quote) {
|
1394
1362
|
var str = '';
|
1395
|
-
while (currentCharacter(state) !== '' && currentCharacter(state) !==
|
1363
|
+
while (currentCharacter(state) !== '' && currentCharacter(state) !== quote) {
|
1396
1364
|
if (currentCharacter(state) === '\\') {
|
1397
|
-
|
1398
|
-
|
1399
|
-
var
|
1400
|
-
if (
|
1401
|
-
|
1365
|
+
next(state);
|
1366
|
+
var _char = currentCharacter(state);
|
1367
|
+
var escapeChar = ESCAPE_CHARACTERS[_char];
|
1368
|
+
if (escapeChar !== undefined) {
|
1369
|
+
// an escaped control character like \" or \n
|
1370
|
+
str += escapeChar;
|
1371
|
+
state.index += 1;
|
1372
|
+
} else if (_char === 'u') {
|
1373
|
+
// escaped unicode character
|
1374
|
+
var unicode = state.expression.slice(state.index + 1, state.index + 5);
|
1375
|
+
if (/^[0-9A-Fa-f]{4}$/.test(unicode)) {
|
1376
|
+
// test whether the string holds four hexadecimal values
|
1377
|
+
str += String.fromCharCode(parseInt(unicode, 16));
|
1378
|
+
state.index += 5;
|
1379
|
+
} else {
|
1380
|
+
throw createSyntaxError(state, "Invalid unicode character \\u".concat(unicode));
|
1381
|
+
}
|
1382
|
+
} else {
|
1383
|
+
throw createSyntaxError(state, "Bad escape character \\".concat(_char));
|
1402
1384
|
}
|
1385
|
+
} else {
|
1386
|
+
// any regular character
|
1387
|
+
str += currentCharacter(state);
|
1403
1388
|
next(state);
|
1404
1389
|
}
|
1405
|
-
str += currentCharacter(state);
|
1406
|
-
next(state);
|
1407
1390
|
}
|
1408
1391
|
getToken(state);
|
1409
|
-
if (state.token !==
|
1410
|
-
throw createSyntaxError(state,
|
1392
|
+
if (state.token !== quote) {
|
1393
|
+
throw createSyntaxError(state, "End of string ".concat(quote, " expected"));
|
1411
1394
|
}
|
1412
1395
|
getToken(state);
|
1413
|
-
return
|
1396
|
+
return str;
|
1414
1397
|
}
|
1415
1398
|
|
1416
1399
|
/**
|
@@ -1505,10 +1488,8 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
|
|
1505
1488
|
getToken(state);
|
1506
1489
|
if (state.token !== '}') {
|
1507
1490
|
// parse key
|
1508
|
-
if (state.token === '"') {
|
1509
|
-
key =
|
1510
|
-
} else if (state.token === '\'') {
|
1511
|
-
key = parseSingleQuotesStringToken(state);
|
1491
|
+
if (state.token === '"' || state.token === "'") {
|
1492
|
+
key = parseStringToken(state, state.token);
|
1512
1493
|
} else if (state.tokenType === TOKENTYPE.SYMBOL || state.tokenType === TOKENTYPE.DELIMITER && state.token in NAMED_DELIMITERS) {
|
1513
1494
|
key = state.token;
|
1514
1495
|
getToken(state);
|
package/lib/cjs/factoriesAny.js
CHANGED
@@ -1929,6 +1929,12 @@ Object.defineProperty(exports, "createSubtract", {
|
|
1929
1929
|
return _subtract.createSubtract;
|
1930
1930
|
}
|
1931
1931
|
});
|
1932
|
+
Object.defineProperty(exports, "createSubtractScalar", {
|
1933
|
+
enumerable: true,
|
1934
|
+
get: function get() {
|
1935
|
+
return _subtractScalar.createSubtractScalar;
|
1936
|
+
}
|
1937
|
+
});
|
1932
1938
|
Object.defineProperty(exports, "createSum", {
|
1933
1939
|
enumerable: true,
|
1934
1940
|
get: function get() {
|
@@ -2174,6 +2180,7 @@ var _unaryPlus = require("./function/arithmetic/unaryPlus.js");
|
|
2174
2180
|
var _abs = require("./function/arithmetic/abs.js");
|
2175
2181
|
var _apply = require("./function/matrix/apply.js");
|
2176
2182
|
var _addScalar = require("./function/arithmetic/addScalar.js");
|
2183
|
+
var _subtractScalar = require("./function/arithmetic/subtractScalar.js");
|
2177
2184
|
var _cbrt = require("./function/arithmetic/cbrt.js");
|
2178
2185
|
var _ceil = require("./function/arithmetic/ceil.js");
|
2179
2186
|
var _cube = require("./function/arithmetic/cube.js");
|
@@ -702,7 +702,7 @@ Object.defineProperty(exports, "createString", {
|
|
702
702
|
return _string.createString;
|
703
703
|
}
|
704
704
|
});
|
705
|
-
exports.createSubtract = exports.createSubsetTransform = exports.createSubset = void 0;
|
705
|
+
exports.createSubtractScalar = exports.createSubtract = exports.createSubsetTransform = exports.createSubset = void 0;
|
706
706
|
Object.defineProperty(exports, "createSum", {
|
707
707
|
enumerable: true,
|
708
708
|
get: function get() {
|
@@ -916,6 +916,7 @@ var createUnaryMinus = exports.createUnaryMinus = /* #__PURE__ */createNumberFac
|
|
916
916
|
var createUnaryPlus = exports.createUnaryPlus = /* #__PURE__ */createNumberFactory('unaryPlus', _index.unaryPlusNumber);
|
917
917
|
var createAbs = exports.createAbs = /* #__PURE__ */createNumberFactory('abs', _index.absNumber);
|
918
918
|
var createAddScalar = exports.createAddScalar = /* #__PURE__ */createNumberFactory('addScalar', _index.addNumber);
|
919
|
+
var createSubtractScalar = exports.createSubtractScalar = /* #__PURE__ */createNumberFactory('subtractScalar', _index.subtractNumber);
|
919
920
|
var createCbrt = exports.createCbrt = /* #__PURE__ */createNumberFactory('cbrt', _index.cbrtNumber);
|
920
921
|
var createCube = exports.createCube = /* #__PURE__ */createNumberFactory('cube', _index.cubeNumber);
|
921
922
|
var createExp = exports.createExp = /* #__PURE__ */createNumberFactory('exp', _index.expNumber);
|
@@ -7,7 +7,7 @@ exports.createLup = void 0;
|
|
7
7
|
var _object = require("../../../utils/object.js");
|
8
8
|
var _factory = require("../../../utils/factory.js");
|
9
9
|
var name = 'lup';
|
10
|
-
var dependencies = ['typed', 'matrix', 'abs', 'addScalar', 'divideScalar', 'multiplyScalar', '
|
10
|
+
var dependencies = ['typed', 'matrix', 'abs', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'larger', 'equalScalar', 'unaryMinus', 'DenseMatrix', 'SparseMatrix', 'Spa'];
|
11
11
|
var createLup = exports.createLup = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
|
12
12
|
var typed = _ref.typed,
|
13
13
|
matrix = _ref.matrix,
|
@@ -15,7 +15,7 @@ var createLup = exports.createLup = /* #__PURE__ */(0, _factory.factory)(name, d
|
|
15
15
|
addScalar = _ref.addScalar,
|
16
16
|
divideScalar = _ref.divideScalar,
|
17
17
|
multiplyScalar = _ref.multiplyScalar,
|
18
|
-
|
18
|
+
subtractScalar = _ref.subtractScalar,
|
19
19
|
larger = _ref.larger,
|
20
20
|
equalScalar = _ref.equalScalar,
|
21
21
|
unaryMinus = _ref.unaryMinus,
|
@@ -104,7 +104,7 @@ var createLup = exports.createLup = /* #__PURE__ */(0, _factory.factory)(name, d
|
|
104
104
|
// s = l[i, k] - data[k, j]
|
105
105
|
s = addScalar(s, multiplyScalar(data[i][k], data[k][j]));
|
106
106
|
}
|
107
|
-
data[i][j] =
|
107
|
+
data[i][j] = subtractScalar(data[i][j], s);
|
108
108
|
}
|
109
109
|
}
|
110
110
|
// row with larger value in cvector, row >= j
|
@@ -8,7 +8,7 @@ exports.createQr = void 0;
|
|
8
8
|
var _extends2 = _interopRequireDefault(require("@babel/runtime/helpers/extends"));
|
9
9
|
var _factory = require("../../../utils/factory.js");
|
10
10
|
var name = 'qr';
|
11
|
-
var dependencies = ['typed', 'matrix', 'zeros', 'identity', 'isZero', 'equal', 'sign', 'sqrt', 'conj', 'unaryMinus', 'addScalar', 'divideScalar', 'multiplyScalar', '
|
11
|
+
var dependencies = ['typed', 'matrix', 'zeros', 'identity', 'isZero', 'equal', 'sign', 'sqrt', 'conj', 'unaryMinus', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'complex'];
|
12
12
|
var createQr = exports.createQr = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
|
13
13
|
var typed = _ref.typed,
|
14
14
|
matrix = _ref.matrix,
|
@@ -23,7 +23,7 @@ var createQr = exports.createQr = /* #__PURE__ */(0, _factory.factory)(name, dep
|
|
23
23
|
addScalar = _ref.addScalar,
|
24
24
|
divideScalar = _ref.divideScalar,
|
25
25
|
multiplyScalar = _ref.multiplyScalar,
|
26
|
-
|
26
|
+
subtractScalar = _ref.subtractScalar,
|
27
27
|
complex = _ref.complex;
|
28
28
|
/**
|
29
29
|
* Calculate the Matrix QR decomposition. Matrix `A` is decomposed in
|
@@ -137,7 +137,7 @@ var createQr = exports.createQr = /* #__PURE__ */(0, _factory.factory)(name, dep
|
|
137
137
|
var alpha = multiplyScalar(sgn, sqrt(alphaSquared));
|
138
138
|
if (!isZero(alpha)) {
|
139
139
|
// first element in vector u
|
140
|
-
var u1 =
|
140
|
+
var u1 = subtractScalar(pivot, alpha);
|
141
141
|
|
142
142
|
// w = v * u1 / |u| (only elements k to (rows-1) are used)
|
143
143
|
w[k] = 1;
|
@@ -172,7 +172,7 @@ var createQr = exports.createQr = /* #__PURE__ */(0, _factory.factory)(name, dep
|
|
172
172
|
// calculate the jth element of [tau * transpose(w) * R]
|
173
173
|
s = multiplyScalar(s, tau);
|
174
174
|
for (i = k; i < rows; i++) {
|
175
|
-
Rdata[i][j] = multiplyScalar(
|
175
|
+
Rdata[i][j] = multiplyScalar(subtractScalar(Rdata[i][j], multiplyScalar(w[i], s)), conjSgn);
|
176
176
|
}
|
177
177
|
}
|
178
178
|
/*
|
@@ -193,7 +193,7 @@ var createQr = exports.createQr = /* #__PURE__ */(0, _factory.factory)(name, dep
|
|
193
193
|
// calculate the ith element of [tau * Q * w]
|
194
194
|
s = multiplyScalar(s, tau);
|
195
195
|
for (j = k; j < rows; ++j) {
|
196
|
-
Qdata[i][j] = divideScalar(
|
196
|
+
Qdata[i][j] = divideScalar(subtractScalar(Qdata[i][j], multiplyScalar(s, conj(w[j]))), conjSgn);
|
197
197
|
}
|
198
198
|
}
|
199
199
|
}
|
@@ -7,13 +7,13 @@ exports.createLsolve = void 0;
|
|
7
7
|
var _factory = require("../../../utils/factory.js");
|
8
8
|
var _solveValidation = require("./utils/solveValidation.js");
|
9
9
|
var name = 'lsolve';
|
10
|
-
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', '
|
10
|
+
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
|
11
11
|
var createLsolve = exports.createLsolve = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
|
12
12
|
var typed = _ref.typed,
|
13
13
|
matrix = _ref.matrix,
|
14
14
|
divideScalar = _ref.divideScalar,
|
15
15
|
multiplyScalar = _ref.multiplyScalar,
|
16
|
-
|
16
|
+
subtractScalar = _ref.subtractScalar,
|
17
17
|
equalScalar = _ref.equalScalar,
|
18
18
|
DenseMatrix = _ref.DenseMatrix;
|
19
19
|
var solveValidation = (0, _solveValidation.createSolveValidation)({
|
@@ -83,7 +83,7 @@ var createLsolve = exports.createLsolve = /* #__PURE__ */(0, _factory.factory)(n
|
|
83
83
|
|
84
84
|
// loop rows
|
85
85
|
for (var i = j + 1; i < rows; i++) {
|
86
|
-
bdata[i] = [
|
86
|
+
bdata[i] = [subtractScalar(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
|
87
87
|
}
|
88
88
|
} else {
|
89
89
|
// degenerate row, we can choose any value
|
@@ -145,7 +145,7 @@ var createLsolve = exports.createLsolve = /* #__PURE__ */(0, _factory.factory)(n
|
|
145
145
|
var xj = divideScalar(bj, vjj);
|
146
146
|
for (var _k = 0, l = jIndices.length; _k < l; _k++) {
|
147
147
|
var _i = jIndices[_k];
|
148
|
-
bdata[_i] = [
|
148
|
+
bdata[_i] = [subtractScalar(bdata[_i][0] || 0, multiplyScalar(xj, jValues[_k]))];
|
149
149
|
}
|
150
150
|
x[j] = [xj];
|
151
151
|
} else {
|
@@ -9,13 +9,13 @@ var _toConsumableArray2 = _interopRequireDefault(require("@babel/runtime/helpers
|
|
9
9
|
var _factory = require("../../../utils/factory.js");
|
10
10
|
var _solveValidation = require("./utils/solveValidation.js");
|
11
11
|
var name = 'lsolveAll';
|
12
|
-
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', '
|
12
|
+
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
|
13
13
|
var createLsolveAll = exports.createLsolveAll = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
|
14
14
|
var typed = _ref.typed,
|
15
15
|
matrix = _ref.matrix,
|
16
16
|
divideScalar = _ref.divideScalar,
|
17
17
|
multiplyScalar = _ref.multiplyScalar,
|
18
|
-
|
18
|
+
subtractScalar = _ref.subtractScalar,
|
19
19
|
equalScalar = _ref.equalScalar,
|
20
20
|
DenseMatrix = _ref.DenseMatrix;
|
21
21
|
var solveValidation = (0, _solveValidation.createSolveValidation)({
|
@@ -86,7 +86,7 @@ var createLsolveAll = exports.createLsolveAll = /* #__PURE__ */(0, _factory.fact
|
|
86
86
|
b[i] = divideScalar(b[i], M[i][i]);
|
87
87
|
for (var j = i + 1; j < columns; j++) {
|
88
88
|
// b[j] -= b[i] * M[j,i]
|
89
|
-
b[j] =
|
89
|
+
b[j] = subtractScalar(b[j], multiplyScalar(b[i], M[j][i]));
|
90
90
|
}
|
91
91
|
} else if (!equalScalar(b[i], 0)) {
|
92
92
|
// singular row, nonzero RHS
|
@@ -106,7 +106,7 @@ var createLsolveAll = exports.createLsolveAll = /* #__PURE__ */(0, _factory.fact
|
|
106
106
|
var bNew = (0, _toConsumableArray2["default"])(b);
|
107
107
|
bNew[i] = 1;
|
108
108
|
for (var _j = i + 1; _j < columns; _j++) {
|
109
|
-
bNew[_j] =
|
109
|
+
bNew[_j] = subtractScalar(bNew[_j], M[_j][i]);
|
110
110
|
}
|
111
111
|
B.push(bNew);
|
112
112
|
}
|
@@ -167,7 +167,7 @@ var createLsolveAll = exports.createLsolveAll = /* #__PURE__ */(0, _factory.fact
|
|
167
167
|
b[i] = divideScalar(b[i], Mii);
|
168
168
|
for (var _j2 = 0, _lastIndex = iIndices.length; _j2 < _lastIndex; _j2++) {
|
169
169
|
var _J = iIndices[_j2];
|
170
|
-
b[_J] =
|
170
|
+
b[_J] = subtractScalar(b[_J], multiplyScalar(b[i], iValues[_j2]));
|
171
171
|
}
|
172
172
|
} else if (!equalScalar(b[i], 0)) {
|
173
173
|
// singular row, nonzero RHS
|
@@ -188,7 +188,7 @@ var createLsolveAll = exports.createLsolveAll = /* #__PURE__ */(0, _factory.fact
|
|
188
188
|
bNew[i] = 1;
|
189
189
|
for (var _j3 = 0, _lastIndex2 = iIndices.length; _j3 < _lastIndex2; _j3++) {
|
190
190
|
var _J2 = iIndices[_j3];
|
191
|
-
bNew[_J2] =
|
191
|
+
bNew[_J2] = subtractScalar(bNew[_J2], iValues[_j3]);
|
192
192
|
}
|
193
193
|
B.push(bNew);
|
194
194
|
}
|
@@ -7,13 +7,13 @@ exports.createUsolve = void 0;
|
|
7
7
|
var _factory = require("../../../utils/factory.js");
|
8
8
|
var _solveValidation = require("./utils/solveValidation.js");
|
9
9
|
var name = 'usolve';
|
10
|
-
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', '
|
10
|
+
var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
|
11
11
|
var createUsolve = exports.createUsolve = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
|
12
12
|
var typed = _ref.typed,
|
13
13
|
matrix = _ref.matrix,
|
14
14
|
divideScalar = _ref.divideScalar,
|
15
15
|
multiplyScalar = _ref.multiplyScalar,
|
16
|
-
|
16
|
+
subtractScalar = _ref.subtractScalar,
|
17
17
|
equalScalar = _ref.equalScalar,
|
18
18
|
DenseMatrix = _ref.DenseMatrix;
|
19
19
|
var solveValidation = (0, _solveValidation.createSolveValidation)({
|
@@ -85,7 +85,7 @@ var createUsolve = exports.createUsolve = /* #__PURE__ */(0, _factory.factory)(n
|
|
85
85
|
// loop rows
|
86
86
|
for (var i = j - 1; i >= 0; i--) {
|
87
87
|
// update copy of b
|
88
|
-
bdata[i] = [
|
88
|
+
bdata[i] = [subtractScalar(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
|
89
89
|
}
|
90
90
|
} else {
|
91
91
|
// zero value at j
|
@@ -149,7 +149,7 @@ var createUsolve = exports.createUsolve = /* #__PURE__ */(0, _factory.factory)(n
|
|
149
149
|
var xj = divideScalar(bj, vjj);
|
150
150
|
for (var _k = 0, _lastIndex = jIndices.length; _k < _lastIndex; _k++) {
|
151
151
|
var _i = jIndices[_k];
|
152
|
-
bdata[_i] = [
|
152
|
+
bdata[_i] = [subtractScalar(bdata[_i][0], multiplyScalar(xj, jValues[_k]))];
|
153
153
|
}
|
154
154
|
x[j] = [xj];
|
155
155
|
} else {
|