/Users/garethloy/Musimathics/Musimat1.2/include/MathFuns.h

Go to the documentation of this file.

/* $Revision: 1.4 $ $Date: 2006/09/07 08:38:33 $ $Author: dgl $ $Name:  $ $Id: MathFuns.h,v 1.4 2006/09/07 08:38:33 dgl Exp $ */
#ifndef MATHFUNS_H
#define MATHFUNS_H

// The Musimat Tutorial © 2006 Gareth Loy
// Derived from Chapter 9 and Appendix B of "Musimathics Vol. 1" © 2006 Gareth Loy 
// and published exclusively by The MIT Press.
// This program is released WITHOUT ANY WARRANTY; without even the implied 
// warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. 
// For information on usage and redistribution, and for a DISCLAIMER OF ALL
// WARRANTIES, see the file, "LICENSE.txt," in this distribution.
// "Musimathics" is available here:     http://mitpress.mit.edu/catalog/item/default.asp?ttype=2&tid=10916
// Gareth Loy's Musimathics website:    http://www.musimathics.com/
// The Musimat website:                 http://www.musimat.com/
// This program is released under the terms of the GNU General Public License
// available here:                      http://www.gnu.org/licenses/gpl.txt

#include "Musimat.h"
#include <math.h>


Const Real Pi = 4.0 * atan( 1.0 );


template<class Type> inline Type Abs( Type x ) { Return x < 0 ? -x : x; }

template<class Type> inline Type Mod(Type j, Type k) {
        While ( j >= k ) { j = j - k; }
        While ( j <= -k ) { j = j + k; }
        Return( j );
}

template<class Type> inline Type PosMod( Type j, Type k )
{
        While ( j >= k ) {
                j = j - k; 
        }
        While ( j < 0 ) {
                j = j + k; 
        }
        Return( j );
}

inline Real Floor( Real x ) { Return floor( x ); }

inline Real Ceiling( Real x ) { Return ceil( x ); }

inline Real Sqrt( Real x ) { Return sqrt( x ); }

inline Real Pow( Real x, Real y ) { Return pow( x, y ); }

template<class Type> inline Type Log( Type x ) { Return log( x ); }

template<class Type> inline Type Log10( Type x ) { Return log( x ) / log (10.0); }

template<class Type> inline Type Atan( Type x ) { Return atan( x ); }

template<class Type> inline Type Atan2( Type x, Type y ) { Return atan2( x, y ); }

template<class Type> inline Type Sin( Type x ) { Return sin( x ); }

template<class Type> inline Type Cos( Type x ) { Return cos( x ); }

inline Real unitInterp( Real f, Integer L, Integer U ) {
        Return( f * ( U - L ) + L );
}

inline Real linearInterpolate( Integer x, Integer xMin, Integer xMax, Integer yMin, Integer yMax ) {
        Real a = Real( x - xMin) / Real(xMax - xMin); // Real division
        Real b = Real( yMax - yMin );
        Return( a * b + yMin );
}

inline Real Round( Real x ) { 
        Return( Floor( x + 0.5 ) );
}

Integer RealToRational( Real f, Integer Reference a, Integer Reference b);



template <class Type> Type Plus(Type i, Type j) { Return i + j; }

template <class Type> Type Minus(Type i, Type j) { Return i - j; }

template <class Type> Type Divide(Type i, Type j) { Return i / j; }

template <class Type> Type Multiply(Type i, Type j) { Return i * j; }

#endif // MATHFUNS_H