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

FreeFem++ version 3.7-1 (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, Jacques Morice .

with the help of Antoine Le Hyaric for the ide interface (freefem++-cs)

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

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

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 synatic.
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.
The FreeFem++-cs (IDE: an integrated development environment) by Antoine Le Hyaric is a excellent tool to develop FreeFem++ scripts.

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) (6826.9 Kb Wed Jan 13 23:42:01 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 14-15/09/2009 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. 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 mesjh 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, GPL) Full List of Downloads The current version of Freefem++ is 3.7-1. You can get the latest source from an anonymous CVS server with the following unix shell commands : cvs -d :pserver:anonymous@idared.ann.jussieu.fr:/Users/pubcvs/cvs login (password = freefem++) cvs -d :pserver:anonymous@idared.ann.jussieu.fr:/Users/pubcvs/cvs co freefem++ Self-contained archives for all other systems : Setup file for generic hardware and Windows 95,98,NT,2000, XP FreeFem++-3.7-1.exe 11143.3 Kb Thu Jan 14 16:33:52 2010 All the sources (tar+gzip) freefem++-3.7-1.tar.gz 18206.3 Kb Wed Jan 13 22:51:10 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.7-Universal_MacOsX.tgz 16772.2 Kb Tue Dec 8 13:33:19 2009 -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.7-1-snow-leopard_MacOsX.tgz 12147.1 Kb Wed Jan 13 22:57:10 2010 INNOVATION INNOVATION 28 Kb Wed Jan 13 23:48:41 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 6826.9 Kb Wed Jan 13 23:42:01 2010 This directory contains all the different versions of FreeFem. Page written by Frédéric Hecht et Antoine Le Hyaric Last modified : Thu Jan 14 16:34:04 2010