~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to demo/ode/aliev-panfilov/cpp/main.cpp

  • Committer: Garth N. Wells
  • Date: 2008-03-29 09:34:25 UTC
  • Revision ID: gnw20@cam.ac.uk-20080329093425-3ea2vhjvccq56zvi
Add some basic build & install instructions.

Show diffs side-by-side

added added

removed removed

Lines of Context:
27
27
  
28
28
  void u0(uBlasVector& u)
29
29
  {
30
 
    u[0] = 0.2;
31
 
    u[1] = 0.0;
 
30
    u(0) = 0.2;
 
31
    u(1) = 0.0;
32
32
  }
33
33
  
34
34
  void f(const uBlasVector& u, real t, uBlasVector& y)
35
35
  {
36
 
    const real eps = eps0 + mu1*u[1] / (u[0] + mu2);
 
36
    const real eps = eps0 + mu1*u(1) / (u(0) + mu2);
37
37
 
38
 
    y[0] = -k*u[0]*(u[0] - a)*(u[0] - 1.0) - u[0]*u[1];
39
 
    y[1] = eps*(-u[1] - k*u[0]*(u[0] - a - 1.0));
 
38
    y(0) = -k*u(0)*(u(0) - a)*(u(0) - 1.0) - u(0)*u(1);
 
39
    y(1) = eps*(-u(1) - k*u(0)*(u(0) - a - 1.0));
40
40
  }
41
41
  
42
42
private: