diff --git a/Addons/FmmApi/CMakeLists.txt b/Addons/FmmApi/CMakeLists.txt deleted file mode 100644 index 249a19cb455c3e3765b9e2f52e4ca666d0350233..0000000000000000000000000000000000000000 --- a/Addons/FmmApi/CMakeLists.txt +++ /dev/null @@ -1,84 +0,0 @@ -# check if compiling into source directories -STRING(COMPARE EQUAL "${CMAKE_SOURCE_DIR}" "${CMAKE_BINARY_DIR}" insource) -if(insource) - MESSAGE(FATAL_ERROR "${PROJECT_NAME} requires an out of source build. Goto ./Build and tapes cmake ../") -endif(insource) - -project(Addons_fmmapi_scalfmm CXX) -ADD_DEFINITIONS( ${ScaLFMM_CXX_FLAGS}) - -# Active language -# ----------------------- -ENABLE_LANGUAGE(CXX ) -MESSAGE(STATUS " CXX ${CMAKE_CXX_COMPILER_ID}" ) - -# Options -OPTION( SCALFMM_ADDON_FMMAPI "Set to ON to build ScaFMM FMM API interface" OFF ) - -# if ask to build addon -if(SCALFMM_ADDON_FMMAPI) - # first build lib scalfmmapi - set(LIBRARY_OUTPUT_PATH ../lib/${CMAKE_BUILD_TYPE}) - - # Searching all cpp file - file( GLOB_RECURSE source_lib_files Src/*.cpp ) - - # Adding cpp files to project - add_library( scalfmmapi ${source_lib_files} ) - - # Add blas library (even if it is set to off) - target_link_libraries( scalfmmapi scalfmm) - - # Adding the entire project dir as an include dir - INCLUDE_DIRECTORIES( - ${SCALFMM_BINARY_DIR}/Src - ${SCALFMM_SOURCE_DIR}/Src - ${SCALFMM_INCLUDES} - ) - - # Install lib - install( TARGETS scalfmmapi ARCHIVE DESTINATION lib LIBRARY DESTINATION lib ) - - # Install header - SET(my_include_dirs "Src") - - file( GLOB hpp_in_dir Src/*.hpp ) - INSTALL( FILES ${hpp_in_dir} DESTINATION include/ScalFmm/FmmApi ) - - file( GLOB_RECURSE source_tests_files Tests/*.cpp ) - INCLUDE_DIRECTORIES( ${SCALFMM_BINARY_DIR}/Src ) - - # Then build test files - foreach(exec ${source_tests_files}) - get_filename_component( - execname ${exec} - NAME_WE - ) - - set(compile_exec "TRUE") - - foreach(fuse_key ${FUSE_LIST}) - file(STRINGS "${exec}" lines_fuse REGEX "@FUSE_${fuse_key}") - if(lines_fuse) - if( NOT SCALFMM_USE_${fuse_key} ) - MESSAGE( STATUS "This needs ${fuse_key} = ${exec}" ) - set(compile_exec "FALSE") - endif() - endif() - endforeach() - - # Dependency are OK - if( compile_exec ) - add_executable( ${execname} ${exec} ) - # link to scalfmm and scalfmmapi - target_link_libraries( - ${execname} - ${scalfmm_lib} - scalfmmapi - ${BLAS_LIBRARIES} - ${LAPACK_LIBRARIES} - ${SCALFMM_LIBRARIES} - ) - endif() - endforeach(exec) -endif() diff --git a/Addons/FmmApi/Src/FmmApi.h b/Addons/FmmApi/Src/FmmApi.h deleted file mode 100644 index 3550ebb6ea823257aadd50ea405bc4586fc3185b..0000000000000000000000000000000000000000 --- a/Addons/FmmApi/Src/FmmApi.h +++ /dev/null @@ -1,114 +0,0 @@ -// =================================================================================== -// Logiciel initial: ScalFmm Version 0.5 -// Co-auteurs : Olivier Coulaud, Bérenger Bramas. -// Propriétaires : INRIA. -// Copyright © 2011-2012, diffusé sous les termes et conditions d’une licence propriétaire. -// Initial software: ScalFmm Version 0.5 -// Co-authors: Olivier Coulaud, Bérenger Bramas. -// Owners: INRIA. -// Copyright © 2011-2012, spread under the terms and conditions of a proprietary license. -// =================================================================================== -#ifndef FMMAPI_H -#define FMMAPI_H - -enum FmmApiErrors { - FMMAPI_NO_ERROR, - FMMAPI_SUPPORTED_PARAMETER, - FMMAPI_UNSUPPORTED_PARAMETER, - FMMAPI_UNKNOWN_PARAMETER -}; - -////////////////////// Opérateurs FMM Core : ////////////////////////// - -enum FmmApiCoreParameters { - FMMCORE_TREE_HEIGHT, // hombre de niveaux de l'arbre (int) - FMMCORE_ROOT_BOX_WIDTH, // taille de la boîte racine (FReal) - FMMCORE_ROOT_BOX_CENTER, // position du centre de la boîte racine (FReal[3]) - FMMCORE_LEAF_BOX_WIDTH, // taille des boîtes feuilles (FReal) - FMMCORE_POINTS_PER_LEAF, // nombre moyen de points par feuille (FReal) - FMMCORE_MPI_COMMUNICATOR, // communicateur MPI (MPI_Comm) - FMMCORE_THREADS_NUMBER, // nombre de threads (int) - FMMCORE_THREAD_ID, // id du thread (int) - FMMCORE_RHS_NUMBER, // nombre de seconds membres (int) - //paramètres en lecture seule : - FMMCORE_HANDLES_P2P // renvoie 0 ou 1 pour dire si le FmmCore gère ou pas le P2P. -}; - -int FmmCore_init(void **fmmCore) ; /*alloue et initialise le FmmCore*/ -int FmmCore_free(void *fmmCore) ; /*libère le FmmCore*/ -int FmmCore_isParameterUsed(void */*fmmCore*/, int *name, int *flag); -int FmmCore_setParameter(void *fmmCore, int *name, void*value); -int FmmCore_setParameter(void *fmmCore, int name, void*value); -int FmmCore_getParameter(void *fmmCore, int *name, void*value); -int FmmCore_getParameter(void *fmmCore, int name, void*value); - - -int FmmCore_getRadius(void*fmmCore, void *boxId, FReal *radius);/*Renvoie le rayon de la boîte*/ -int FmmCore_getCentre(void*/*fmmCore*/, void *boxId, FReal **centre); /*Renvoie dans le FReal[3] centre les coordonnées du centre*/ -int FmmCore_getLevel(void*/*fmmCore*/, void *boxId, int *level); /*Renvoie dans level le niveau de la boîte boxId dans son arbre*/ -int FmmCore_getMultipoleArray(void* /*fmmCore*/, void *boxId, void **F); /*Renvoie dans F l'adresse où stocker l'expansion multipôle associée à la boîte boxId*/ -int FmmCore_getLocalArray(void* /*fmmCore*/, void *boxId, void **F); /*Renvoie dans F l'adresse où stocker l'expansion locale associée à la boîte boxId*/ -int FmmCore_getCoord(void*/*fmmCore*/, void *boxId, int *coord); /*Renvoie dans coord la position dans l'arbre*/ - - -/* Données potentiel/champ */ -int FmmCore_getSource(void* /*fmmCore*/, void *boxId, FReal** position, void** potential, int *number); /* Appelé par P2P et P2M pour obtenir le nombre, la position et le potentiel des sources. -Les différents tableaux sont (éventuellement) alloués par le FmmCore. */ -int FmmCore_releaseSource(void*fmmCore, void *boxId, void* potential, FReal* position); /* si le core veut libérer ces tableaux potentiel et position.*/ -int FmmCore_getTargetPoints(void* /*fmmCore*/, void *boxId, FReal** position, int *number) ; /* : Appelé par P2P et L2P pour obtenir le nombre et la position des points cibles.*/ -int FmmCore_releaseTargetPoints(void*fmmCore, void *boxId, FReal* position); /* si le core veut libérer ce tableau "position".*/ -int FmmCore_getTargetField(void* /*fmmCore*/, void *boxId, FReal* F); /* obtient dans un tableau F alloué/libéré par L2P/P2P les valeurs des champs aux points cible -(pour le cas où P2P et L2P doivent sommer leurs résultats).*/ -int FmmCore_setTargetField(void* /*fmmCore*/, void *boxId, FReal* F); /* transmets au FmmCore dans F les valeurs des champs aux points cibles mis à jour.*/ - -/* Entrée/sortie principale */ -int FmmCore_setKernelData(void *fmmCore, void *fmmKernel); /* stocke l'identifiant du FmmKernel dans le FmmCore. Cela permet par la suite aux -opérateurs FMM d'accéder au FmmKernel, et donc d'accéder aux données spécifiques au kernel -(p.ex. fréquence dans le cas Helmholtz, …)*/ -int FmmCore_getKernelData(void*fmmCore, void **fmmKernel); /* récupère l'identifiant du FmmKernel. */ -int FmmCore_setPositions(void *fmmCore, int *nb, FReal *position) ; /* transmet au FmmCore les potentiels associés aux points sources. -Le tableau potential est alloué et libéré par la routine appelant, le FmmCore doit donc en faire une copie.*/ -int FmmCore_setPotentials(void *fmmCore, void *potentials); /* transmet au FmmCore les potentiels associés aux points sources. -Le tableau potential est alloué et libéré par la routine appelant, le FmmCore doit donc en faire une copie.*/ -int FmmCore_doComputation(void *fmmCore) ; /* réalise le produit multipôle. */ -/* !!! Warning use *filed and not **field */ -int FmmCore_getField(void *fmmCore, void *fields) ;/* récupère après le produit multipôle la valeur des champs en chaque point. -Le tableau field doit être alloué et libéré par la routine appelante. */ - -////////////////////// Opérateurs FMM Kernel : ////////////////////////// - -enum FmmApiKernelParameters { - FMMKERNEL_ACCURACY, // précision demandée à la FMM : 1e-3, 1e-6, ... (FReal) - FMMKERNEL_POTENTIAL_DATA_SIZE, // taille en octet de la donnée "potentiel" pour 1 point source (int) - FMMKERNEL_FIELD_DATA_SIZE, // taille en octet de la donnée "field" pour 1 point cible (int) - //paramètres en lecture seule : - FMMKERNEL_HANDLES_P2P // renvoie 0 ou 1 pour dire si le FmmKernel gère ou pas le P2P. -}; - -/******* Allocation : ******/ -int FmmKernel_init(void *fmmCore, void **fmmKernel);/* : alloue et initialise le FmmKernel */ -int FmmKernel_free(void *fmmKernel); /* libére le FmmKernel */ - -/******* Configuration : ***/ - -int FmmKernel_isParameterUsed(void * /*fmm*/, int *name, int *flag); -int FmmKernel_setParameter(void *fmmKernel, int *name, void*value); -int FmmKernel_setParameter(void *fmmKernel, int name, void*value); -int FmmKernel_getParameter(void *fmmKernel, int *name, void*value); -int FmmKernel_getParameter(void *fmmKernel, int name, void*value); - -/****** Données FMM : *****/ -int FmmKernel_getMultipoleArraySize(void *fmmCore, int *size); /* Renvoie dans size la taille (en octets) de l'expansion multipôle associée à la boîte boxId */ -int FmmKernel_getLocalArraySize(void *fmmCore, int *size); /* Renvoie dans size la taille (en octets) de l'expansion locale associée à la boîte boxId*/ - -/******* Opérateurs FMM : **/ -int FmmKernel_P2M(void *fmmCore, void* boxId); -int FmmKernel_L2P(void *fmmCore, void* boxId); -int FmmKernel_M2M(void *fmmCore, void *boxIdFather, void *boxIdSon); -int FmmKernel_L2L(void *fmmCore, void *boxIdFather, void *boxIdSon); -int FmmKernel_M2L(void *fmmCore, void *boxIdSrc, void *boxIdDest); -int FmmKernel_P2P_inner(void *fmmCore, void *boxIdSrcDest); -int FmmKernel_P2P(void *fmmCore, void *boxIdSrc, void *boxIdDest); /* pas mutuel, i.e. on fait seulement dans 1 sens. */ - - -#endif // FMMAPI_H diff --git a/Addons/FmmApi/Src/ScalfmmApiCore.cpp b/Addons/FmmApi/Src/ScalfmmApiCore.cpp deleted file mode 100644 index bd0a87cc7b5cbc13ed31148fc5108cff04a88165..0000000000000000000000000000000000000000 --- a/Addons/FmmApi/Src/ScalfmmApiCore.cpp +++ /dev/null @@ -1,662 +0,0 @@ -// =================================================================================== -// Logiciel initial: ScalFmm Version 0.5 -// Co-auteurs : Olivier Coulaud, Bérenger Bramas. -// Propriétaires : INRIA. -// Copyright © 2011-2012, diffusé sous les termes et conditions d’une licence propriétaire. -// Initial software: ScalFmm Version 0.5 -// Co-authors: Olivier Coulaud, Bérenger Bramas. -// Owners: INRIA. -// Copyright © 2011-2012, spread under the terms and conditions of a proprietary license. -// =================================================================================== -#include "../../Src/Containers/FOctree.hpp" -#include "../../Src/Containers/FVector.hpp" - -#include "../../Src/Components/FSimpleLeaf.hpp" -#include "../../Src/Components/FBasicCell.hpp" - -#include "../../Src/Utils/FPoint.hpp" - -#include "../../Src/Core/FFmmAlgorithm.hpp" -#include "../../Src/Core/FFmmAlgorithmThread.hpp" -#include "../../Src/Core/FFmmAlgorithmTask.hpp" - -#include "../../Src/Components/FBasicKernels.hpp" - -#include "../../Src/Components/FBasicParticleContainer.hpp" - -#ifdef SCALFMM_USE_MPI -#include "../../Src/Utils/FMpi.hpp" -#endif - -#include "FmmApi.h" - -template -class CoreCell : public FBasicCell { - char* multipole; - char* local; - int level; - FReal positions[3]; - mutable ContainerClass* container; - -public: - CoreCell() : multipole(nullptr), local(nullptr), level(0), container(nullptr) { - } - void createArrays(const int multipoleSize, const int localSize){ - multipole = new char[multipoleSize]; - local = new char[localSize]; - - memset(multipole, 0, multipoleSize); - memset(local, 0, localSize); - } - ~CoreCell(){ - delete[] multipole; - delete[] local; - } - - const void* getMultipole() const{ - return multipole; - } - - const void* getLocal() const{ - return local; - } - - void* getMultipole(){ - return multipole; - } - - void* getLocal(){ - return local; - } - - void setLevel(const int inLevel){ - level = inLevel; - } - - int getLevel() const { - return level; - } - - void setPosition(const FReal inPositions[3]){ - positions[0] = inPositions[0]; - positions[1] = inPositions[1]; - positions[2] = inPositions[2]; - } - - const FReal* getPosition() const { - return positions; - } - - void setContainer(ContainerClass* inContainer) const { - container = inContainer; - } - - ContainerClass* getContainer() const { - return container; - } -}; - - - -template< class CellClass, class ContainerClass> -class CoreKernel : public FAbstractKernels { - void* fmmCore; - -public: - CoreKernel(void* inFmmCore): fmmCore(inFmmCore){ - } - - /** Default destructor */ - virtual ~CoreKernel(){ - } - - /** Do nothing */ - virtual void P2M(CellClass* const cell, const ContainerClass* const container) { - cell->setContainer(const_cast(container)); - FmmKernel_P2M(fmmCore,(void*) cell); - } - - /** Do nothing */ - virtual void M2M(CellClass* const FRestrict cell, const CellClass*const FRestrict *const FRestrict children, const int ) { - for(int idx = 0 ; idx < 8 ; ++idx){ - if( children[idx] ){ - FmmKernel_M2M(fmmCore, (void*)cell, (void*)children[idx]); - } - } - } - - /** Do nothing */ - virtual void M2L(CellClass* const FRestrict cell, const CellClass* interactions[], const int , const int ) { - for(int idx = 0 ; idx < 343 ; ++idx){ - if( interactions[idx] ){ - FmmKernel_M2L(fmmCore, (void*)cell, (void*)interactions[idx]); - } - } - } - - /** Do nothing */ - virtual void L2L(const CellClass* const FRestrict cell, CellClass* FRestrict *const FRestrict children, const int ) { - for(int idx = 0 ; idx < 8 ; ++idx){ - if( children[idx] ){ - FmmKernel_L2L(fmmCore, (void*)cell, (void*)children[idx]); - } - } - } - - /** Do nothing */ - virtual void L2P(const CellClass* const cell, ContainerClass* const container){ - cell->setContainer((ContainerClass*)container); - FmmKernel_L2P(fmmCore, (void*)cell); - } - - - /** Do nothing */ - virtual void P2P(const FTreeCoordinate& , - ContainerClass* const FRestrict targets, const ContainerClass* const FRestrict /*sources*/, - ContainerClass* const neighbors[27], const int ){ - CellClass* cell = (CellClass*)targets->getParentCell(); - FmmKernel_P2P_inner(fmmCore, cell); - - for(int idx = 0 ; idx < 27 ; ++idx){ - if( neighbors[idx] ){ - FmmKernel_P2P(fmmCore, (CellClass*)neighbors[idx]->getParentCell(), cell); - } - } - } - - /** Do nothing */ - virtual void P2PRemote(const FTreeCoordinate& , - ContainerClass* const FRestrict , const ContainerClass* const FRestrict , - ContainerClass* const [27], const int ){ - printf("Error remote not implemented!!!!\n"); - } - -}; - - -class CoreVector { - mutable void* parentCell; - void* fields; - int sizeOfField; - void* potentials; - int sizeOfPential; - FVector positions; - FVector indexes; -public: - CoreVector() : parentCell(nullptr), fields(nullptr), sizeOfField(0), - potentials(nullptr), sizeOfPential(0){ - } - - ~CoreVector(){ - delete[] (char*)fields; - delete[] (char*)potentials; - } - - void setParentCell(void* inCell) const{ - parentCell = inCell; - } - - void* getParentCell() const{ - return parentCell; - } - - void allocateFields(const int memSizePerParticles){ - if(fields) delete[](char*)fields; - fields = new char[getNbParticles() * memSizePerParticles]; - sizeOfField = memSizePerParticles; - memset(fields, 0, getNbParticles() * memSizePerParticles); - } - - void* getFields() const{ - return fields; - } - - int getSizeOfField() const{ - return sizeOfField; - } - - void allocatePotentials(const int memSizePerParticles){ - if(potentials) delete[](char*)potentials; - potentials = new char[getNbParticles() * memSizePerParticles]; - sizeOfPential = memSizePerParticles; - memset(potentials, 0, getNbParticles() * memSizePerParticles); - } - - void* getPotentials() const{ - return potentials; - } - - int getSizeOfPotential() const{ - return sizeOfPential; - } - - const int* getIndexes() const { - return indexes.data(); - } - - const FReal* getPositions() const { - return positions.data(); - } - - int getNbParticles() { - return indexes.getSize(); - } - - void push(const FPoint& partPosition, const int partIndex){ - positions.push(partPosition.getX()); - positions.push(partPosition.getY()); - positions.push(partPosition.getZ()); - indexes.push(partIndex); - } -}; - - -typedef CoreVector CoreContainerClass; - -typedef CoreCell CoreCellClass; -typedef FSimpleLeaf LeafClass; -typedef FOctree OctreeClass; -typedef CoreKernel CoreKernelClass; - -typedef FFmmAlgorithm FmmClass; -typedef FFmmAlgorithmThread FmmClassThread; - - -struct ScalFmmCoreHandle { - struct ScalFmmCoreConfig { - //paramètres en lecture/écriture : - int treeHeight; // hombre de niveaux de l'arbre (int) - FReal boxWidth; // taille de la boîte racine (FReal) - FReal boxCenter[3]; // position du centre de la boîte racine (FReal[3]) -#ifdef SCALFMM_USE_MPI - MPI_Comm mpiCom; // communicateur MPI (MPI_Comm) -#endif - int nbThreads; // nombre de threads (int) - int rhsNumber; // nombre de seconds membres (int) - }; - - ScalFmmCoreConfig config; - OctreeClass* octree; - void *kernelHandle; -}; - -int FmmCore_init(void **fmmCore) { - ScalFmmCoreHandle* corehandle = new ScalFmmCoreHandle; - memset(corehandle, 0, sizeof(corehandle)); - - corehandle->config.nbThreads = omp_get_max_threads(); - - *fmmCore = corehandle; - - return FMMAPI_NO_ERROR; -} /*alloue et initialise le FmmCore*/ - - -int FmmCore_free(void *fmmCore) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - if(corehandle->octree) delete corehandle->octree; - delete corehandle; - return FMMAPI_NO_ERROR; -} /*libère le FmmCore*/ - -int FmmCore_isParameterUsed(void */*fmmCore*/, int *name, int *flag){ - switch( *name ){ - case FMMCORE_ROOT_BOX_WIDTH : - case FMMCORE_ROOT_BOX_CENTER : - case FMMCORE_TREE_HEIGHT : -#ifdef SCALFMM_USE_MPI - case FMMCORE_MPI_COMMUNICATOR: -#endif - case FMMCORE_THREADS_NUMBER: - case FMMCORE_THREAD_ID: - case FMMCORE_RHS_NUMBER: - case FMMCORE_HANDLES_P2P: - *flag = FMMAPI_SUPPORTED_PARAMETER; - break; - case FMMCORE_LEAF_BOX_WIDTH : - case FMMCORE_POINTS_PER_LEAF : - *flag = FMMAPI_UNSUPPORTED_PARAMETER; - break; - default: - *flag = FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmCore_setParameter(void *fmmCore, int *name, void*value){ - int flag; - - FmmCore_isParameterUsed(fmmCore, name, &flag); - if( flag != FMMAPI_SUPPORTED_PARAMETER){ - return flag; - } - - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - switch( *name ){ - case FMMCORE_TREE_HEIGHT : - corehandle->config.treeHeight = *(int*)value ; - break; - case FMMCORE_ROOT_BOX_WIDTH : - corehandle->config.boxWidth = *(FReal*)value; - break; - case FMMCORE_ROOT_BOX_CENTER : - memcpy(corehandle->config.boxCenter, value, sizeof(FReal)*3); - break; -#ifdef SCALFMM_USE_MPI - case FMMCORE_MPI_COMMUNICATOR: - corehandle->config.mpiCom = *(MPI_Comm*)value; - break; -#endif - case FMMCORE_THREADS_NUMBER: - corehandle->config.nbThreads = *(int*)value; - break; - case FMMCORE_RHS_NUMBER: - corehandle->config.rhsNumber = *(int*)value; - break; - case FMMCORE_HANDLES_P2P: - case FMMCORE_THREAD_ID: - return FMMAPI_UNSUPPORTED_PARAMETER; - default: - return FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmCore_setParameter(void *fmmCore, int name, void*value){ - return FmmCore_setParameter(fmmCore, &name, value); -} - - -int FmmCore_getParameter(void *fmmCore, int *name, void*value){ - int flag; - - FmmCore_isParameterUsed(fmmCore, name, &flag); - if( flag != FMMAPI_SUPPORTED_PARAMETER){ - return flag; - } - - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - switch( *name ){ - case FMMCORE_TREE_HEIGHT : - *(int*)value = corehandle->config.treeHeight; - break; - case FMMCORE_ROOT_BOX_WIDTH : - *(FReal*)value = corehandle->config.boxWidth; - break; - case FMMCORE_ROOT_BOX_CENTER : - memcpy(value,corehandle->config.boxCenter, sizeof(FReal)*3); - break; -#ifdef SCALFMM_USE_MPI - case FMMCORE_MPI_COMMUNICATOR: - *(MPI_Comm*)value = corehandle->config.mpiCom; - break; -#endif - case FMMCORE_THREADS_NUMBER: - *(int*)value = corehandle->config.nbThreads; - break; - case FMMCORE_THREAD_ID: - *(int*)value = omp_get_thread_num(); - break; - case FMMCORE_RHS_NUMBER: - *(int*)value = corehandle->config.rhsNumber; - break; - case FMMCORE_HANDLES_P2P: - *(int*)value = true; - break; - default: - return FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmCore_getParameter(void *fmmCore, int name, void*value){ - return FmmCore_getParameter(fmmCore, &name, value); -} - - -int FmmCore_getRadius(void*fmmCore, void *boxId, FReal *radius) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - *radius = corehandle->octree->getBoxWidth() / FReal(1<getLevel()); - return FMMAPI_NO_ERROR; -} /*Renvoie le rayon de la boîte*/ - -int FmmCore_getCentre(void*/*fmmCore*/, void *boxId, FReal **centre) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - - (*centre)[0] = boxhandle->getPosition()[0]; - (*centre)[1] = boxhandle->getPosition()[1]; - (*centre)[2] = boxhandle->getPosition()[2]; - - return FMMAPI_NO_ERROR; -} /*Renvoie dans le FReal[3] centre les coordonnées du centre*/ - -int FmmCore_getLevel(void*/*fmmCore*/, void *boxId, int *level) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - *level = boxhandle->getLevel(); - return FMMAPI_NO_ERROR; -} /*Renvoie dans level le niveau de la boîte boxId dans son arbre*/ - -int FmmCore_getMultipoleArray(void* /*fmmCore*/, void *boxId, void **F) { - CoreCellClass* cell = (CoreCellClass*)boxId; - *F = cell->getMultipole(); - return FMMAPI_NO_ERROR; -} /*Renvoie dans F l'adresse où stocker l'expansion multipôle associée à la boîte boxId*/ - -int FmmCore_getLocalArray(void* /*fmmCore*/, void *boxId, void **F) { - CoreCellClass* cell = (CoreCellClass*)boxId; - *F = cell->getLocal(); - return FMMAPI_NO_ERROR; -} /*Renvoie dans F l'adresse où stocker l'expansion locale associée à la boîte boxId*/ - -int FmmCore_getCoord(void*/*fmmCore*/, void *boxId, int *coord) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - coord[0] = boxhandle->getCoordinate().getX(); - coord[1] = boxhandle->getCoordinate().getY(); - coord[2] = boxhandle->getCoordinate().getZ(); - return FMMAPI_NO_ERROR; -} /*Renvoie dans coord la position dans l'arbre*/ - - -/* Données potentiel/champ */ -int FmmCore_getSource(void* /*fmmCore*/, void *boxId, FReal** position, void** potential, int *number) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - CoreContainerClass* sources = boxhandle->getContainer(); - - *number = sources->getNbParticles(); - *position = const_cast(sources->getPositions()); - *potential = sources->getPotentials(); - - return FMMAPI_NO_ERROR; -} /* Appelé par P2P et P2M pour obtenir le nombre, la position et le potentiel des sources. -Les différents tableaux sont (éventuellement) alloués par le FmmCore. */ - -int FmmCore_releaseSource(void*fmmCore, void *boxId, void* potential, FReal* position) { - return FMMAPI_NO_ERROR; -} /* si le core veut libérer ces tableaux potentiel et position.*/ - -int FmmCore_getTargetPoints(void* /*fmmCore*/, void *boxId, FReal** position, int *number) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - CoreContainerClass* targets = boxhandle->getContainer(); - - *number = targets->getNbParticles(); - *position = const_cast(targets->getPositions()); - - return FMMAPI_NO_ERROR; -} /* : Appelé par P2P et L2P pour obtenir le nombre et la position des points cibles.*/ - -int FmmCore_releaseTargetPoints(void*fmmCore, void *boxId, FReal* position) { - return FMMAPI_NO_ERROR; -} /* si le core veut libérer ce tableau "position".*/ - -int FmmCore_getTargetField(void* /*fmmCore*/, void *boxId, FReal* F) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - CoreContainerClass* targets = boxhandle->getContainer(); - - memcpy(F, targets->getFields(), targets->getNbParticles() * targets->getSizeOfField()); - - return FMMAPI_NO_ERROR; -} /* obtient dans un tableau F alloué/libéré par L2P/P2P les valeurs des champs aux points cible -(pour le cas où P2P et L2P doivent sommer leurs résultats).*/ - -int FmmCore_setTargetField(void* /*fmmCore*/, void *boxId, FReal* F) { - CoreCellClass* boxhandle = (CoreCellClass*)boxId; - CoreContainerClass* targets = boxhandle->getContainer(); - - memcpy(targets->getFields(), F, targets->getNbParticles() * targets->getSizeOfField()); - - return FMMAPI_NO_ERROR; -} /* transmets au FmmCore dans F les valeurs des champs aux points cibles mis à jour.*/ - - - -/* Entrée/sortie principale */ -int FmmCore_setKernelData(void *fmmCore, void *fmmKernel) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - corehandle->kernelHandle = fmmKernel; - return FMMAPI_NO_ERROR; -} /* stocke l'identifiant du FmmKernel dans le FmmCore. Cela permet par la suite aux -opérateurs FMM d'accéder au FmmKernel, et donc d'accéder aux données spécifiques au kernel -(p.ex. fréquence dans le cas Helmholtz, …)*/ - -int FmmCore_getKernelData(void*fmmCore, void **fmmKernel) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - *fmmKernel = corehandle->kernelHandle; - return FMMAPI_NO_ERROR; -} /* récupère l'identifiant du FmmKernel. */ - - -int FmmCore_setPositions(void *fmmCore, int *nb, FReal *position) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - - corehandle->octree = new OctreeClass(corehandle->config.treeHeight, FMath::Min(3,corehandle->config.treeHeight-1), - corehandle->config.boxWidth, FPoint(corehandle->config.boxCenter)); - - if( corehandle->config.nbThreads != 0){ - omp_set_num_threads(corehandle->config.nbThreads); - } - - for(FSize idxPart = 0 ; idxPart < (*nb) ; ++idxPart){ - const FReal* pos = &position[idxPart * 3]; - corehandle->octree->insert(FPoint(pos[0], pos[1], pos[2]),idxPart); - } - - return FMMAPI_NO_ERROR; -} /* transmet au FmmCore les potentiels associés aux points sources. -Le tableau potential est alloué et libéré par la routine appelant, le FmmCore doit donc en faire une copie.*/ - - -int FmmCore_setPotentials(void *fmmCore, void *potentials) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - OctreeClass* octree = corehandle->octree; - - void* fmmKernel; - FmmCore_getKernelData(fmmCore, &fmmKernel); - - int sizeOfPotentials; - FmmKernel_getParameter(fmmKernel, FMMKERNEL_POTENTIAL_DATA_SIZE, &sizeOfPotentials); - - typename OctreeClass::Iterator octreeIterator(octree); - octreeIterator.gotoBottomLeft(); - do{ - CoreContainerClass* particles = octreeIterator.getCurrentLeaf()->getSrc(); - particles->allocatePotentials(sizeOfPotentials); - - for(int idx = 0 ; idx < particles->getNbParticles() ; ++idx){ - memcpy(((char*)particles->getPotentials())+idx*sizeOfPotentials, - ((char*)potentials)+particles->getIndexes()[idx]*sizeOfPotentials, - sizeOfPotentials); - } - - } while( octreeIterator.moveRight() ); - - return FMMAPI_NO_ERROR; -} /* transmet au FmmCore les potentiels associés aux points sources. -Le tableau potential est alloué et libéré par la routine appelant, le FmmCore doit donc en faire une copie.*/ - -int FmmCore_doComputation(void *fmmCore) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - - { // Ceck if there is number of NbPart summed at level 1 - FPoint corner = corehandle->octree->getBoxCenter(); - corner -= (corehandle->octree->getBoxWidth()/2); - - FReal boxWidth = corehandle->octree->getBoxWidth()/( 1<< (corehandle->config.treeHeight-1) ); - - typename OctreeClass::Iterator octreeIterator(corehandle->octree); - octreeIterator.gotoBottomLeft(); - for(int idxLevel = corehandle->config.treeHeight - 1 ; idxLevel >= 1 ; --idxLevel ){ - int multipoleSize; - int localSize; - - FmmKernel_getMultipoleArraySize(fmmCore, &multipoleSize); - FmmKernel_getLocalArraySize(fmmCore, &localSize); - - do{ - octreeIterator.getCurrentCell()->createArrays(multipoleSize, localSize); - octreeIterator.getCurrentCell()->setLevel(idxLevel); - const FTreeCoordinate coord = octreeIterator.getCurrentCell()->getCoordinate(); - FPoint position( coord.getX()*boxWidth + boxWidth/2.0 + corner.getX(), - coord.getY()*boxWidth + boxWidth/2.0 + corner.getY(), - coord.getZ()*boxWidth + boxWidth/2.0 + corner.getZ()); - octreeIterator.getCurrentCell()->setPosition(position.getDataValue()); - } while( octreeIterator.moveRight() ); - - octreeIterator.moveUp(); - octreeIterator.gotoLeft(); - - boxWidth *= FReal(2.0); - } - } - { // Ceck if there is number of NbPart summed at level 1 - void* fmmKernel; - FmmCore_getKernelData(fmmCore, &fmmKernel); - - int sizeOfFields; - FmmKernel_getParameter(fmmKernel, FMMKERNEL_FIELD_DATA_SIZE, &sizeOfFields); - - typename OctreeClass::Iterator octreeIterator(corehandle->octree); - octreeIterator.gotoBottomLeft(); - do{ - octreeIterator.getCurrentLeaf()->getTargets()->setParentCell(octreeIterator.getCurrentCell()); - octreeIterator.getCurrentLeaf()->getTargets()->allocateFields(sizeOfFields); - } while( octreeIterator.moveRight() ); - } - - if( corehandle->config.nbThreads <= 1){ - CoreKernelClass kernels(fmmCore ); - FmmClass algo(corehandle->octree,&kernels); - algo.execute(); - } - else{ - CoreKernelClass kernels(fmmCore ); - FmmClassThread algo(corehandle->octree,&kernels); - algo.execute(); - } - - return FMMAPI_NO_ERROR; -} /* réalise le produit multipôle. */ - -/* !!! Warning use *filed and not **field */ -int FmmCore_getField(void *fmmCore, void *fields) { - ScalFmmCoreHandle* corehandle = (ScalFmmCoreHandle*)fmmCore; - - typename OctreeClass::Iterator octreeIterator(corehandle->octree); - octreeIterator.gotoBottomLeft(); - do{ - CoreContainerClass* particles = octreeIterator.getCurrentLeaf()->getTargets(); - for(int idx = 0 ; idx < particles->getNbParticles() ; ++idx){ - memcpy(((char*)fields)+particles->getIndexes()[idx]*particles->getSizeOfField(), - ((char*)particles->getFields())+idx*particles->getSizeOfField(), - particles->getSizeOfField()); - } - } while(octreeIterator.moveRight()); - - return FMMAPI_NO_ERROR; -} /* récupère après le produit multipôle la valeur des champs en chaque point. -Le tableau field doit être alloué et libéré par la routine appelante. */ - - diff --git a/Addons/FmmApi/Src/ScalfmmApiRotation.cpp b/Addons/FmmApi/Src/ScalfmmApiRotation.cpp deleted file mode 100644 index b6a55a85a2490575cc5d158e752d35f24cb3d6c0..0000000000000000000000000000000000000000 --- a/Addons/FmmApi/Src/ScalfmmApiRotation.cpp +++ /dev/null @@ -1,508 +0,0 @@ -// =================================================================================== -// Logiciel initial: ScalFmm Version 0.5 -// Co-auteurs : Olivier Coulaud, Bérenger Bramas. -// Propriétaires : INRIA. -// Copyright © 2011-2012, diffusé sous les termes et conditions d’une licence propriétaire. -// Initial software: ScalFmm Version 0.5 -// Co-authors: Olivier Coulaud, Bérenger Bramas. -// Owners: INRIA. -// Copyright © 2011-2012, spread under the terms and conditions of a proprietary license. -// =================================================================================== -#include "../../Src/Utils/FGlobal.hpp" -#include "../../Src/Utils/FAssert.hpp" - -#include "../../Src/Containers/FVector.hpp" - -#include "../../Src/Components/FSimpleLeaf.hpp" - -#include "../../Src/Utils/FPoint.hpp" - -#include "../../Src/Kernels/Rotation/FRotationCell.hpp" -#include "../../Src/Kernels/Rotation/FRotationKernel.hpp" -#include "../../Src/Kernels/P2P/FP2PParticleContainer.hpp" - -#include "FmmApi.h" - -////////////////////// Opérateurs FMM Kernel : ////////////////////////// - -class KernelCell : public FBasicCell { - FComplex* multipole; - FComplex* local; -public: - KernelCell() : multipole(nullptr), local(nullptr){ - } - void attachArrays(FComplex inMultipole[], FComplex inLocal[]){ - multipole = inMultipole; - local = inLocal; - } - - const FComplex* getMultipole() const{ - return multipole; - } - - const FComplex* getLocal() const{ - return local; - } - - FComplex* getMultipole(){ - return multipole; - } - - FComplex* getLocal(){ - return local; - } -}; - - -static const int P = 5; - -typedef KernelCell KernelCellClass; -typedef FP2PParticleContainer ContainerClass; -typedef FSimpleLeaf LeafClass; -typedef FRotationKernel< KernelCellClass, ContainerClass , P> KernelClass; - -struct ScalFmmKernelHandle { - KernelClass** kernel; - int potentialDataSize; - int fieldDataSize; - int nbthread; -}; - -/******* Allocation : ******/ -int FmmKernel_init(void *fmmCore, void **fmmKernel){ - ScalFmmKernelHandle* kernelhandle = new ScalFmmKernelHandle; - memset(kernelhandle, 0, sizeof(ScalFmmKernelHandle)); - - int NbLevels; - FmmCore_getParameter(fmmCore, FMMCORE_TREE_HEIGHT, &NbLevels); - FReal boxWidth; - FmmCore_getParameter(fmmCore, FMMCORE_ROOT_BOX_WIDTH, &boxWidth); - FReal centerOfBox[3]; - FmmCore_getParameter(fmmCore, FMMCORE_ROOT_BOX_CENTER, centerOfBox); - - FmmCore_getParameter(fmmCore, FMMCORE_THREADS_NUMBER, &kernelhandle->nbthread); - - KernelClass original( NbLevels, boxWidth, FPoint(centerOfBox) ); - kernelhandle->kernel = new KernelClass*[kernelhandle->nbthread]; - for(int idxThread = 0 ; idxThread < kernelhandle->nbthread ; ++idxThread){ - kernelhandle->kernel[idxThread] = new KernelClass(original); - } - - kernelhandle->potentialDataSize = 1; - kernelhandle->fieldDataSize = 4; - - *fmmKernel = kernelhandle; - - return FMMAPI_NO_ERROR; -}/* : alloue et initialise le FmmKernel */ -int FmmKernel_free(void *fmmKernel){ - ScalFmmKernelHandle* kernelhandle = (ScalFmmKernelHandle*) fmmKernel; - - for(int idxThread = 0 ; idxThread < kernelhandle->nbthread ; ++idxThread){ - delete kernelhandle->kernel[idxThread]; - } - - delete[] kernelhandle->kernel; - delete kernelhandle; - - return FMMAPI_NO_ERROR; -} /* libére le FmmKernel */ - - - - -/******* Configuration : ***/ - -int FmmKernel_isParameterUsed(void * /*fmm*/, int *name, int *flag){ - switch( *name ){ - case FMMKERNEL_POTENTIAL_DATA_SIZE: - case FMMKERNEL_FIELD_DATA_SIZE: - case FMMKERNEL_HANDLES_P2P: - *flag = FMMAPI_SUPPORTED_PARAMETER; - break; - case FMMKERNEL_ACCURACY : - *flag = FMMAPI_UNSUPPORTED_PARAMETER; - break; - default: - *flag = FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_setParameter(void *fmmKernel, int *name, void*value){ - /*ScalFmmKernelHandle* kernelhandle = (ScalFmmKernelHandle*) fmmKernel;*/ - int flag; - - FmmKernel_isParameterUsed(fmmKernel, name, &flag); - if( flag != FMMAPI_SUPPORTED_PARAMETER){ - return flag; - } - - switch( *name ){ - case FMMKERNEL_POTENTIAL_DATA_SIZE : - case FMMKERNEL_FIELD_DATA_SIZE : - return FMMAPI_SUPPORTED_PARAMETER; - default: - return FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_setParameter(void *fmmKernel, int name, void*value){ - return FmmKernel_setParameter( fmmKernel, &name, value); -} - -int FmmKernel_getParameter(void *fmmKernel, int *name, void*value){ - ScalFmmKernelHandle* kernelhandle = (ScalFmmKernelHandle*) fmmKernel; - int flag; - - FmmKernel_isParameterUsed(fmmKernel, name, &flag); - if( flag != FMMAPI_SUPPORTED_PARAMETER){ - return flag; - } - - switch( *name ){ - case FMMKERNEL_POTENTIAL_DATA_SIZE : - *(int*)value = kernelhandle->potentialDataSize*sizeof(FReal); - break; - case FMMKERNEL_FIELD_DATA_SIZE : - *(int*)value = kernelhandle->fieldDataSize*sizeof(FReal); - break; - default: - return FMMAPI_UNKNOWN_PARAMETER; - } - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_getParameter(void *fmmKernel, int name, void*value){ - return FmmKernel_getParameter(fmmKernel, &name, value); -} - - -/****** Données FMM : *****/ -int FmmKernel_getMultipoleArraySize(void */*fmmCore*/, int *size) { - *size = ((P+2)*(P+1))/2 * sizeof(FComplex); - return FMMAPI_NO_ERROR; -} /* Renvoie dans size la taille (en octets) de l'expansion multipôle associée à la boîte boxId */ - -int FmmKernel_getLocalArraySize(void */*fmmCore*/, int *size){ - *size = ((P+2)*(P+1))/2 * sizeof(FComplex); - return FMMAPI_NO_ERROR; -} /* Renvoie dans size la taille (en octets) de l'expansion locale associée à la boîte boxId*/ - - -/******* Opérateurs FMM : **/ -int FmmKernel_P2M(void *fmmCore, void* boxId){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FComplex* multipole; - FmmCore_getMultipoleArray(fmmCore, boxId, (void**)&multipole); - - KernelCellClass cell; - cell.attachArrays(multipole, nullptr); - int coord[3]; - FmmCore_getCoord(fmmCore, boxId, coord); - cell.setCoordinate(coord[0], coord[1], coord[2]); - - FReal* positions; - FReal* physicalValues; - int number; - FmmCore_getSource(fmmCore, boxId, &positions, (void**)&physicalValues, &number); - - FP2PParticleContainer sources; - for(FSize idxPart = 0 ; idxPart < number ; ++idxPart){ - sources.push(FPoint(positions[idxPart*3],positions[idxPart*3+1],positions[idxPart*3+2]),physicalValues[idxPart]); - } - - kernelhandle->kernel[threadId]->P2M(&cell, &sources); - - FmmCore_releaseSource(fmmCore, boxId, physicalValues, positions); - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_L2P(void *fmmCore, void* boxId){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FComplex* local; - FmmCore_getLocalArray(fmmCore, boxId, (void**)&local); - - KernelCellClass cell; - cell.attachArrays(nullptr,local); - int coord[3]; - FmmCore_getCoord(fmmCore, boxId, coord); - cell.setCoordinate(coord[0], coord[1], coord[2]); - - FReal* physicalValues; - FReal* positions; - int number; - FmmCore_getSource(fmmCore, boxId, &positions, (void**)&physicalValues, &number); - - FReal* fields = new FReal[number*kernelhandle->fieldDataSize]; - FmmCore_getTargetField(fmmCore, boxId, fields); - - FP2PParticleContainer targets; - for(FSize idxPart = 0 ; idxPart < number ; ++idxPart){ - targets.push(FPoint(&positions[idxPart*3]),physicalValues[idxPart], - fields[idxPart*kernelhandle->fieldDataSize], - fields[idxPart*kernelhandle->fieldDataSize+1], - fields[idxPart*kernelhandle->fieldDataSize+2], - fields[idxPart*kernelhandle->fieldDataSize+3]); - } - - kernelhandle->kernel[threadId]->L2P(&cell, &targets); - - const FReal*const potentials = targets.getPotentials(); - const FReal*const forcesX = targets.getForcesX(); - const FReal*const forcesY = targets.getForcesY(); - const FReal*const forcesZ = targets.getForcesZ(); - - for(FSize idxPart = 0 ; idxPart < number ; ++idxPart){ - fields[idxPart*kernelhandle->fieldDataSize] += potentials[idxPart]; - fields[idxPart*kernelhandle->fieldDataSize+1] += forcesX[idxPart]; - fields[idxPart*kernelhandle->fieldDataSize+2] += forcesY[idxPart]; - fields[idxPart*kernelhandle->fieldDataSize+3] += forcesZ[idxPart]; - } - - FmmCore_releaseTargetPoints(fmmCore, boxId, positions); - FmmCore_setTargetField(fmmCore, boxId, fields); - delete[] fields; - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_M2M(void *fmmCore, void *boxIdFather, void *boxIdSon){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FComplex* multipole; - FmmCore_getMultipoleArray(fmmCore, boxIdFather, (void**)&multipole); - - KernelCellClass cellFather; - cellFather.attachArrays(multipole, nullptr); - int coordFather[3]; - FmmCore_getCoord(fmmCore, boxIdFather, coordFather); - cellFather.setCoordinate(coordFather[0], coordFather[1], coordFather[2]); - - FmmCore_getMultipoleArray(fmmCore, boxIdSon, (void**)&multipole); - - KernelCellClass cellSon; - cellSon.attachArrays(multipole, nullptr); - int coordChild[3]; - FmmCore_getCoord(fmmCore, boxIdSon, coordChild); - cellSon.setCoordinate(coordChild[0], coordChild[1], coordChild[2]); - - int level; - FmmCore_getLevel(fmmCore,boxIdFather, &level); - - const KernelCellClass* children[8]; - memset(children, 0, sizeof(KernelCellClass*)*8); - const int mindex = ((coordChild[0]&1) * 2 + (coordChild[1]&1)) * 2 + (coordChild[2]&1); - children[mindex] = &cellSon; - - kernelhandle->kernel[threadId]->M2M(&cellFather, children, level); - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_L2L(void *fmmCore, void *boxIdFather, void *boxIdSon){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FComplex* local; - FmmCore_getLocalArray(fmmCore, boxIdFather, (void**)&local); - - KernelCellClass cellFather; - cellFather.attachArrays(nullptr, local); - int coordFather[3]; - FmmCore_getCoord(fmmCore, boxIdFather, coordFather); - cellFather.setCoordinate(coordFather[0], coordFather[1], coordFather[2]); - - FmmCore_getLocalArray(fmmCore, boxIdSon, (void**)&local); - - KernelCellClass cellSon; - cellSon.attachArrays(nullptr, local); - int coordChild[3]; - FmmCore_getCoord(fmmCore, boxIdSon, coordChild); - cellSon.setCoordinate(coordChild[0], coordChild[1], coordChild[2]); - - int level; - FmmCore_getLevel(fmmCore,boxIdFather, &level); - - KernelCellClass* children[8]; - memset(children, 0, sizeof(KernelCellClass*)*8); - const int mindex = ((coordChild[0]&1) * 2 + (coordChild[1]&1)) * 2 + (coordChild[2]&1); - children[mindex] = &cellSon; - - kernelhandle->kernel[threadId]->L2L(&cellFather, children, level); - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_M2L(void *fmmCore, void *boxIdSrc, void *boxIdDest){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FComplex* multipole; - FmmCore_getMultipoleArray(fmmCore, boxIdSrc, (void**)&multipole); - KernelCellClass cellSrc; - cellSrc.attachArrays(multipole,nullptr); - int coord[3]; - FmmCore_getCoord(fmmCore, boxIdSrc, coord); - cellSrc.setCoordinate(coord[0], coord[1], coord[2]); - - FComplex* local; - FmmCore_getLocalArray(fmmCore, boxIdDest, (void**)&local); - KernelCellClass cellDst; - cellDst.attachArrays(nullptr, local); - FmmCore_getCoord(fmmCore, boxIdDest, coord); - cellDst.setCoordinate(coord[0], coord[1], coord[2]); - - int level; - FmmCore_getLevel(fmmCore, boxIdDest, &level); - - const int xdiff = cellSrc.getCoordinate().getX() - cellDst.getCoordinate().getX(); - const int ydiff = cellSrc.getCoordinate().getY() - cellDst.getCoordinate().getY(); - const int zdiff = cellSrc.getCoordinate().getZ() - cellDst.getCoordinate().getZ(); - const int index = (((xdiff+3) * 7) + (ydiff+3)) * 7 + zdiff + 3; - - const KernelCellClass* inter[343]; - memset(inter, 0, sizeof(KernelCellClass*)*343); - inter[index] = &cellSrc; - - kernelhandle->kernel[threadId]->M2L(&cellDst, inter, 1, level); - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_P2P_inner(void *fmmCore, void *boxIdSrcDest){ - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FReal* positionsTargets; - FReal* potentialsTargets; - int numberTargets; - FmmCore_getSource(fmmCore, boxIdSrcDest, &positionsTargets, (void**)&potentialsTargets, &numberTargets); - - FReal* fieldsTargets = new FReal[numberTargets*kernelhandle->fieldDataSize]; - FmmCore_getTargetField(fmmCore, boxIdSrcDest, fieldsTargets); - - int coordTargets[3]; - FmmCore_getCoord(fmmCore, boxIdSrcDest, coordTargets); - FTreeCoordinate treecoord(coordTargets[0], coordTargets[1], coordTargets[2]); - - FP2PParticleContainer targets; - for(FSize idxPart = 0 ; idxPart < numberTargets ; ++idxPart){ - targets.push(FPoint(positionsTargets[idxPart*3],positionsTargets[idxPart*3+1],positionsTargets[idxPart*3+2]), - potentialsTargets[idxPart], - fieldsTargets[idxPart*kernelhandle->fieldDataSize], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+1], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+2], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+3]); - } - - - FP2PParticleContainer* ptrs[27]; - memset(ptrs, 0, sizeof(FP2PParticleContainer*) * 27); - kernelhandle->kernel[threadId]->P2P(treecoord, &targets, &targets, ptrs, 0); - - const FReal*const potentials = targets.getPotentials(); - const FReal*const forcesX = targets.getForcesX(); - const FReal*const forcesY = targets.getForcesY(); - const FReal*const forcesZ = targets.getForcesZ(); - - for(FSize idxPart = 0 ; idxPart < numberTargets ; ++idxPart){ - fieldsTargets[idxPart*kernelhandle->fieldDataSize] = potentials[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+1] = forcesX[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+2] = forcesY[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+3] = forcesZ[idxPart]; - } - - FmmCore_releaseTargetPoints(fmmCore, boxIdSrcDest, positionsTargets); - FmmCore_setTargetField(fmmCore, boxIdSrcDest, fieldsTargets); - delete[] fieldsTargets; - - return FMMAPI_NO_ERROR; -} - -int FmmKernel_P2P(void *fmmCore, void *boxIdSrc, void *boxIdDest){ //return FMMAPI_NO_ERROR; - ScalFmmKernelHandle* kernelhandle; - FmmCore_getKernelData(fmmCore, (void**)&kernelhandle); - int threadId; - FmmCore_getParameter(fmmCore, FMMCORE_THREAD_ID, &threadId); - - FReal* positionsTargets; - FReal* potentialsTargets; - int numberTargets; - FmmCore_getSource(fmmCore, boxIdDest, &positionsTargets, (void**)&potentialsTargets, &numberTargets); - - FReal* fieldsTargets = new FReal[numberTargets*kernelhandle->fieldDataSize]; - FmmCore_getTargetField(fmmCore, boxIdDest, fieldsTargets); - - int coordTargets[3]; - FmmCore_getCoord(fmmCore, boxIdDest, coordTargets); - FTreeCoordinate treecoord(coordTargets[0], coordTargets[1], coordTargets[2]); - - FP2PParticleContainer targets; - for(FSize idxPart = 0 ; idxPart < numberTargets ; ++idxPart){ - targets.push(FPoint(positionsTargets[idxPart*3],positionsTargets[idxPart*3+1],positionsTargets[idxPart*3+2]), - potentialsTargets[idxPart], - fieldsTargets[idxPart*kernelhandle->fieldDataSize], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+1], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+2], - fieldsTargets[idxPart*kernelhandle->fieldDataSize+3]); - } - - FReal* positionsSrc; - FReal* potentialsSrc; - int numberSources; - FmmCore_getSource(fmmCore, boxIdSrc, &positionsSrc, (void**)&potentialsSrc, &numberSources); - - FP2PParticleContainer sources; - for(FSize idxPart = 0 ; idxPart < numberSources ; ++idxPart){ - sources.push(FPoint(positionsSrc[idxPart*3],positionsSrc[idxPart*3+1],positionsSrc[idxPart*3+2]),potentialsSrc[idxPart]); - } - - FP2PParticleContainer* ptrs[27]; - memset(ptrs, 0, sizeof(FP2PParticleContainer*) * 27); - ptrs[0] = &sources; - kernelhandle->kernel[threadId]->P2PRemote(treecoord, &targets, &targets, ptrs, 1); - - const FReal*const potentials = targets.getPotentials(); - const FReal*const forcesX = targets.getForcesX(); - const FReal*const forcesY = targets.getForcesY(); - const FReal*const forcesZ = targets.getForcesZ(); - - for(FSize idxPart = 0 ; idxPart < numberTargets ; ++idxPart){ - fieldsTargets[idxPart*kernelhandle->fieldDataSize] = potentials[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+1] = forcesX[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+2] = forcesY[idxPart]; - fieldsTargets[idxPart*kernelhandle->fieldDataSize+3] = forcesZ[idxPart]; - } - - FmmCore_releaseSource(fmmCore, boxIdDest, potentialsSrc, positionsSrc); - FmmCore_releaseTargetPoints(fmmCore, boxIdDest, positionsTargets); - FmmCore_setTargetField(fmmCore, boxIdDest, fieldsTargets); - delete[] fieldsTargets; - - return FMMAPI_NO_ERROR; -} /* pas mutuel, i.e. on fait seulement dans 1 sens. */ - - diff --git a/Addons/FmmApi/Tests/testFmmApi.cpp b/Addons/FmmApi/Tests/testFmmApi.cpp deleted file mode 100644 index ae699f5d368873395da7d38c4763d5546535cafd..0000000000000000000000000000000000000000 --- a/Addons/FmmApi/Tests/testFmmApi.cpp +++ /dev/null @@ -1,193 +0,0 @@ -// =================================================================================== -// Logiciel initial: ScalFmm Version 0.5 -// Co-auteurs : Olivier Coulaud, Bérenger Bramas. -// Propriétaires : INRIA. -// Copyright © 2011-2012, diffusé sous les termes et conditions d’une licence propriétaire. -// Initial software: ScalFmm Version 0.5 -// Co-authors: Olivier Coulaud, Bérenger Bramas. -// Owners: INRIA. -// Copyright © 2011-2012, spread under the terms and conditions of a proprietary license. -// =================================================================================== - -#include -#include - -#include "../../Src/Utils/FParameters.hpp" -#include "../../Src/Utils/FTic.hpp" -#include "../../Src/Utils/FMemUtils.hpp" - -#include "../../Src/Files/FRandomLoader.hpp" - -#include "../Src/FmmApi.h" - -/** This program show an example of use of the fmm api - */ - -// Simply create particles and try the kernels -int main(int argc, char ** argv){ - ///////////////////////What we do///////////////////////////// - std::cout << ">> This executable has to be used to test the FMM API\n"; - ////////////////////////////////////////////////////////////// - - int NbLevels = FParameters::getValue(argc,argv,"-h", 7); - int SizeSubLevels = FParameters::getValue(argc,argv,"-sh", 3); - const int NbPart = FParameters::getValue(argc,argv,"-nb", 2000000); - FTic counter; - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - FRandomLoader loader(NbPart, 1, FPoint(0.5,0.5,0.5), 1); - - void* FmmCoreHandle; - FmmCore_init(&FmmCoreHandle); - - FmmCore_setParameter(FmmCoreHandle, FMMCORE_TREE_HEIGHT, &NbLevels); - FReal boxWidth = loader.getBoxWidth(); - FmmCore_setParameter(FmmCoreHandle, FMMCORE_ROOT_BOX_WIDTH, &boxWidth); - FmmCore_setParameter(FmmCoreHandle, FMMCORE_ROOT_BOX_CENTER, loader.getCenterOfBox().getDataValue()); - - void* FmmKernelHandle; - FmmKernel_init(FmmCoreHandle, &FmmKernelHandle); - - FmmCore_setKernelData(FmmCoreHandle, FmmKernelHandle); - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - std::cout << "Creating & Inserting " << NbPart << " particles ..." << std::endl; - std::cout << "\tHeight : " << NbLevels << " \t sub-height : " << SizeSubLevels << std::endl; - counter.tic(); - - FSize r.getNumberOfParticles(); - FReal* potentials = new FReal[nbPart]; - FReal* positions = new FReal[nbPart*3]; - - { - FPoint part; - const FReal physicalValue = 0.1; - - for(int idx = 0 ; idx < nbPart ; ++idx){ - loader.fillParticle(&part); - - potentials[idx] = physicalValue; - positions[3*idx] = part.getX(); - positions[3*idx+1] = part.getY(); - positions[3*idx+2] = part.getZ(); - } - - FmmCore_setPositions(FmmCoreHandle, &nbPart, positions); - FmmCore_setPotentials(FmmCoreHandle, potentials); - } - - counter.tac(); - std::cout << "Done " << "(@Creating and Inserting Particles = " << counter.elapsed() << "s)." << std::endl; - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - std::cout << "Working on particles ..." << std::endl; - counter.tic(); - - FmmCore_doComputation(FmmCoreHandle); - - counter.tac(); - std::cout << "Done " << "(@Algorithm = " << counter.elapsed() << "s)." << std::endl; - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - std::cout << "Getting results ..." << std::endl; - counter.tic(); - - FReal* fields = new FReal[4*nbPart]; - - FmmCore_getField(FmmCoreHandle, fields); - - counter.tac(); - std::cout << "Done " << "(@retrieve = " << counter.elapsed() << "s)." << std::endl; - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - std::cout << "Computing direct ..." << std::endl; - counter.tic(); - - FReal* fieldsdirect = new FReal[4*nbPart]; - FMemUtils::setall(fieldsdirect, 0.0, nbPart * 4); - - for(int idxTarget = 0 ; idxTarget < nbPart ; ++idxTarget){ - - for(int idxSource = idxTarget + 1 ; idxSource < nbPart ; ++idxSource){ - FReal dx = positions[idxSource*3] - positions[idxTarget*3]; - FReal dy = positions[idxSource*3+1] - positions[idxTarget*3+1]; - FReal dz = positions[idxSource*3+2] - positions[idxTarget*3+2]; - - FReal inv_square_distance = FReal(1.0) / (dx*dx + dy*dy + dz*dz); - FReal inv_distance = FMath::Sqrt(inv_square_distance); - - inv_square_distance *= inv_distance; - inv_square_distance *= potentials[idxTarget] * potentials[idxSource]; - - dx *= inv_square_distance; - dy *= inv_square_distance; - dz *= inv_square_distance; - - fieldsdirect[4*idxTarget+1] += dx; - fieldsdirect[4*idxTarget+2] += dy; - fieldsdirect[4*idxTarget+3] += dz; - fieldsdirect[4*idxTarget+0] += inv_distance * potentials[idxSource]; - - fieldsdirect[4*idxSource+1] += -dx; - fieldsdirect[4*idxSource+2] += -dy; - fieldsdirect[4*idxSource+3] += -dz; - fieldsdirect[4*idxSource+0] += inv_distance * potentials[idxTarget]; - - } - } - - - counter.tac(); - std::cout << "Done " << "(@direct = " << counter.elapsed() << "s)." << std::endl; - - ////////////////////////////////////////////////////////////////////////////////// - ////////////////////////////////////////////////////////////////////////////////// - - - std::cout << "Comparing results ..." << std::endl; - counter.tic(); - - FMath::FAccurater forces; - FMath::FAccurater potential; - - for(int idx = 0 ; idx < nbPart ; ++idx){ - if(idx < 10){ - std::cout << fieldsdirect[4*idx] << " fmm " << fields[4*idx] << std::endl; - } - forces.add( fieldsdirect[4*idx] ,fields[4*idx]); - forces.add( fieldsdirect[4*idx+1] ,fields[4*idx+1]); - forces.add( fieldsdirect[4*idx+2] ,fields[4*idx+2]); - potential.add( fieldsdirect[4*idx+3] ,fields[4*idx+3]); - } - - counter.tac(); - std::cout << "Done " << "(@comparing = " << counter.elapsed() << "s)." << std::endl; - - std::cout << "\tForces inf " << forces.getInfNorm() << " normL2 " << forces.getL2Norm() << std::endl; - std::cout << "\tPotential inf " << potential.getInfNorm() << " normL2 " << potential.getL2Norm() << std::endl; - - ////////////////////////////////////////////////////////////////////////////////// - - delete[] fieldsdirect; - delete[] fields; - delete[] potentials; - delete[] positions; - - FmmCore_free(FmmCoreHandle); - FmmKernel_free(FmmKernelHandle); - - return 0; -} - - -