integral_ex1.cpp File Reference

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

Go to the source code of this file.

Functions

void check_integral (ex &f)
int main ()


Function Documentation

void check_integral ( ex &  f  ) 

Definition at line 8 of file integral_ex1.cpp.

00008                           {
00009 }

int main (  ) 

Definition at line 11 of file integral_ex1.cpp.

References SyFi::Lagrange::compute_basis_functions(), SyFi::EQUAL_OR_DIE(), demos::crouzeixraviart::fe, SyFi::initSyFi(), demos::poisson1::integrand, SyFi::ReferenceLine::integrate(), SyFi::StandardFE::N(), SyFi::StandardFE::set_order(), and SyFi::StandardFE::set_polygon().

00011            {
00012 
00013     initSyFi(1); 
00014 
00015     ReferenceLine line; 
00016     Lagrange fe; 
00017     fe.set_polygon(line);
00018     fe.set_order(3); 
00019     fe.compute_basis_functions(); 
00020     ex integrand = fe.N(0)*fe.N(1);  
00021     cout <<integrand<<endl; 
00022     ex integral1 = line.integrate(integrand); 
00023     cout <<"Integral " <<integral1<<endl; 
00024     EQUAL_OR_DIE(integral1, "33/560"); 
00025 
00026 
00027 }


Generated on Mon Aug 31 16:16:49 2009 for SyFi by  doxygen 1.5.9