Ejemplo n.º 1
0
  /** ************************************************************************* LAPACK */
  public static DoubleMatrix gesv(DoubleMatrix a, int[] ipiv, DoubleMatrix b) {
    int info = NativeBlas.dgesv(a.rows, b.columns, a.data, 0, a.rows, ipiv, 0, b.data, 0, b.rows);
    checkInfo("DGESV", info);

    if (info > 0)
      throw new LapackException(
          "DGESV", "Linear equation cannot be solved because the matrix was singular.");

    return b;
  }