Subject: Re: [Syfi-dev] (no subject)

From: [email-masked]
Date: Thu, 5 Jul 2007 17:56:55 +0200 (CEST)

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