gravitation_common.hpp

Go to the documentation of this file.

/*************************************************************************
 * Copyright (C) 2010 by Saleh Dindar and the Swarm-NG Development Team  *
 *                                                                       *
 * 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 3 of the License.        *
 *                                                                       *
 * 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.             *
 ************************************************************************/

#pragma once

// For using CoalescedStructArray in the gravitation implementations
#include "../types/coalescedstructarray.hpp"

#include "bppt.hpp"

namespace swarm { namespace gpu { namespace bppt {

template<int W>
struct GravitationAccScalars {
        static const int CHUNK_SIZE = W;
        typedef double scalar_t; 

        double _acc[CHUNK_SIZE];

        // Accessors
        GENERIC double& acc() { return _acc[0];  }
};


template<int W>
struct GravitationAccJerkScalars {
        static const int CHUNK_SIZE = W;
        typedef double scalar_t; 

        double _acc[CHUNK_SIZE];
        double _jerk[CHUNK_SIZE];

        // Accessors
        GENERIC double& acc() { return _acc[0];  }
        GENERIC double& jerk() { return _jerk[0];  }
};


GENERIC double inner_product(const double a[3],const double b[3]){
        return a[0]*b[0]+a[1]*b[1]+a[2]*b[2];
}

#include "pair_calculation.hpp"


} } } // end of namespace bppt :: gpu :: swarm