mathjs 11.11.2 → 12.0.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/HISTORY.md +48 -0
- package/README.md +1 -1
- 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/dependenciesEigs.generated.js +4 -0
- 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 +43 -37
- package/lib/cjs/entry/pureFunctionsNumber.generated.js +16 -13
- package/lib/cjs/expression/embeddedDocs/function/matrix/eigs.js +2 -2
- package/lib/cjs/expression/node/AssignmentNode.js +1 -1
- package/lib/cjs/expression/node/FunctionAssignmentNode.js +1 -1
- 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/function/matrix/eigs/complexEigs.js +73 -68
- package/lib/cjs/function/matrix/eigs/{realSymetric.js → realSymmetric.js} +57 -51
- package/lib/cjs/function/matrix/eigs.js +118 -45
- package/lib/cjs/header.js +2 -2
- package/lib/cjs/type/unit/Unit.js +5 -5
- package/lib/cjs/utils/number.js +1 -1
- 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/dependenciesEigs.generated.js +4 -0
- 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 +42 -37
- package/lib/esm/entry/pureFunctionsNumber.generated.js +16 -13
- package/lib/esm/expression/embeddedDocs/function/matrix/eigs.js +2 -2
- package/lib/esm/expression/node/AssignmentNode.js +1 -1
- package/lib/esm/expression/node/FunctionAssignmentNode.js +1 -1
- 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/function/matrix/eigs/complexEigs.js +73 -68
- package/lib/esm/function/matrix/eigs/{realSymetric.js → realSymmetric.js} +55 -51
- package/lib/esm/function/matrix/eigs.js +119 -47
- package/lib/esm/type/unit/Unit.js +5 -5
- package/lib/esm/utils/number.js +1 -1
- 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 +14 -14
- package/types/EXPLANATION.md +54 -0
- package/types/index.d.ts +6797 -6446
|
@@ -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)({
|
|
@@ -1586,7 +1590,9 @@ var eigs = exports.eigs = /* #__PURE__ */(0, _factoriesAny.createEigs)({
|
|
|
1586
1590
|
number: number,
|
|
1587
1591
|
qr: qr,
|
|
1588
1592
|
re: re,
|
|
1593
|
+
reshape: reshape,
|
|
1589
1594
|
sin: sin,
|
|
1595
|
+
size: size,
|
|
1590
1596
|
smaller: smaller,
|
|
1591
1597
|
sqrt: sqrt,
|
|
1592
1598
|
subtract: subtract,
|
|
@@ -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,
|
|
@@ -8,7 +8,7 @@ var eigsDocs = exports.eigsDocs = {
|
|
|
8
8
|
name: 'eigs',
|
|
9
9
|
category: 'Matrix',
|
|
10
10
|
syntax: ['eigs(x)'],
|
|
11
|
-
description: 'Calculate the eigenvalues and eigenvectors of a
|
|
12
|
-
examples: ['eigs([[5, 2.3], [2.3, 1]])'],
|
|
11
|
+
description: 'Calculate the eigenvalues and optionally eigenvectors of a square matrix',
|
|
12
|
+
examples: ['eigs([[5, 2.3], [2.3, 1]])', 'eigs([[1, 2, 3], [4, 5, 6], [7, 8, 9]], { precision: 1e-6, eigenvectors: false }'],
|
|
13
13
|
seealso: ['inv']
|
|
14
14
|
};
|
|
@@ -331,7 +331,7 @@ var createAssignmentNode = exports.createAssignmentNode = /* #__PURE__ */(0, _fa
|
|
|
331
331
|
if (needParenthesis(this, options && options.parenthesis, options && options.implicit)) {
|
|
332
332
|
value = "\\left(".concat(value, "\\right)");
|
|
333
333
|
}
|
|
334
|
-
return object + index + '
|
|
334
|
+
return object + index + '=' + value;
|
|
335
335
|
}
|
|
336
336
|
}], [{
|
|
337
337
|
key: "fromJSON",
|
|
@@ -270,7 +270,7 @@ var createFunctionAssignmentNode = exports.createFunctionAssignmentNode = /* #__
|
|
|
270
270
|
if (needParenthesis(this, parenthesis, options && options.implicit)) {
|
|
271
271
|
expr = "\\left(".concat(expr, "\\right)");
|
|
272
272
|
}
|
|
273
|
-
return '\\mathrm{' + this.name + '}\\left(' + this.params.map(_latex.toSymbol).join(',') + '\\right)
|
|
273
|
+
return '\\mathrm{' + this.name + '}\\left(' + this.params.map(_latex.toSymbol).join(',') + '\\right)=' + expr;
|
|
274
274
|
}
|
|
275
275
|
}], [{
|
|
276
276
|
key: "fromJSON",
|
|
@@ -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
|
}
|