Package jline.util
Class ComplexOps
java.lang.Object
jline.util.ComplexOps
The few complex operations the transform layer needs.
WHY THIS EXISTS. A Laplace-Stieltjes transform is evaluated OFF the real axis
by everything that inverts it or locates its roots: the Abate-Whitt Euler sum
walks the line Re(s) = A/(2t), and the matrix transform int exp(Ut) dF(t) is
read off the spectrum of U, which is complex in general.
ComplexMatrix
carries complex data but offers no solve, and pulling in a full complex linear
algebra stack for an n x n system of the size a phase-type law has would be
out of proportion. Gaussian elimination with partial pivoting is enough and is
what this class provides.-
Method Summary
Modifier and TypeMethodDescriptionstatic org.apache.commons.math3.complex.Complex[]solve(org.apache.commons.math3.complex.Complex[][] A, org.apache.commons.math3.complex.Complex[] b) Solves A x = b by Gaussian elimination with partial pivoting.
-
Method Details
-
solve
public static org.apache.commons.math3.complex.Complex[] solve(org.apache.commons.math3.complex.Complex[][] A, org.apache.commons.math3.complex.Complex[] b) Solves A x = b by Gaussian elimination with partial pivoting.- Parameters:
A- n x n coefficient matrix, consumed by value (a copy is taken)b- right-hand side of length n- Returns:
- the solution x
-