LinkList.hpp

Go to the documentation of this file.

/*
 *  This file is part of AQUAgpusph, a free CFD program based on SPH.
 *  Copyright (C) 2012  Jose Luis Cercos Pita <jl.cercos@upm.es>
 *
 *  AQUAgpusph 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.
 *
 *  AQUAgpusph 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 AQUAgpusph.  If not, see <http://www.gnu.org/licenses/>.
 */

 
#ifndef LINKLIST_H_INCLUDED
#define LINKLIST_H_INCLUDED
 
#include "aquagpusph/sphPrerequisites.hpp"
#include <vector>
#include "Tool.hpp"
#include "Reduction.hpp"
 
namespace Aqua {
namespace CalcServer {

class LinkList : public Aqua::CalcServer::Tool
{
  public:
        LinkList(const std::string tool_name,
                 const std::string input = "r",
                 const std::string input_min = "r_min",
                 const std::string input_max = "r_max",
                 const std::string ihoc = "ihoc",
                 const std::string icell = "icell",
                 const std::string n_cells = "n_cells",
                 const std::string permutations = "id_unsorted",
                 const std::string inv_permutations = "id_sorted",
                 bool recompute_grid = true,
                 const std::string sorter = "radix-sort",
                 bool once = false);

        ~LinkList();

        void setup();

        inline cl_event getUserEvent() const { return _user_event; }
 
  protected:
        template<typename T>
        void nCells();

        void allocate();

        void setVariables();

        cl_event _execute(const std::vector<cl_event> events);
 
  private:
        inline uivec4 nCells32() const {
                uivec4 n_cells;
                for (unsigned int i = 0; i < 4; i++)
                        n_cells.s[i] = narrow_cast<uicl>(_n_cells.s[i]);
                return n_cells;
        }

        void setupOpenCL();

        std::string _input_name;

        std::string _input_min_name;

        std::string _input_max_name;

        std::string _ihoc_name;

        std::string _icell_name;

        std::string _n_cells_name;

        float _cell_length;

        ulvec4 _n_cells;

        Reduction* _min_pos;

        Reduction* _max_pos;

        bool _recompute_minmax;

        Tool* _sort;

        cl_kernel _ihoc;
        size_t _ihoc_lws;
        size_t _ihoc_gws;
        std::vector<void*> _ihoc_args;

        cl_kernel _icell;
        size_t _icell_lws;
        size_t _icell_gws;
        std::vector<void*> _icell_args;

        cl_kernel _ll;
        size_t _ll_lws;
        size_t _ll_gws;
        std::vector<void*> _ll_args;

        cl_event _user_event;
};
 
}
} // namespace
 
#endif // LINKLIST_H_INCLUDED