GContainers.h

Go to the documentation of this file.

/*
 * GContainers.h
 *
 * Copyright 2016 Antonio Augusto Alves Junior
 *
 * Created on : Feb 25, 2016
 *      Author: Antonio Augusto Alves Junior
 */

/*
 *   This file is part of MCBooster.
 *
 *   MCBooster 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 3 of the License, or
 *   (at your option) any later version.
 *
 *   MCBooster 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 MCBooster.  If not, see <http://www.gnu.org/licenses/>.
 */

#ifndef GCONTAINERS_H_
#define GCONTAINERS_H_


#include <mcbooster/Config.h>
#include <mcbooster/Vector3R.h>
#include <mcbooster/Vector4R.h>
#include <vector>
#include <mcbooster/GTypes.h>
#include <thrust/device_vector.h>
#include <thrust/host_vector.h>
#include <thrust/complex.h>

#if !(THRUST_DEVICE_SYSTEM==THRUST_DEVICE_BACKEND_OMP || THRUST_DEVICE_SYSTEM==THRUST_DEVICE_BACKEND_TBB)
#include <thrust/system/cuda/experimental/pinned_allocator.h>
#endif

using namespace std;

namespace MCBooster
{

#if (THRUST_DEVICE_SYSTEM==THRUST_DEVICE_BACKEND_OMP || THRUST_DEVICE_SYSTEM==THRUST_DEVICE_BACKEND_TBB)

    template <typename T>
        using  mc_device_vector = thrust::host_vector<T>;
    template <typename T>
        using  mc_host_vector   = thrust::host_vector<T>;

#else

    template <typename T>
        using  mc_device_vector = thrust::device_vector<T>;
    template <typename T>
        using  mc_host_vector = thrust::host_vector<T,
                thrust::cuda::experimental::pinned_allocator<T>>;

#endif


//-----------------------------------------------------------------------
//complex number container
typedef thrust::complex<GReal_t> GComplex_t; 
//-----------------------------------------------------------------------

typedef mc_host_vector<Vector4R> FourVectors_h; 
typedef mc_host_vector<Vector3R> ThreeVectors_h; 
//-----------------------------------------------------------------------
//basic containers on host

typedef mc_host_vector<GBool_t>    BoolVector_h; 
typedef mc_host_vector<GReal_t>    RealVector_h; 
typedef mc_host_vector<GComplex_t> ComplexVector_h; 
typedef mc_host_vector<Vector4R>   Particles_h;
typedef vector<Particles_h*>       ParticlesSet_h;
typedef vector<RealVector_h*>      VariableSet_h;
//-----------------------------------------------------------------------
//basic containers on device
typedef mc_device_vector<GBool_t>    BoolVector_d; 
typedef mc_device_vector<GReal_t>    RealVector_d; 
typedef mc_device_vector<GComplex_t> ComplexVector_d; 
typedef mc_device_vector<Vector4R>   Particles_d; 
typedef vector<Particles_d*>  ParticlesSet_d; 
typedef vector<RealVector_d*> VariableSet_d; 
typedef thrust::detail::tuple_of_iterator_references<Vector4R &,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type> GT1;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type> GT2;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GT3;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GT4;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GT5;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type> GT6;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        thrust::null_type, thrust::null_type, thrust::null_type> GT7;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        thrust::null_type, thrust::null_type> GT8;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        Vector4R &, thrust::null_type> GT9;

typedef thrust::detail::tuple_of_iterator_references<Vector4R &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        Vector4R &, Vector4R &> GT10;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type> GTR2;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GTR3;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, thrust::null_type, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GTR4;
typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type,
        thrust::null_type> GTR5;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, thrust::null_type,
        thrust::null_type, thrust::null_type, thrust::null_type> GTR6;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        thrust::null_type, thrust::null_type, thrust::null_type> GTR7;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        thrust::null_type, thrust::null_type> GTR8;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        Vector4R &, thrust::null_type> GTR9;

typedef thrust::detail::tuple_of_iterator_references<GReal_t &, Vector4R &,
        Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &, Vector4R &,
        Vector4R &, Vector4R &> GTR10;

}
#endif /* GCONTAINERS_H_ */