Function MultTriangularUR

basiclinalg/triangular.cpp:587–598

  }

  void MultTriangularUR (SliceMatrix<double> X, BareSliceMatrix<double> U)  
  {
    // MultTriangularUR3 (X, L);
    // return;
    
    size_t i = 0;
    constexpr size_t bw = 128;
    for ( ; i+bw <= X.Height(); i += bw)
      MultTriangularUR2 (X.Rows(i,i+bw), U);
    if (i < X.Height())
      MultTriangularUR2 (X.Rows(i,X.Height()), U);      
  }