Modification of source code not working

Submitted by aishk on Wed, 08/28/2019 - 20:28

Hi,
I need to extend the source code on LIGGGHTS so I decided to start with looking at adding a few fprintf statements for debugging but any fprintf statements I add to the source code doesn't seem to be showing up on the terminal screen or on the logfile.
This is how I am remaking LIGGGHTS after any modifications:
make clean-auto
make stubs
make auto

Is there something I am doing wrong or not doing?
Any help is appreciated!
Thanks

mschramm | Thu, 09/12/2019 - 15:18

Hello,
what change did you make and are you sure that you are calling that change?

Also, if you simply make a small change, you do not need to use the clean command.
This could also check if you are making a correct change.
Like in the following
make auto
After building, make a change to particleToInsert.cpp
make auto
I should now see itself, and files associated with particleToInsert.cpp, being recompiled.

Kashminder | Sat, 10/09/2021 - 21:18

I have an analytical solution of my velocity profile of the fluid , so I need not to use openfoam and is just using the analytical solution for the fluid part and the LIGGGHTS for the particles part. My question is that I have defined the whole analytical solution in the input script (having a large number of variables) but it is not giving the accurate results. So, shall I define a different fix for the velocity profile? or it is advisable to use it in this way? below is my analytical solution.

variable ueq_r atom (3/8)*(1/(1-v_teq))*(1/v_req)*((1-v_req^2)-(1-v_req^2)^(-v_lmbda))*((v_zeq^2)/(v_heq^2)-2*v_zeq/v_heq)+((v_req*(v_h0^2)*v_heq/v_R^2)*(v_Jeq*v_lmbda*(1-v_req^2)^(-v_lmbda-1)+1)*(v_zeq/v_heq-(3/2)*((v_zeq^2)/(v_heq^2))))

variable ur atom v_ueq_r*v_R/v_tf

variable ux atom (v_ur*v_pic1/v_pic)*1000

variable ueq_z atom (3/4)*(1/(1-v_teq))*(1+v_lmbda*(1-v_req^2)^(-v_lmbda-1))*((v_zeq^3)/(3*v_heq^2)-(v_zeq^2)/(v_heq))+(3/2)*(1/(1-v_teq))*((1-v_req^2)-(1-v_req^2)^(-v_lmbda))*((v_zeq^2)/(2*v_heq^2)-(v_zeq^3)/(3*v_heq^3))*v_heq0-((v_h0^2/v_R^2)*(v_Jeq*v_lmbda*(1-v_req^2)^(-v_lmbda-1)+1)*(v_zeq^2-(v_zeq^3)/v_heq)+v_req*((v_h0^2)/v_R^2)*v_Jeq*v_lmbda*(v_lmbda+1)*((1-v_req^2)^(-v_lmbda-2))*(v_zeq^2-(v_zeq^3)/v_heq)-v_req*((v_h0^2)/v_R^2)*(v_Jeq*v_lmbda*((1-v_req^2)^(-v_lmbda-1))+1)*((v_zeq^3)/(v_heq^2))*v_heq0)

variable uz atom (v_ueq_z*v_h0/v_tf)*1000

variable uy atom (v_ur*v_nic1/v_pic)*1000

I have already defined all these variables above in the script.