Changeset d98327

Timestamp:

Jul 7, 2009, 8:17:53 AM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

a811a0

Parents:

54de04

Message:

BUGFIX: LinkedCell list had to created with 2.*RADIUS of sphere

Some points were missed during the checking for the third point in tesselation due to the linked cell edge length being too small. We have to check TWICE the radius (i.e. the diameter) of the rolling sphere. NOTE: This was even suggested by a comment in that very line ... argh

File:

• molecuilder/src/builder.cpp (modified) (1 diff)

molecuilder/src/builder.cpp

@@ -1491 +1491 @@
                                                                 cout << Verbose(0) << "Evaluating non-convex envelope.";
                                                                 cout << Verbose(1) << "Using rolling ball of radius " << atof(argv[argptr]) << " and storing tecplot data in " << argv[argptr+1] << "." << endl;
-                                                                LinkedCell LCList(mol, atof(argv[argptr]));             // \NOTE not twice the radius??
+                                                                LinkedCell LCList(mol, atof(argv[argptr])*2.);
                                                                 Find_non_convex_border((ofstream *)&cout, mol, &T, &LCList, argv[argptr+1], atof(argv[argptr]));
                                                                 //FindDistributionOfEllipsoids((ofstream *)&cout, &T, &LCList, N, number, filename.c_str());