code: public BandLU factor(BandMatrix A, boolean inplace) {
          if (inplace)
              return factor(A);
          else
              return factor(new BandMatrix(A, kl, kl + ku));
      }
docs: Creates an LU decomposition of the given matrix
  
  @param A
  Matrix to decompose. If the decomposition is in-place, its
  number of superdiagonals must equal <code>kl+ku</code>
  @param inplace
  Wheter or not the decomposition should overwrite the passed
  matrix
  @return The current decomposition