Freefem++ Home Page (7-november-2005)

FreeFem++ version 3.9-2 (2d and 3d)

Université Pierre et Marie Curie Laboratoire Jacques-Louis Lions

Financé par l'ANR / FF2A3, Réf. : ANR-07-CIS7-002-01

FreeFem++ team : Olivier Pironneau, Frédéric Hecht, Antoine Le Hyaric, Jacques Morice ,

Please send questions, comments and bug reports to Frédéric Hecht or to the FreeFem++ mailing list.

The second tutorial and Workshop on FreeFem++ , held September 1st and 2nd, 2010, in Paris at Université Pierre et Marie Curie, Barre 16-15, 3ieme, 4 place Jussieu, Paris

Introduction

FreeFem++ is an implementation of a language dedicated to the finite element method. It enables you to solve Partial Differential Equations (PDE) easily.
Problems involving PDE (2d, 3d) from several branches of physics such as fluid-structure interactions require interpolations of data on several meshes and their manipulation within one program. FreeFem++ includes a fast 2^d-tree-based interpolation algorithm and a language for the manipulation of data on multiple meshes (as a follow up of bamg).
FreeFem++ is written in C++ and the FreeFem++ language is a C++ idiom. It runs on any Unix-like OS (with g++ version 3 or higher, X11R6 or OpenGL with GLUT) Linux, FreeBSD, Solaris 10, Microsoft Windows (95, 98, 2000, NT, XP, Vista) and MacOS X (native version using OpenGL). FreeFem++ replaces the older freefem and freefem+.

Quick Start

The FreeFem++-cs (IDE: an integrated development environment) by Antoine Le Hyaric is a excellent tool to develop FreeFem++ scripts.

Unix Linux, FreeBSD, NetBSD, Solaris 10, ...	Windows 98, NT, 2000, XP	MacOS X Universal ( 10.3 or + )	Documentation	Online Demonstration

Download this sources archive, decompress, compile; , install and run	Download this EXE and run it	Download Universal (Intel, G3), decompress, install (see INSTALL-MacOSX) and run	Full documentation : freefem++doc.pdf	Try FreeFem++ online at : LJLL Software Online

I have built a precompiled version on ubuntu 32bits, to install do under root in /usr/local (with openmpi )
tar zxvf freefem++-v3.5-1-ubuntu.tar.gz -C /
If some lib is missing, use the -k option (existing files are kept missing ones are installed in /usr/lib):
tar zxvf freefem++-v3.2-usr-lib.tar.gz -C / -k

Bug under UBUNTU/GLUT: missing graphics window title , problem with compiz package and freeglut, you can try to kill compiz with killall compiz.real && killall compiz && metacity --replace --display :0 & or To disable/enable it graphically, install the compiz-swith package with synaptic.
To compile FreeFem++, go into the FreeFem++ directory and type : ./configure --enable-download make make install (as root).
For all other platforms and versions, please refer to the full list of downloads.

Documentation To receive the manual of FreeFem++ 2.21 in book form, send three blank DVDs and your postal address to : Dr FREEFEM++ Laboratoire JL Lions University of Paris VI 175 rue du Chevaleret 75013 Paris France A new reedition of the book is now available New Manual (freefem++doc.pdf) (8515.8 Kb Tue Aug 31 13:27:12 2010). This is a complete rewrite with lots of examples. INNOVATION : last news. HISTORY : life of the software from 2005. HISTORY_BEFORE_2005 : life of the software before 2005. Old obsolete but shorter version of the full manual still available as manual.pdf (0 Kb Thu Jan 1 01:00:00 1970). News 1/10/2005 Read or participate in general discussions about FreeFEM++ with the FreeFEM++ mailing list 23/9/2005 FreeFem++ and OpenFEM presentation day at IHP (Paris) : details, slides and examples 16/10/2008 Ateliers de Simualtion Numérique, (Setif,Algérie) : details, par Olivier Pantz The first tutorial and Workshop on FreeFem++ , held September 14th and 15th, 2009, in Paris at IHP, amphi Hermite ,11 Rue Pierre et Marie Curie. 1-2/09/2010 The second tutorial and Workshop on FreeFem++ , held September 1st and 2nd, 2010, in Paris at Université Pierre et Marie Curie, Barre 16-15, 3ieme, 4 place Jussieu, Paris Examples in 2d and 3d. A small movie (340Kb) : Cool air (green) comes from the lower left and mix with hot air (magenta), the right boundary is free. This is Navier-Stokes-Boussinesq integrated with P1-bubblle P1 mixte finite element. A very small example 2d of how to solve the Poisson equation on a L shape : border aaa(t=0,1){x=t;y=0;}; border bbb(t=0,0.5){x=1;y=t;}; border ccc(t=0,0.5){x=1-t;y=0.5;}; border ddd(t=0.5,1){x=0.5;y=t;}; border eee(t=0.5,1){x=1-t;y=1;}; border fff(t=0,1){x=0;y=1-t;}; mesh Th = buildmesh (aaa(6) + bbb(4) + ccc(4) +ddd(4) + eee(4) + fff(6)); fespace Vh(Th,P1); // to change P1 in P2 to make P2 finite element. Vh u=0,v; func f= 1; func g= 0; int i=0; real error=0.1, coef= 0.1^(1./5.); problem Probem1(u,v,solver=CG,eps=-1.0e-6) = int2d(Th)( dx(u)*dx(v) + dy(u)*dy(v)) + int2d(Th) ( v*f ) + on(aaa,bbb,ccc,ddd,eee,fff,u=g) ; for (i=0;i< 10;i++) { real d = clock(); Probem1; // solve the problem plot(u,Th,wait=1); Th=adaptmesh(Th,u,inquire=1,err=error); error = error * coef; } ; Solution on adapted mesh and associated mesh. A 3d example A very small example of how to solve the Stokes equation 3d on cube shape : load "msh3" load "medit" // dynamics load tools for 3d. int nn=8; mesh Th2=square(nn,nn); fespace Vh2(Th2,P2); Vh2 ux,uz,p2; int[int] rup=[0,2], rdown=[0,1], rmid=[1,1,2,1,3,1,4,1]; real zmin=0,zmax=1; mesh3 Th=buildlayers(Th2,nn, zbound=[zmin,zmax], reffacemid=rmid, reffaceup = rup, reffacelow = rdown); medit("c10x10x10",Th); // see the 3d mesh with medit software fespace VVh(Th,[P23d,P23d,P23d,P13d]); macro Grad(u) [dx(u),dy(u),dz(u)] // EOM macro div(u1,u2,u3) (dx(u1)+dy(u2)+dz(u3)) //EOM VVh [u1,u2,u3,p]; VVh [v1,v2,v3,q]; solve vStokes([u1,u2,u3,p],[v1,v2,v3,q]) = int3d(Th,qforder=3)( Grad(u1)'*Grad(v1) + Grad(u2)'*Grad(v2) + Grad(u3)'*Grad(v3) - div(u1,u2,u3)*q - div(v1,v2,v3)*p + 1e-10*q*p ) + on(2,u1=1.,u2=0,u3=0) + on(1,u1=0,u2=0,u3=0) ; plot(p,wait=1, nbiso=5); // a 3d plot of iso pressure. in progress... march 2009 // to see the 10 cut plan in 2d for(int i=1;i<10;i++) { real yy=i/10.; // compute yy. // do 3d -> 2d interpolation. ux= u1(x,yy,y); uz= u3(x,yy,y); p2= p(x,yy,y); plot([ux,uz],p2,cmm=" cut y = "+yy,wait= 1); } Solution on cup plan y=0.5 and mesh 10x10x10 and associated mesh. Ongoing Work 3D implementation: new solver, new mesh tools, new kind of finite element There will be some changes in the general policy for FreeFEM (GNU Public License, LGPL) Full List of Downloads The current version of Freefem++ is 3.9-2. You can get the latest source from an anonymous Mercurial SCM copy with the following unix shell commands : hg clone http://www.freefem.org/ff++/ff++ Self-contained archives for all other systems : Setup file for generic hardware and Windows 95,98,NT,2000, XP FreeFem++-3.9-02.exe 13058.5 Kb Tue Aug 31 15:12:15 2010 All the sources (tar+gzip) freefem++-3.9-2.tar.gz 21425.4 Kb Tue Aug 31 13:00:48 2010 All the sources (zip file) freefem++-3.5.zip 18730.3 Kb Mon Aug 31 11:39:04 2009 MacOs 10.3 generic PPC Cocoa OpenGL (tar+gzip) FreeFem++v2.22_MacOsX.tgz 8039.3 Kb Wed Nov 19 16:12:51 2008 MacOs 10.5 MacIntel Cocoa OpenGL (tar+gzip) FreeFem++v3.5-1-MacIntel_MacOsX.tgz 10286.7 Kb Fri Sep 25 11:26:32 2009 -Universal MacOsX (tar+ gzip file) FreeFem++v3.9-2-Universal_MacOsX.tgz 18951.7 Kb Tue Aug 31 13:25:05 2010 -usr-lib (tar+ gzip file) freefem++-v3.2-usr-lib.tar.gz 3900.6 Kb Thu May 14 21:49:33 2009 -ubuntu (tar+ gzip file) freefem++-v3.6-ubuntu.tar.gz 22350.7 Kb Thu Nov 5 13:56:09 2009 -snow-leopard MacOsX (tar+ gzip file) FreeFem++v3.9-2-snow-leopard_MacOsX.tgz 14746.4 Kb Tue Aug 31 13:26:58 2010 INNOVATION INNOVATION 31.1 Kb Tue Aug 31 13:27:11 2010 HISTORY HISTORY 64.8 Kb Wed Nov 19 16:12:44 2008 HISTORY BEFORE 2005 HISTORY_BEFORE_2005 22 Kb Wed Nov 19 16:13:01 2008 New manual (pdf file) freefem++doc.pdf 8515.8 Kb Tue Aug 31 13:27:12 2010 This directory contains all the different versions of FreeFem. Page written by Frédéric Hecht et Antoine Le Hyaric Last modified : Tue Aug 31 15:12:28 2010