switched all scalar und performance uncritical floating point operations from fptype to double precision (this was long overdue and solving the Sa=f system for SR was actually problematic in single precision)