FlexPDE takes the user's boundary and feature descriptions and forms an initial triangular grid. If the specified solution accuracy cannot be achieved with this initial grid, FlexPDE will automatically subdivide those patches where the error exceeds the user specification and reiterate its solutions. By subdividing only those elements where the error is high, FlexPDE keeps the number of nodes to a minimum while still insuring the accuracy desired. Gridding of two small wells in a large field.
Zoomed plot shows fine grid around one of the wells.