syfi team mailing list archive

[kent-and@xxxxxxxxx]

>
> Hello,
> My name is Mehdi. I am a French student
>
> Can you say me please how i can
>
>        *derivate fe.N(i) rwt x (first composant, for example)
>        *integrate on surface domaine
>
> Think you very much
>


In C++ you can do like this:


#include <ginac/ginac.h>
#include <iostream>
#include <SyFi/SyFi.h>

using namespace GiNaC;
using namespace SyFi;

int main () {

  initSyFi(3);
  ReferenceTetrahedron t;
  Lagrange fe(t,1);

  std::cout << "dN(0)/dx "<<diff(fe.N(0), x);
  std::cout << t.triangle(0).integrate(fe.N(0));
}