Example usage for org.apache.commons.math.linear BigMatrixImpl BigMatrixImpl

Introduction

In this page you can find the example usage for org.apache.commons.math.linear BigMatrixImpl BigMatrixImpl.

Prototype

public BigMatrixImpl(BigDecimal[] v) 

Source Link

Document

Create a new (column) BigMatrix using v as the data for the unique column of the v.length x 1 matrix created.

Usage

From source file:rmc.mat.RmcNewtonSolver.java

public void iterateUntilMSVlessThan(double maxMSV, long maxIterations) throws Exception {
    while (currentMSV > maxMSV && currentIterations < maxIterations) {
        System.out.println(currentMSV);
        currentMSV = 0;/*from w w  w. j  a v  a2s .c  o m*/
        currentIterations++;
        BigDecimal jacobianMatrix[][] = new BigDecimal[inputVariables.size()][inputVariables.size()];
        for (int row = 0; row < inputVariables.size(); row++) {
            for (int col = 0; col < inputVariables.size(); col++) {
                jacobianMatrix[row][col] = jacobianEquationMatrix[row][col].calculate();
            }
        }
        BigMatrix invertedJacobianMatrix = new BigMatrixImpl(jacobianMatrix).inverse();
        for (int row = 0; row < inputVariables.size(); row++) {
            BigDecimal diffValue = new BigDecimal(0);
            for (int col = 0; col < inputVariables.size(); col++) {
                BigDecimal functionToOptimizeValue = functionToOptimize[col].calculate();
                diffValue = diffValue
                        .add(invertedJacobianMatrix.getEntry(row, col).multiply(functionToOptimizeValue));
            }
            resultVariables.set(row,
                    inputVariables.getVariableValueFromId(row).subtract(diffValue).stripTrailingZeros());
        }
        for (RmcEquation rmcEquation : rmcEquations) {
            currentMSV = Math.max(currentMSV,
                    resultVariables.get(rmcEquation.variableId)
                            .subtract(inputVariables.getVariableValueFromId(rmcEquation.variableId)).abs()
                            .doubleValue());
            inputVariables.setVariableValue(rmcEquation.variableName,
                    resultVariables.get(rmcEquation.variableId));
        }
    }
}

---