aeslaughter / mfem Goto Github PK

Locating the element is done in the get_element, should this be another function?
The algorithm used to locate the element should be improved.

Create an overloaded version of point_value and point_gradient that search for the element, then call the element based version.

Compatible to free matlab clone octave?
No license is needed for octave.

If yes, please extend your info.
if unknown, untested with octave, info about this also nice.
if no, it is also ok. Some commands are not compatible in octave.

So you are the expert. Perhaps you can make a fast test run.
Here some last infos about octave.
https://wiki.octave.org/Differences_between_Octave_and_Matlab
https://wiki.octave.org/Release_History
https://octave.org/NEWS-7.html

Then edit example problems that use the gradient to use these functions.

Currently all boundaries and subdomains are labeled with numbers, this should be changed to accept text or number labels. Also, the labels left, right, ... should be automatically created thus available directly in the System.add_ members.

The FEmesh may be better suited as two classes, FE and Mesh. The FEmesh class could still exist but simply be a combination of the two.

The naming of the methods in the System class needs to be improved, especially with the equation interpretation functions. It would also be nice to separate the direct and standard assembly routines into separate methods to make thinks easier to read. Also, the application of functions, vectors, and gradients should be placed into a single method.