SymMatrixSolver.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef ALICEO2_GPUCOMMON_TPCFASTTRANSFORMATION_SYMMATRIXSOLVER_H
#define ALICEO2_GPUCOMMON_TPCFASTTRANSFORMATION_SYMMATRIXSOLVER_H
 
#include "GPUCommonDef.h"
#include "GPUCommonRtypes.h"
#include <vector>
#include <cassert>
#include <algorithm>
 
namespace o2
{
namespace gpu
{
 
class SymMatrixSolver
{
 public:
  SymMatrixSolver(int32_t N, int32_t M) : mN(N), mM(M), mShift(mN + mM)
  {
    assert(N > 0 && M > 0);
    mA.resize(mN * mShift, 0.);
  }
 
  double& A(int32_t i, int32_t j)
  {
    auto ij = std::minmax(i, j);
    assert(ij.first >= 0 && ij.second < mN);
    return mA[ij.first * mShift + ij.second];
  }
 
  double& B(int32_t i, int32_t j)
  {
    assert(i >= 0 && i < mN && j >= 0 && j < mM);
    return mA[i * mShift + mN + j];
  }
 
  void solve();
 
  void print();
 
  static int32_t test(bool prn = 0);
 
 private:
 private:
  int32_t mN = 0;
  int32_t mM = 0;
  int32_t mShift = 0;
  std::vector<double> mA;
 
  ClassDefNV(SymMatrixSolver, 0);
};
 
} // namespace gpu
} // namespace o2
 
#endif