tamaas/35db0e0d30c5archive/legacy_mpi
archive/legacy_mpi vs master
Commit | Author | Details | Committed | ||||
---|---|---|---|---|---|---|---|
b7fa37a3b7a3 | anciaux | code in place but not working | Dec 18 2017 | ||||
56c73d50b3ec | anciaux | test ok but not for dimension 1 | Dec 15 2017 | ||||
1c6e308391c0 | anciaux | fully working distribution and loop/reduce | Dec 15 2017 | ||||
3a32a987f380 | anciaux | renaming a file | Dec 15 2017 | ||||
d0e969ce80ae | anciaux | generalize test (not working at the moment) | Dec 15 2017 | ||||
9614f93afd40 | anciaux | reductions are working in mpi | Dec 14 2017 | ||||
5d4d5c48c3e0 | anciaux | adding the test for mpi loops and the gitignore | Dec 14 2017 | ||||
63f2a04e3716 | anciaux | first commit for the mpi version | Dec 14 2017 |
README
#+TITLE: Welcome to Tamaas, the Rough surface toolkit of LSMS
- Build
The build system uses SCons. In order to construct the library you should hit:
scons
And to speedup the process you can do:
scons -j 6
In order to clean the build
scons -c
In order to compile in debug
scons build_type=debug
Indeed the default was
scons build_type=release
In order to make the compilation more verbose
scons verbose=true
- Coding
- C++
All classes and functions in Tamaas live in the tamaas namespace. Here is an example main:
#+NAME: Example main #+BEGIN_SRC c++ using namespace tamaas; int main (int argc, char * argv[]) {
// Initializing tamaas initialize(); // Surface parameters const UInt grid_size = 1024; const UInt k0 = 4, k1 = 4, k2 = 32; const Real hurst = 0.8; const Real rms = 0.1; // Surface generator SurfaceGeneratorFilterFFT sg; sg.getGridSize() = grid_size; sg.getQ0() = k0; sg.getQ1() = k1; sg.getQ2() = k2; sg.getRMS() = rms; sg.getHurst() = hurst; sg.getRandomSeed() = 56; sg.Init(); Surface<Real> & surface = sg.buildSurface(); // BEM object BemPolonski bem(surface); bem.setEffectiveModulus(1.); // Solving equilibrium Real load = 0.01; bem.computeEquilibrium(1e-12, load); Surface<Real> & tractions = bem.getTractions(); Surface<Real> & displacements = bem.getDisplacements(); // Cleaning up tamaas finalize();
} #+END_SRC
- Python
Tamaas also has a python interface. Here is an example:
#+NAME: Python example #+BEGIN_SRC python import numpy as np import tamaas as tm
tm.initialize()
Surface generation
SG = tm.SurfaceGeneratorFilterFFT() SG.getGridSize().assign(512) SG.getHurst().assign(0.8) SG.getRMS().assign(1.); SG.getQ0().assign(4); SG.getQ1().assign(4); SG.getQ2().assign(32); SG.getRandomSeed().assign(156); SG.Init() surface = SG.buildSurface()
bem = tm.BemPolonski(surface) bem.setEffectiveModulus(1.)
load = 0.01 bem.computeEquilibrium(1e-12, load)
tractions = bem.getTractions() displacements = bem.getDisplacements()
tm.finalize() #+END_SRC
- Paralellism
Tamaas features shared-memory paralellism with OpenMP. The number of threads can be controlled via the OMP_NUM_THREADS environement variable or the omp_set_num_thread() function in the OpenMP API.