basicmath.h

//SPDX-License-Identifier:         LGPL-3.0-or-later
/* * eisgenerator - a shared library and application to generate EIS spectra
 * Copyright (C) 2022-2024 Carl Philipp Klemm <carl@uvos.xyz>
 *
 * This file is part of eisgenerator.
 *
 * eisgenerator is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * eisgenerator 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 Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with eisgenerator.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#pragma once
#include <vector>
#include "eistype.h"
 
namespace eis
{
        std::complex<fvalue> absGrad(const std::vector<eis::DataPoint>& data, size_t index);
 
        fvalue grad(const std::vector<fvalue>& data, const std::vector<fvalue>& omega, size_t index);
 
        std::complex<fvalue> grad(const std::vector<eis::DataPoint>& data, size_t index);
 
        fvalue mean(const std::vector<fvalue>& data);
 
        std::complex<fvalue> mean(const std::vector<eis::DataPoint>& data);
 
        fvalue median(std::vector<fvalue> data);
 
        std::complex<fvalue> median(const std::vector<eis::DataPoint>& data);
 
 
        void difference(std::vector<eis::DataPoint>& a, const std::vector<eis::DataPoint>& b);
 
        std::vector<eis::DataPoint> rescale(const std::vector<eis::DataPoint>& data, size_t outputSize);
 
        fvalue pearsonCorrelation(const std::vector<eis::DataPoint>& data);
 
        fvalue nyquistAreaVariance(const std::vector<eis::DataPoint>& data, eis::DataPoint* centroid = nullptr);
 
        fvalue maximumNyquistJump(const std::vector<eis::DataPoint>& data);
 
        void noise(std::vector<eis::DataPoint>& data, double amplitude, bool relative);
 
        void removeDuplicates(std::vector<eis::DataPoint>& data);
 
        bool fvalueEq(fvalue a, fvalue b, unsigned int ulp = 4);
 
        std::vector<eis::DataPoint> fitToFrequencies(std::vector<fvalue> omegas,
                                                     const std::vector<eis::DataPoint>& data,
                                                     bool linearExtrapolation = false);
 
}