Well, this is a basic test.
But, it does not work.
Here is the in file:
#particle distributions and insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 8000 radius constant 0.005
fix pts2 all particletemplate/sphere 15485867 atom_type 1 density constant 8000 radius constant 0.0025
fix pts3 all particletemplate/sphere 15485867 atom_type 1 density constant 8000 radius constant 0.001
fix pdd1 all particledistribution/discrete/numberbased 32452843 3 pts1 0.3 pts2 0.5 pts3 0.2
fix ins all insert/pack seed 32452867 distributiontemplate pdd1 vel constant 0. 0. -0.3 &
insert_every once overlapcheck yes all_in yes particles_in_region 500 region bc
#apply nve integration to all particles
fix integr all nve/sphere
#output settings
compute rke all erotate/sphere
thermo_style custom step atoms ke vol
thermo 1000
~~~~~~~~~~~~~~~~~~~~~~~~
Error: Fix %s, ID %s: Random number generation: It is required that all the random seeds of this fix insert/*, \n"
" the random seed of particle distribution fix (id %s) template and all random seeds of the \n"
" fix particletemplate/* commands used by particle distribution fix (id %s) are different\n"
" Hint: possible valid (different) seeds would be the following numbers:\n"
" 15485863, 15485867, 32452843, 32452867, 49979687, 49979693, 67867967, 67867979, 86028121, 86028157",
~~~~~~~~~~~~~~~~~~~~~~~~
Honestly, I have no idea of picking a good seed number.
I also check the fix_insert.cpp.
I got this:
if(std::unique(seeds.begin(),seeds.end()) !=seeds.end() )
{
char errstr[1024];
sprintf(errstr,"Fix %s, ID %s: Random number generation: It is required that all the random seeds of this fix insert/*, \n"
" the random seed of particle distribution fix (id %s) template and all random seeds of the \n"
" fix particletemplate/* commands used by particle distribution fix (id %s) are different\n"
" Hint: possible valid (different) seeds would be the following numbers:\n"
" 15485863, 15485867, 32452843, 32452867, 49979687, 49979693, 67867967, 67867979, 86028121, 86028157",
style,id,fix_distribution->id,fix_distribution->id);
if(input->seed_check_throw_error())
error->one(FLERR,errstr);
else
error->warning(FLERR,errstr);
}
still very confused,,
mschramm | Tue, 03/21/2017 - 03:07
Can't use the same seed
You used the same seed for both pts2 and pts3. Try using 32452843 for pts3
adam05 | Tue, 03/21/2017 - 11:09
thx
thx
I change that and it works.
But I still feel confused about the seed number.
LIGGGHTS only gives out 10 option for the seed numbers.
So we can only add 8 different particles inside? (one for particle distribution, another one for insert/pack)
jsams | Tue, 03/21/2017 - 03:24
Unique seed needed
Your fix pts2 and pts3 have the same seed number. As stated by the error message they need to be different. The message also gives possible valid seed numbers, e.g. 15485863, 15485867, 32452843, 32452867, 49979687, 49979693, 67867967, 67867979, 86028121, 86028157
Choose any of those making sure you don't have any duplicates and try to run your simulations again. I hope this helps. Cheers, J
adam05 | Tue, 03/21/2017 - 11:10
thx
thx
It works.
But I still feel confused about the seed number.
LIGGGHTS only gives out 10 option for the seed numbers?
So we can only add 8 different particles inside? (one for particle distribution, another one for insert/pack)
mschramm | Tue, 03/21/2017 - 15:44
Just Examples
The numbers that are listed are examples of numbers that you can use.
The requirements for a number to be acceptable are;
1) The number must be greater than 10,000.
2) The number must be prime.
As an example, you can use any of the following
10007, 10009, 10037, 10039, 10061, 10067, 10069, 10079, 10091, 10093
As for why the change from using any number to having to use a prime number greater than 10,000, I have no clue. That would have to be answered by someone from CFDEM.
adam05 | Tue, 03/21/2017 - 16:04
Thx
Thx
That is what I need, I will try.