RD2

Bulletin:

2 dimensional simulation of the pattern on 2-variable Activator-Inhibitor Reaction Diffusion system [see rd2 on ALBERTA].

The CG solver of cusp library makes it faster.

* Use "double" for fast conversion of solver.

Specs:

  • Intel(R) Core(TM) i5-2400 CPU @ 3.10GHz, 16GB memory

deviceQuery,:

  • CUDA Driver = CUDART, CUDA Driver Version = 4.0, CUDA Runtime Version = 4.0, NumDevs = 1, Device = GeForce GTX 460

Compile options:

[single]

[double] -gencode arch=compute_20,code=sm_20

Average time for 1 step

  1. (0,1) x (0,1) [mesh:100x100, 1000step]

    1. (iteration limit)=10000, (relative tolerance)=1.0e-20

convergence of the iteration of the cg-method on cusp-library

    1. (0,10) x (0,10) [mesh:1000x1000, 20000step]

ε =1.0e-2, τ =0.1, dv=1.0, α =0.25, γ =10.0, vbias=0.0, tmax=20

ε =1.0e-2, τ =1.0, dv=1.0, α =0.125, γ =2.0, vbias=0.0, tmax=200

Acknowledgement :

Prof. T.Katagiri(u-tokyo) introduced me cusp-library. Thanks.

This work was supported by KAKENHI 00268666.