sph.cuh

Go to the documentation of this file.

 
#ifndef MILUPHPC_SPH_CUH
#define MILUPHPC_SPH_CUH
 
#include "../particles.cuh"
#include "../subdomain_key_tree/subdomain.cuh"
#include "../parameter.h"
#include "../helper.cuh"
#include "../materials/material.cuh"
#include <float.h>
 
#include <boost/mpi.hpp>
#include <assert.h>
 
namespace SPH {
 
    void exchangeParticleEntry(SubDomainKeyTree *subDomainKeyTree, real *entry, real *toSend, integer *sendLengths,
                               integer *receiveLengths, integer numParticlesLocal);
 
    namespace Kernel {
 
        __global__ void fixedRadiusNN_bruteForce(Tree *tree, Particles *particles, integer *interactions,
                                                 integer numParticlesLocal, integer numParticles, integer numNodes);
 
        __global__ void fixedRadiusNN(Tree *tree, Particles *particles, integer *interactions, real radius,
                                      integer numParticlesLocal, integer numParticles, integer numNodes);
 
        __global__ void fixedRadiusNN_withinBox(Tree *tree, Particles *particles, integer *interactions,
                                                integer numParticlesLocal, integer numParticles, integer numNodes);
 
        __global__ void fixedRadiusNN_sharedMemory(Tree *tree, Particles *particles, integer *interactions,
                                                   integer numParticlesLocal, integer numParticles, integer numNodes);
 
        __global__ void fixedRadiusNN_variableSML(Material *materials, Tree *tree, Particles *particles, integer *interactions,
                                                  integer numParticlesLocal, integer numParticles,
                                                  integer numNodes);
 
        __device__ void redoNeighborSearch(Tree *tree, Particles *particles, int particleId,
                                             int *interactions, real radius, integer numParticles, integer numNodes);
 
        __global__ void compTheta(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                  DomainList *lowestDomainList, Curve::Type curveType);
 
        __global__ void symbolicForce(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                      DomainList *lowestDomainList, integer *sendIndices, real searchRadius,
                                      integer n, integer m, integer relevantIndex,
                                      Curve::Type curveType);
 
        __global__ void symbolicForce_test(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                           DomainList *lowestDomainList, integer *sendIndices, real searchRadius,
                                           integer n, integer m, integer relevantProc, integer relevantIndicesCounter,
                                           integer relevantIndexOld, Curve::Type curveType);
 
        __global__ void symbolicForce_test2(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                            DomainList *domainList, integer *sendIndices, real searchRadius,
                                            integer n, integer m, integer relevantProc, integer relevantIndicesCounter,
                                            Curve::Type curveType);
 
        __global__ void collectSendIndices(Tree *tree, Particles *particles, integer *sendIndices,
                                           integer *particles2Send, integer *particlesCount,
                                           integer n, integer length, Curve::Type curveType);
 
        __global__ void collectSendIndices_test2(Tree *tree, Particles *particles, integer *sendIndices,
                                                 integer *particles2Send, integer *particlesCount,
                                                 integer numParticlesLocal, integer numParticles,
                                                 integer treeIndex, int currentProc, Curve::Type curveType);
 
        __global__ void particles2Send(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                       DomainList *domainList, DomainList *lowestDomainList, integer maxLevel,
                                       integer *toSend, integer *sendCount, integer *alreadyInserted,
                                       integer insertOffset,
                                       integer numParticlesLocal, integer numParticles, integer numNodes, real radius,
                                       Curve::Type curveType = Curve::lebesgue);
 
        __global__ void collectSendIndicesBackup(integer *toSend, integer *toSendCollected, integer count);
 
        __global__ void collectSendEntriesBackup(SubDomainKeyTree *subDomainKeyTree, real *entry, real *toSend,
                                           integer *sendIndices, integer *sendCount, integer totalSendCount,
                                           integer insertOffset);
 
        __global__ void insertReceivedParticles(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                                DomainList *domainList, DomainList *lowestDomainList, int n, int m);
 
 
        __global__ void calculateCentersOfMass(Tree *tree, Particles *particles, integer level);
 
        __global__ void determineSearchRadii(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                             DomainList *domainList, DomainList *lowestDomainList, real *searchRadii,
                                             int n, int m, Curve::Type curveType);
 
        __global__ void info(Tree *tree, Particles *particles, Helper *helper,
                             integer numParticlesLocal, integer numParticles, integer numNodes);
 
        namespace Launch {
 
            real fixedRadiusNN_bruteForce(Tree *tree, Particles *particles, integer *interactions,
                                          integer numParticlesLocal, integer numParticles, integer numNodes);
 
            real fixedRadiusNN(Tree *tree, Particles *particles, integer *interactions, real radius,
                               integer numParticlesLocal, integer numParticles, integer numNodes);
 
            real fixedRadiusNN_sharedMemory(Tree *tree, Particles *particles, integer *interactions,
                                            integer numParticlesLocal, integer numParticles, integer numNodes);
 
            real fixedRadiusNN_withinBox(Tree *tree, Particles *particles, integer *interactions,
                                         integer numParticlesLocal, integer numParticles, integer numNodes);
 
            real fixedRadiusNN_variableSML(Material *materials, Tree *tree, Particles *particles, integer *interactions,
                                           integer numParticlesLocal, integer numParticles,
                                           integer numNodes);
 
            real compTheta(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                           DomainList *lowestDomainList, Curve::Type curveType);
 
            real symbolicForce(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                              DomainList *lowestDomainList, integer *sendIndices, real searchRadius,
                              integer n, integer m, integer relevantIndex,
                              Curve::Type curveType);
 
            real symbolicForce_test(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                               DomainList *lowestDomainList, integer *sendIndices, real searchRadius,
                                               integer n, integer m, integer relevantProc, integer relevantIndicesCounter,
                                               integer relevantIndexOld, Curve::Type curveType);
 
            real symbolicForce_test2(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                     DomainList *domainList, integer *sendIndices, real searchRadius,
                                     integer n, integer m, integer relevantProc, integer relevantIndicesCounter,
                                     Curve::Type curveType);
 
            real collectSendIndices(Tree *tree, Particles *particles, integer *sendIndices,
                                               integer *particles2Send, integer *particlesCount,
                                               integer n, integer length, Curve::Type curveType);
 
            real collectSendIndices_test2(Tree *tree, Particles *particles, integer *sendIndices,
                                          integer *particles2Send, integer *particlesCount,
                                          integer numParticlesLocal, integer numParticles,
                                          integer treeIndex, int currentProc, Curve::Type curveType);
 
            real particles2Send(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                           DomainList *domainList, DomainList *lowestDomainList, integer maxLevel,
                                           integer *toSend, integer *sendCount, integer *alreadyInserted,
                                           integer insertOffset,
                                           integer numParticlesLocal, integer numParticles, integer numNodes, real radius,
                                           Curve::Type curveType = Curve::lebesgue);
 
            real collectSendIndicesBackup(integer *toSend, integer *toSendCollected, integer count);
 
            real collectSendEntriesBackup(SubDomainKeyTree *subDomainKeyTree, real *entry, real *toSend, integer *sendIndices,
                                    integer *sendCount, integer totalSendCount, integer insertOffset);
 
            real insertReceivedParticles(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                         DomainList *domainList, DomainList *lowestDomainList, int n, int m);
 
            real info(Tree *tree, Particles *particles, Helper *helper,
                                 integer numParticlesLocal, integer numParticles, integer numNodes);
 
            real calculateCentersOfMass(Tree *tree, Particles *particles, integer level);
 
            real determineSearchRadii(SubDomainKeyTree *subDomainKeyTree, Tree *tree, Particles *particles,
                                      DomainList *domainList, DomainList *lowestDomainList, real *searchRadii,
                                      int n, int m, Curve::Type curveType);
        }
    }
 
}
 
#endif //MILUPHPC_SPH_CUH