CAMLprim value c_arraydensematrix_potrf(value va) { CAMLparam1(va); struct caml_ba_array *ba = ARRAY2_DATA(va); intnat m = ba->dim[1]; #if SUNDIALS_ML_SAFE == 1 intnat n = ba->dim[0]; if (m != n) caml_invalid_argument("ArrayDenseMatrix.potrf: matrix not square"); #endif densePOTRF(ARRAY2_ACOLS(va), m); CAMLreturn (Val_unit); }
long int DensePOTRF(DlsMat A) { return(densePOTRF(A->cols, A->M)); }