You can do operator overloading in fortran to add a new number type, don't remember the details but I wrote a quaternion library once.