default random seed
fixed power function for R semantics, now also using libmath (via JNI) as it is faster