inlcude niput files not working with lammps as lib

Hi everyone,

I'm calling lammps as a shared lib from OpenFOAM. I would like to use
the "include" command as described here"
http://lammps.sandia.gov/doc/include.html " to include some variables
and fixes.

When I run the input script with lammps as standalone it includes the
files but when I call the script with lammps as a shared lib it doesn't
work. I've attached the input script and include files.

Any thoughts or tips? Thanks in advance.

in.lammps (875 Bytes)

variables (1.41 KB)

fixes (522 Bytes)

Hi everyone,

I'm calling lammps as a shared lib from OpenFOAM. I would like to use
the "include" command as described here"
http://lammps.sandia.gov/doc/include.html " to include some variables
and fixes.

When I run the input script with lammps as standalone it includes the
files but when I call the script with lammps as a shared lib it doesn't
work. I've attached the input script and include files.

Any thoughts or tips? Thanks in advance.

​i cannot reproduce it with the current version of LAMMPS (10 March 2017).
if i comment out the 'read_restart' command (for lack of having your
restart file), i can see "include" working via the library interface. e.g.
when using this little python script:

​[[email protected] Downloads]$ cat test.py
from lammps import *
lmp=lammps()
lmp.command('echo both')
lmp.file('in.lammps')

​which produces the following screen output indicating that the "include"
directive is properly processed.

axel.

[[email protected] Downloads]$ python test.py
LAMMPS (10 Mar 2017)
  using 1 OpenMP thread(s) per MPI task
# input script lammpsFoam simulations
package omp 0
using multi-threaded neighbor list subroutines
atom_style sphere
atom_modify map array
boundary f f f
newton off
processors * * *
units cgs
comm_modify vel yes

# Read restart and import variables
variable path getenv "PWD"
print "path is \{path\}" path is /home/akohlmey/Downloads \# read\_restart restart\.grow include {path}/variables
include /home/akohlmey/Downloads/variables
# Variables used in multiple scripts

# domain variables
variable H0 equal 11.5 # [cm]
variable EF equal 3.0 # [-]
variable hb equal v_H0*v_EF # [cm]
variable height equal 1.1*v_H0*v_EF # [cm], add 10% extra height to domain
variable width equal 1.0 # [cm]
variable depth equal 1.0 # [cm]
variable eps equal 1.-$((1.-0.4)/v_EF) # average vol fraction, needed for
particle insertion eqn
variable eps equal 1.-0.19999999999999998335
variable rhop equal 2.8 # [g/cm3]

# particle insertion variables
variable dslice equal 0.2 # thicknes of a single slice
variable slices equal floor(\{hb\}/{dslice}) # nr of insertion slices
variable slices equal floor(34.5/\{dslice\}\) variable slices equal floor\(34\.5/0\.2\) variable vslice equal (v_width)*\(v\_width\)\*v\_dslice \# volume of a single slice variable vslice equal 1\*(v_width)*v_dslice
variable vslice equal 1*1*v_dslice
variable dred equal 0.5 # reduction factor for insertion size
variable CG equal 2. # particle coarse graining factor

# Growth variables
variable dred equal 0.1 # reduction factor for insertion size
variable growts equal 25000
variable growevery equal 1
variable relaxts1 equal 10000
variable relaxts2 equal 100000

# granular variables (based on eqns from Landry 2003 )
variable dpav equal 0.025
variable mpav equal $(v_rhop*PI/6.*v_dpav*v_dpav*v_dpav)
variable mpav equal 2.2907446432425579132e-05
variable g equal 981.

variable kn equal \(2\.e5\*v\_mpav\*v\_g/v\_dpav\) variable kn equal 179777\.63960167608457 variable kt equal (2./7.*v_kn)
variable kt equal 51365.039886193138955
variable gamman equal 50*\(sqrt\(v\_g/v\_dpav\)\) variable gamman equal 50\*198\.09088823063012796 variable gammat equal {gamman}
variable gammat equal 9904.54441153151
variable mu equal 0.15

print "included variables"
included variables
include \{path\}/fixes include /home/akohlmey/Downloads/fixes \# pair\_style none pair\_style gran/hertz/history/omp {kn} \{kt\} {gamman} \{gammat\} {mu} 1
pair_style gran/hertz/history/omp 179777.639601676 \{kt\} {gamman}
\{gammat\} {mu} 1
pair_style gran/hertz/history/omp 179777.639601676 51365.0398861931
\{gamman\} {gammat} \{mu\} 1 pair\_style gran/hertz/history/omp 179777\.639601676 51365\.0398861931 9904\.54441153151 {gammat} \{mu\} 1 pair\_style gran/hertz/history/omp 179777\.639601676 51365\.0398861931 9904\.54441153151 9904\.54441153151 {mu} 1
pair_style gran/hertz/history/omp 179777.639601676 51365.0398861931
9904.54441153151 9904.54441153151 0.15 1
pair_coeff * *
ERROR: Pair_coeff command before simulation box is defined
(../input.cpp:1762)