Following the recent postings regarding tochnog and contact issues, I remembered having similar issues before and reviewed the options and problems.
I realized that contact between a mesh and a point works, as it is shown in example 9. Contact between a mesh and a circle worked also (Please see attached problem)
contact with a line works perfectly well (see example 14) but contact between an ellipse and a mesh had issues and did not work.
After long hours looking at the problem, I have edited and corrected the errors in the code in file "geometry.cc" and it works very well, as it was supposed to.
I am attaching a copy of the corrected file, so that those building their own versions can update their code.
I will be updating the Mac OSX version binary very soon.
I would ask those of you who have access to the 64 bit versions of linux or windows and would like to contribute with binaries, please send them to me.
I will post them so that other users can download them. I would ask that you test your binaries before sending them to me.
If you can, for those using Windows, I know that there is a LAPACK version already built for windows. It would be nice if it can be built with LAPACK for
Windows, so that all versions would have the same capabilities.
If you do not have the LAPACK library option in your build please edit the corresponding line in the attached example files. The three small attached examples
solve the same problem for contact with a point 0.3 units in diameters, a circle of 0.3 units radius, and an ellipse with a=0.5 and b=0.3 all having the same center.
For those trying the contact between meshes, remember to include the "node_boundary" line in the problem.
Let me have some feedback please!