double droundup_lwork(int LWORK) { var result = LWORK.toDouble(); if (result.toInt() < LWORK) { // Force round up of LWORK result *= (1.0 + epsilon(0.0)); } return result; }