code: @Override
      public void spotrf(byte uplo, int N, INDArray A, INDArray INFO) {
          int status = LAPACKE_spotrf(getColumnOrder(A), uplo, N, 
                          (FloatPointer)A.data().addressPointer(), getLda(A) );
          if( status != 0 ) {
              throw new BlasException( "Failed to execute spotrf", status ) ;
          }
          if( uplo == 'U' ) {			
  			INDArrayIndex ix[] = new INDArrayIndex[ 2 ] ;
  			for( int i=1 ; i<Math.min( A.rows(), A.columns() ) ; i++ ) {
  				ix[0] = NDArrayIndex.point( i ) ;
  				ix[1] = NDArrayIndex.interval( 0, i ) ;				
  				A.put(ix, 0) ;
  			}            
          } else {
              INDArrayIndex ix[] = new INDArrayIndex[ 2 ] ;
              for( int i=0 ; i<Math.min( A.rows(), A.columns()-1 ) ; i++ ) {
                  ix[0] = NDArrayIndex.point( i ) ;
                  ix[1] = NDArrayIndex.interval( i+1, A.columns() ) ;
                  A.put(ix, 0) ;
              }        
          }
      }
docs: CHOLESKY DECOMP