Revision 1488 trunk/src/df1b2-separable/df3fun.cpp

df3fun.cpp (revision 1488)
378 378
                   + ::exp(*x.get_u())* *x.get_udot2();
379 379

  
380 380
    *z.get_udot3() = ::exp(*x.get_u()) * cube(*x.get_udot())
381
                   + 3.0 * ::exp(*x.get_u()) * *x.get_udot() * *x.get_udot2();
381
                   + 3.0 * ::exp(*x.get_u()) * *x.get_udot() * *x.get_udot2()
382 382
                   + ::exp(*x.get_u()) * *x.get_udot3();
383 383
    return z;
384 384
  }
......
399 399
                   + xp * *x.get_udot2();
400 400

  
401 401
    *z.get_udot3() = xp3 * cube(*x.get_udot())
402
                   + 3.0 * xp2 * *x.get_udot() * *x.get_udot2();
402
                   + 3.0 * xp2 * *x.get_udot() * *x.get_udot2()
403 403
                   + xp * *x.get_udot3();
404 404
    return z;
405 405
  }
......
779 779
  if (value(x)>=0)
780 780
  {
781 781
    df3_one_variable u1=p*x;
782
    df3_one_variable u2=1.+p*x;
783
    df3_one_variable u=1./(1.+p*x);
782
    df3_one_variable u2=1.+u1;
783
    df3_one_variable u=1./u2;
784 784
    df3_one_variable y=  ((((b5*u+b4)*u+b3)*u+b2)*u+b1)*u;
785 785
    df3_one_variable tmp1=-0.3989422804*exp(-.5*x*x);
786 786
    z=1.0+tmp1*y;
......
811 811
  if (value(x)>=0)
812 812
  {
813 813
    df3_one_variable u1=p*x;
814
    df3_one_variable u2=1.+p*x;
815
    df3_one_variable u=1./(1.+p*x);
814
    df3_one_variable u2=1.+u1;
815
    df3_one_variable u=1./u2;
816 816
    df3_one_variable y=  ((((b5*u+b4)*u+b3)*u+b2)*u+b1)*u;
817 817
    df3_one_variable tmp1=-0.3989422804*exp(-.5*x*x);
818 818
    z=1.0+tmp1*y;

Also available in: Unified diff