util.h

/*************************************************************************************
 * MechSys - A C++ library to simulate (Continuum) Mechanical Systems                *
 * Copyright (C) 2005 Dorival de Moraes Pedroso <dorival.pedroso at gmail.com>       *
 * Copyright (C) 2005 Raul Dario Durand Farfan  <raul.durand at gmail.com>           *
 *                                                                                   *
 * This file is part of MechSys.                                                     *
 *                                                                                   *
 * MechSys is free software; you can redistribute it and/or modify it under the      *
 * terms of the GNU General Public License as published by the Free Software         *
 * Foundation; either version 2 of the License, or (at your option) any later        *
 * version.                                                                          *
 *                                                                                   *
 * MechSys is distributed in the hope that it will be useful, but WITHOUT ANY        *
 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A   *
 * PARTICULAR PURPOSE. See the GNU General Public License for more details.          *
 *                                                                                   *
 * You should have received a copy of the GNU General Public License along with      *
 * MechSys; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, *
 * Fifth Floor, Boston, MA 02110-1301, USA                                           *
 *************************************************************************************/

#ifndef MECHSYS_SORT_H
#define MECHSYS_SORT_H

#include <cmath>

#ifdef HAVE_CONFIG_H
  #include "config.h"
#else
  #ifndef REAL
    #define REAL double
  #endif
#endif

#include "util/string.h"

namespace Util
{

inline REAL Pi         ()                    { return 3.14159265358979;             }
inline REAL ToRad      (REAL deg_angle)      { return 0.0174532925199433*deg_angle; }
inline REAL ToDeg      (REAL rad_angle)      { return 57.2957795130823*rad_angle;   }
inline REAL Sgn        (REAL Val)            { return (Val>=0.0 ? +1.0 : -1.0);     }
inline REAL Signal     (REAL Val, REAL Zero) { return (fabs(Val)<=Zero ? 0.0 : Sgn(Val)); }
inline REAL Acos       (REAL Val)            { return (Val>=1.0 ?  0.0 : (Val<=-1.0 ? Pi() : acos(Val)) ); }
inline bool Str2Bool   (String const & Str)  { if (Str==String("true") || Str==String("TRUE") || Str==String("True")) return true; else return false; }
inline bool IsNanOrInf (REAL Val)            { int r=std::fpclassify(Val); if (r==FP_NAN) return true; if (r==FP_INFINITE) return true; return false; }
inline REAL Sq2        ()                    { return 1.41421356237310; }
inline REAL Sq3        ()                    { return 1.73205080756888; }
inline REAL Sq6        ()                    { return 2.44948974278318; }
inline REAL Sign       (REAL a, REAL b)      { return (b>=0.0 ? fabs(a) : -fabs(a)); }
template <typename Type>
inline Type Min        (Type a, Type b)      { return (a<b ? a : b); }
template <typename Type>
inline Type Max        (Type a, Type b)      { return (a>b ? a : b); }

inline void Swap(REAL & x, REAL & y)
{ // {{{
    REAL temp = x;
    x = y;
    y = temp;
} // }}}

inline void Sort(REAL A[], int Size) // ascending sort
{ // {{{
    for (int i=0; i<Size-1; ++i)
    {
        int min_index = i;

        // Find the index of the minimum element
        for (int j=i+1; j<Size; ++j)
            if (A[j] < A[min_index])
                min_index = j;

        // Swap if i-th element not already smallest
        if (min_index > i)
            Util::Swap(A[i], A[min_index]);
    }
} // }}}

inline void FindBestSquare (int Size, int & nRow, int & nCol)
{ // {{{
    nRow = -1;  // not found
    nCol = -1;  // not found
    for (int x=1; x<=Size; ++x)
    {
        if ((x*x)>=Size)
        {
            if ((x*x)==Size)
            {
                nRow = x;
                nCol = x;
                return;
            }
            else
            {
                for (int y=x; y>=1; --y)
                {
                    if ((x*y)==Size)
                    {
                        nRow = x;
                        nCol = y;
                        return;
                    }
                }
            }
        }
    }
} // }}}

}; // namespace Util

#endif // MECHSYS_SORT_H

// vim:fdm=marker

---