OSMathUtil.h

Go to the documentation of this file.

/* $Id: OSMathUtil.h 5284 2017-12-08 13:52:50Z stefan $ */
#ifndef MATHUTIL_H
#define MATHUTIL_H

#include "OSParameters.h"
#include "OSdtoa.h"
#include "OSErrorClass.h"
#include "OSGeneral.h"


#include <string>

#ifdef __cplusplus
extern std::string os_dtoa_format(double  x);
#endif

#ifdef __cplusplus
extern "C" {
#endif

    double os_strtod_wrap(const char *str,   char **strEnd);

#ifdef __cplusplus
}
#endif

class MathUtil
{
public:

    MathUtil();

    ~MathUtil();

    /*public static double round (double x, int precision){
        double mask = 0.0, y, result;
        try{
                mask = Math.pow (10.0, -(double)precision);
        }
        catch(ArithmeticException e){
                return x;
        }
        y  = mod(x, mask);
        result  = x - y;
        if (y / mask >= 0.5) result += mask;
        return result;
    }//round
    */
    /*public static double mod (double x, double y){
        return  x - Math.floor(x / y) * y;
    }//mod
    */

    static SparseMatrix* convertLinearConstraintCoefficientMatrixToTheOtherMajor(
        bool isColumnMajor, int startSize, int valueSize, int* start, int* index,
        double* value, int dimension);

    std::string format_os_dtoa( double x);


    //double os_strtod_wrap(const char *str) throw(ErrorClass);

    //double os_strtod_wrap(const char *str,  const char *strEnd) throw(ErrorClass);



};//class MathUtil


inline void getMultIncr(int* i, int *mult, int *incr, int size, int defaultIncr)
{
    int k;

    *mult = 1;
    *incr = defaultIncr;

    if (size == 1) return;

    for (k=1; (k < size) && (i[k] - i[k-1] == defaultIncr); k++)
    {
        (*mult)++;
    }
    if (*mult > 1 || size == 2) return;

    *incr = i[1] - i[0];
    if (i[2] - i[1] != *incr) return;

    *mult = 3;
    for (k=3; (k < size) && (i[k] - i[k-1] == *incr); k++)
    {
        (*mult)++;
    }
    return;
}


inline void getMultIncr(double* a, int *mult, double *incr, int size)
{
    double mark;
    int k;

    *mult = 1;
    *incr = 0;

    if (size == 1) return;

    mark = a[0];
    for (k=1; (k < size) && (a[k] == mark); k++)
    {
        (*mult)++;
    }
    if (*mult > 1 || size == 2) return;

    *incr = a[1] - a[0];
    if (a[2] - a[1] != *incr) return;

    *mult = 3;
    for (k=3; (k < size) && (a[k] - a[k-1] == *incr); k++)
    {
        (*mult)++;
    }
    return;
}



inline int getMult(int* i, int size)
{
    int mark;

    int mult = 1;

    if (size == 1) return mult;

    mark = i[0];
    for (int k=1; (k < size) && (i[k] == mark); k++)
    {
        mult++;
    }
    return mult;
}



inline int getMult(double* a, int size)
{
    double mark;

    int mult = 1;

    if (size == 1) return mult;

    mark = a[0];
    for (int k=1; (k < size) && (a[k] == mark); k++)
    {
        mult++;
    }
    return mult;
}


double OSRand();


double OSiRand(int iMin, int iMax);

#endif