Diffusion tamaas (maxwell_viscoelastic-branch) Browse

tamaas/
266740e73f8amaxwell_viscoelastic…

maxwell_viscoelastic-branch vs master

Commit	Author	Details	Committed
c22af0bd2267	Zichen	Initiation of maxwell_viscoelastic	May 14
43b439ffbc71	Zichen	Add maxwell_viscoelastic.cpp to source files	May 6
276abf143714	Zichen	Merge branch 'master' into maxwell_viscoelastic-branch	May 6
8d3562256686	Zichen	test git	Apr 29
101c853ffd8c	Zichen	test git	Apr 29

/

.gitignore
.gitmodules
COPYING
README.md
README.org
SConstruct
doc/
examples/
gpu/
python/
site_scons/
src/
tests/
version.py

README.md

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

You can customize a few compilation variables:

CXX changes the compiler
CXXFLAGS adds flags to the compilation process

For example:

scons CXX=icpc compiles with the intel C++ compiler
CXXFLAGS=-mavx2 scons adds the -mavx2 flag (which enables some Intel instructions) to the compilation process

If the CXX variable is defined in your environment, scons will take it as a default compiler.

Build the doxygen

scons build_doc=true

Coding

C++

All classes and functions in Tamaas live in the \tamaas\ namespace. Here is an example main:

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();
}

Python

Tamaas also has a python interface. Here is an example:

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()

Paralellism

Tamaas features shared-memory paralellism with OpenMP. The number of threads can be controlled via the OMP_NUM_THREADS environment variable or the omp_set_num_thread() function in the OpenMP API.

tamaas/266740e73f8amaxwell_viscoelastic…