search

huilinye / OpenFSI

A highly efficient and portable fluid-structure simulation package based on immersed-boundary method
GNU General Public License v3.0
32 stars 20 forks source link

Not compatible with the latest stable Lammps #1

Open Yeyke opened 4 years ago

Yeyke commented 4 years ago

Hi,

There is an error when compiling the LAMMPS as a library with the latest stable Lammps version (3 Mar 2020). The error log is shown below:

g++: warning: /usr/include: linker input file unused because linking not done mpic++ -g -O3 -DLAMMPS_GZIP -DLAMMPS_JPEG -DLAMMPS_PNG -DLAMMPS_FFMPEG /lib/x86_64-linux-gnu/openmpi/include/ -DFFT_FFTW3 /usr/include -c ../angle_cosine.cpp ../angle_cosine.cpp: In member function ‘virtual void LAMMPS_NS::AngleCosine::compute(int, int)’: ../angle_cosine.cpp:56:3: error: ‘ev_init’ was not declared in this scope; did you mean ‘init’? 56 | ev_init(eflag,vflag); | ^~~ | init make[1]: [Makefile:120: angle_cosine.o] Error 1 make[1]: Leaving directory '/home/test/lammps/lammps-3Mar20/src/Obj_ubuntu' make: [Makefile:217: ubuntu] Error 2

Can you update the openFSI with the new Lammps? Thanks very much,

Weicai Ye

huilinye commented 4 years ago

Hi,

As some functions used in this package are removed by new Lammps, you need to use an old version like Jan2017. We were also considering adjusting to the new Lammps, but many parts need refreshing. It's a hard work to do that, so please use the old version instead.

Best,

ChiaCheng-Tsai commented 4 years ago

I used the latest stable Lammps version (3 Mar 2020) and it worked well for the 2d example after minor modification. For examples, in the files of in.* "neigh_modify exclude molecule move" should be changed to "neigh_modify exclude molecule/inter move".

ChiaCheng-Tsai commented 4 years ago

However, I cannot figure out how "fix ha all vtk 1000 Hex meshfile octa_ele.data output solid/cy" can work as the system complain "Unrecognized fix style 'vtk'". And I cannot find it out in Jan2017 version, either

ChiaCheng-Tsai commented 4 years ago

For the 2D case, it results in warming as "Warming: Bond/angle/dihedral extent > half of periodic box length, domain.cpp:906" And finally on iT = 109237, the problem stoped by the error as "Bond atom missing in image check, domain.cpp:765"