#include <SyFi.h>
#include <ginac/ginac.h>
Go to the source code of this file.
Functions | |
void | check_integral (ex &f) |
int | main () |
void check_integral | ( | ex & | f | ) |
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 }