MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

RE: FEM contourplotter - conversion completed

  • To: mathgroup at
  • Subject: [mg88244] RE: [mg88188] FEM contourplotter - conversion completed
  • From: "Jose Luis Gomez" <jose.luis.gomez at>
  • Date: Tue, 29 Apr 2008 06:50:39 -0400 (EDT)
  • References: <>

Dear Carlos

I have been using the great IMS (Imtek Mathematica Supplement) for Finite
Element Calculations in Mathematica: 

In case you already know IMS, could you please summarize in some rows what
would be the relationship between your work and IMS? Do you know if there
are calculations that can be done with your work and cannot be done with

Best regards!


-----Mensaje original-----
De: carlos at Colorado.EDU [mailto:carlos at Colorado.EDU] 
Enviado el: Domingo, 27 de Abril de 2008 06:01 a.m.
Para: mathgroup at
Asunto: [mg88188] FEM contourplotter - conversion completed

The translation of the Fortran IV 1966 finite element
contourplot to Mathematica is complete. Those interested in
that specialized topic may download a mini test version from

Click on ContBandPlotterMini.nb  link. Download and execute.
Self contained data. Tested so far under V4.1, 4.2 & 5.2.
Untested under V6.0 as my laptop in Spain lacks it.
Test plots take ~0.5 sec on 2.4 GHz MacBook Pro under 5.2.

To run under 6.0 set DisplayChannel to Print (see top comment).

The original Fortran code was roughly 9000 lines.  Code was
typical "write only" thesis spaghetti, with over 1000 GOTOs
and (believe it or not) even a PAUSE statement! Mathematica
code is about 800 lines. Imperative procedural logic reflects
source, but could simplify later conversion to C or assembly.

Posted mini is about 220 lines including test code. Size
was reduced by deleting 6 element types out of 7, and
removing an alternative divide & conquer contourplotting
method. (D&C produces higher quality plots but is slower.)

Code is slow (about 10K poly/sec) but that is enough for
coursework. A C or assembly version will be needed for FEM
production work (>=1M poly/sec) at the cost of portability.
Worst flaw is the shoddy contour value labeling, done in a
hurry.  To see what's wrong, try resizing top plot with
mouse.  Advice on how to fix that mess welcome.

  • Prev by Date: Re: Does Mathematica really need more printed, introductory documentation?
  • Next by Date: Re: Does Mathematica really need more printed, introductory documentation?
  • Previous by thread: Re: FEM contourplotter - conversion completed
  • Next by thread: how to use ones(A)?