Commit 2eadc337 authored by Berenger Bramas's avatar Berenger Bramas

use old prolate non unif generator

parent fa0d77f3
...@@ -124,7 +124,7 @@ int main(int argc, char* argv[]){ ...@@ -124,7 +124,7 @@ int main(int argc, char* argv[]){
FReal * tmpParticles = new FReal[4*NbParticlesPerNode]; FReal * tmpParticles = new FReal[4*NbParticlesPerNode];
if(FParameters::existParameter(argc, argv, "-ellipsoid")) { if(FParameters::existParameter(argc, argv, "-ellipsoid")) {
std::cout << "ellipsoid\n"; std::cout << "ellipsoid\n";
nonunifRandomPointsOnElipsoid(NbParticlesPerNode, boxWidth/2, boxWidth/4, boxWidth/8, tmpParticles); nonunifRandomPointsOnElipsoid(NbParticlesPerNode, 0.5, 0.1, tmpParticles);
} }
else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) { else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) {
std::cout << "ellipsoidv2\n"; std::cout << "ellipsoidv2\n";
......
...@@ -127,7 +127,7 @@ int main(int argc, char* argv[]){ ...@@ -127,7 +127,7 @@ int main(int argc, char* argv[]){
memset(tmpParticles,0,(unsigned int) (sizeof(FReal)* loader.getNumberOfParticles() * 4)); memset(tmpParticles,0,(unsigned int) (sizeof(FReal)* loader.getNumberOfParticles() * 4));
if(FParameters::existParameter(argc, argv, "-ellipsoid")) { if(FParameters::existParameter(argc, argv, "-ellipsoid")) {
std::cout << "ellipsoid\n"; std::cout << "ellipsoid\n";
nonunifRandomPointsOnElipsoid(loader.getNumberOfParticles(), boxWidth/2, boxWidth/4, boxWidth/8, tmpParticles); nonunifRandomPointsOnElipsoid(loader.getNumberOfParticles(), 0.5, 0.1, tmpParticles);
} }
else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) { else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) {
std::cout << "ellipsoidv2\n"; std::cout << "ellipsoidv2\n";
......
...@@ -103,7 +103,7 @@ int main(int argc, char* argv[]){ ...@@ -103,7 +103,7 @@ int main(int argc, char* argv[]){
memset(tmpParticles,0,(unsigned int) (sizeof(FReal)* loader.getNumberOfParticles() * 4)); memset(tmpParticles,0,(unsigned int) (sizeof(FReal)* loader.getNumberOfParticles() * 4));
if(FParameters::existParameter(argc, argv, "-ellipsoid")) { if(FParameters::existParameter(argc, argv, "-ellipsoid")) {
std::cout << "ellipsoid\n"; std::cout << "ellipsoid\n";
nonunifRandomPointsOnElipsoid(loader.getNumberOfParticles(), boxWidth/2, boxWidth/4, boxWidth/8, tmpParticles); nonunifRandomPointsOnElipsoid(loader.getNumberOfParticles(), 0.5, 0.1, tmpParticles);
} }
else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) { else if(FParameters::existParameter(argc, argv, LocalOptionEllipsoidv2.options)) {
std::cout << "ellipsoidv2\n"; std::cout << "ellipsoidv2\n";
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment