FEAP User Forum

FEAP => Programming => Topic started by: Brian on September 17, 2015, 10:35:45 PM

Title: How to add a new solver
Post by: Brian on September 17, 2015, 10:35:45 PM

How can a new solver be added to FEAP? Please explain it in detail. Thank you.

Title: Re: How to add a new solver
Post by: Prof. R.L. Taylor on September 18, 2015, 06:35:44 AM

There are some routines at the feap web site for adding solvers.  They are a bit out of date but it should be possible to update for the current version.  We will update the site with more recent versions in the near future.

Title: Re: How to add a new solver
Post by: trustetj on September 20, 2015, 04:56:53 PM

I also posted an example to the forum for adding the Intel MKL Pardiso solver through the usolve.f subroutine. http://feap.berkeley.edu/forum/index.php?topic=761.msg2081 (http://feap.berkeley.edu/forum/index.php?topic=761.msg2081). This is a parallel direct solver that uses OpenMP functionality. There are also open source versions of the Pardiso package. The zip file contains the relevant subroutines compatible with FEAP 8.4 as well as a sample input file.

Note that since making the post, I have encountered memory errors when large meshes are run for several time steps, so the allocation of the memory for the stiffness matrix may not have been created in a completely stable fashion.

Title: Re: How to add a new solver
Post by: FEAP_Admin on September 20, 2015, 08:35:20 PM

Prof. Taylor posted an up-to-date user solver using Pardisio from MKL: http://www.ce.berkeley.edu/projects/feap/
Please try it to see if your memory issues go away.

Title: Re: How to add a new solver
Post by: Brian on October 21, 2015, 04:20:59 AM

Could you please explain how to add superlu84 to feap by visual studio? How about a new umacr and usolve? Thank you.

Title: Re: How to add a new solver
Post by: FEAP_Admin on October 21, 2015, 09:46:20 AM

add the files to your project and rebuild.