new_parralelism

test_FMM.config
View Options

			Section MultiScale AtomsUnits
			DIMENSION 2
			UNITCODE AtomsUnits

			#lattice constant for aluminum at 0K first neighbor interaction
			LET sqrt3 = 1.73205080756887729353
			LET sqrt2 = 1.4142135623730950488
			LET r0 = 4.032*sqrt2/2
			LET r1 = r0*sqrt3

			LET LX = 160*r0
			LET LY = 90*r1

			GEOMETRY 1 CUBE BBOX -LX/2 LX/2 -LY/2 LY/2

			COM DISTRIBUTED
			MD_CODE LAMMPS md


			LET when = 0
			LET factor = 1

			FOR k in range 1 160 1
			LET wavelengthX = 160/(2*k)
			LET wavelengthY = 90/(2*k)
			STIMULATION harmonic${when}-${wavelengthX}-${wavelengthY} SIN INPUT md WAVELENGTH wavelengthXr0 wavelengthYr1 0 ONESHOT when DIRECTION 0 INTENSITY 1.0e-3 STAGE PRE_DUMP
			LET when = when+1
			endFOR

			LET order = 3
			LET nlevels = 4
			LET startlevel = 4

			LET fsize = 20*r0
			LET cut = 10*r0
			COMPUTE fmm FMM INPUT md GRIDSPACE 5r0 5r0 0 CUTOFF_WAVELENGTH cut FILTER_SIZE fsize ORDER order FILTER_NPOINT 1000 FUNCTION lowpass RCUT fsize PBC 1 1 0 NLEVELS nlevels START_LEVEL startlevel DENSITY r0*r1/2

			COMPUTE amplitudeMax REDUCE INPUT fmm OPERATOR MAX
			COMPUTE compute_wavelength PYTHON INPUT amplitudeMax FILENAME wl
			DUMPER dump_amplitude-${order}-${startlevel}-${nlevels} TEXT INPUT compute_wavelength PREFIX ./
			COMPUTE test EXTRACT INPUT md FIELD displacement COMPONENT 0

			DUMPER dumper_paraview PARAVIEW INPUT md MASS P0 DISP VEL FORCE FREQ 1 ADDFIELD fmm ADDFIELD test PREFIX ./
			endSection

			Section LAMMPS:md AtomsUnits
			LAMMPS_FILE in.lj
			DOMAIN_GEOMETRY 1
			endSection