Aurora/src/Function1D.h

#ifndef AURORA_FUNCTION1D_H
#define AURORA_FUNCTION1D_H

#include "Matrix.h"

namespace Aurora {

    enum InterpnMethod
    {
        Spline=0,Linear
    };

    Matrix complex(const Matrix& matrix);

    Matrix real(const Matrix& matrix);

    Matrix imag(const Matrix& matrix);

    Matrix ceil(const Matrix& matrix);

    Matrix ceil(const Matrix&& matrix);

    Matrix round(const Matrix& matrix);

    Matrix round(const Matrix&& matrix);

    /**
     * 开根号，暂时只支持正整数
     * @param matrix
     * @return
     */
    Matrix sqrt(const Matrix& matrix);

    Matrix sqrt(Matrix&& matrix);

    Matrix abs(const Matrix& matrix);

    Matrix abs(Matrix&& matrix);

    Matrix sign(const Matrix& matrix);

    Matrix sign(Matrix&& matrix);

    Matrix interp1(const Matrix& aX, const Matrix& aV, const Matrix& aX1, InterpnMethod aMethod);

    Matrix repmat(const Matrix& aMatrix,int aRowTimes, int aColumnTimes);

    Matrix repmat(const Matrix& aMatrix,int aRowTimes, int aColumnTimes, int aSliceTimes);

    Matrix log(const Matrix& aMatrix, int aBaseNum = -1);

    Matrix exp(const Matrix& aMatrix);

    Matrix mod(const Matrix& aMatrix, double aValue);

    Matrix acos(const Matrix& aMatrix);

    Matrix acosd(const Matrix& aMatrix);

    /**
     * 多项式计算
     * @brief 例如p[1 0 1],x[3 2 5],代表对多项式 y = x^2 + 1 求(x=3, x=2, x=5)时所有的y
     * @attention 本函数的数据一律视为向量，在计算时会完全忽视维度，值考虑作为数组输入
     * @param aP 多项式系数，指定为向量
     * @param aX 查询点，指定为向量
     * @return 查询结果
     */
    Matrix polyval(const Matrix& aP, const Matrix& aX);
};


#endif //AURORA_FUNCTION1D_H