raw.h

Go to the documentation of this file.

/* Copyright (C) 2020
 * Richard Larsson <larsson@mps.mpg.de>
 * 
 * This program 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, or (at your option) any
 * later version.
 * 
 * This program 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 this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
 * USA. */

#ifndef RAW_H
#define RAW_H

#include "covariance_matrix.h"
#include "matpackI.h"

constexpr Numeric calibration(Numeric pc, Numeric pa, Numeric ph, Numeric tc, Numeric th) noexcept
{
  return tc + (th - tc) * (pa - pc) / (ph - pc);
}

constexpr Numeric systemtemp(Numeric pc, Numeric ph, Numeric tc, Numeric th) noexcept
{
  return (th*pc - tc*ph) / (ph - pc);
}

namespace linalg {
void avg(VectorView y, const ArrayOfVector& ys, const Index start=0, const Index end=-1);

void std(VectorView std, const Vector& y, const ArrayOfVector& ys, const Index start=0, const Index end=-1);

void var(VectorView var, const Vector& y, const ArrayOfVector& ys, const Index start=0, const Index end=-1);

void cov(MatrixView cov, const Vector& y, const ArrayOfVector& ys, const Index start=0, const Index end=-1);

Numeric median(const ConstVectorView v, const ArrayOfIndex& pos=ArrayOfIndex{});
};  // linalg

#endif  // RAW_H