Add import of NumPy and correct a missing np. for the ceil function