Revision: 55183
http://brlcad.svn.sourceforge.net/brlcad/?rev=55183&view=rev
Author: starseeker
Date: 2013-04-18 20:51:30 +0000 (Thu, 18 Apr 2013)
Log Message:
-----------
Add test case for get_closest_point, using an isolated surface and two 3D
inputs that are unexpectedly returning the same 2D pullback result.
Modified Paths:
--------------
brlcad/trunk/src/librt/CMakeLists.txt
Added Paths:
-----------
brlcad/trunk/src/librt/tests/
brlcad/trunk/src/librt/tests/CMakeLists.txt
brlcad/trunk/src/librt/tests/nurbs_surfaces.g
brlcad/trunk/src/librt/tests/nurbs_tests.cpp
Modified: brlcad/trunk/src/librt/CMakeLists.txt
===================================================================
--- brlcad/trunk/src/librt/CMakeLists.txt 2013-04-18 19:41:41 UTC (rev
55182)
+++ brlcad/trunk/src/librt/CMakeLists.txt 2013-04-18 20:51:30 UTC (rev
55183)
@@ -313,10 +313,13 @@
#BRLCAD_ADDEXEC(test_bot2nurbs test_bot2nurbs.cpp "librt;libwdb;libbrep;libbu"
NOSTRICT NO_INSTALL)
#BRLCAD_ADDEXEC(test_nurbsfit test_nurbsfit.cpp "librt;libwdb;libbrep;libbu"
NOSTRICT NO_INSTALL)
#BRLCAD_ADDEXEC(test_root3-subd test_root3-subd.cpp
"librt;libwdb;libbrep;libbu" NOSTRICT NO_INSTALL)
+#BRLCAD_ADDEXEC(test_surfacetree test_surfacetree.cpp
"librt;libwdb;libbrep;libbu" NOSTRICT NO_INSTALL)
BRLCAD_ADDEXEC(test_dvec test_dvec.cpp "librt" NOSTRICT NO_INSTALL LOCAL)
ADD_TEST(test_dvec test_dvec)
+add_subdirectory(tests)
+
# Local Variables:
# tab-width: 8
# mode: cmake
Added: brlcad/trunk/src/librt/tests/CMakeLists.txt
===================================================================
--- brlcad/trunk/src/librt/tests/CMakeLists.txt (rev 0)
+++ brlcad/trunk/src/librt/tests/CMakeLists.txt 2013-04-18 20:51:30 UTC (rev
55183)
@@ -0,0 +1,12 @@
+BRLCAD_ADDEXEC(nurbs_tester nurbs_tests.cpp "librt;libbrep;libbu" NO_INSTALL
LOCAL)
+
+add_test(get_closest_point-distinct_points nurbs_tester
${CMAKE_CURRENT_SOURCE_DIR}/nurbs_surfaces.g 1)
+
+CMAKEFILES(nurbs_surfaces.g)
+
+# Local Variables:
+# tab-width: 8
+# mode: cmake
+# indent-tabs-mode: t
+# End:
+# ex: shiftwidth=2 tabstop=8 textwidth=0 wrapmargin=0
Property changes on: brlcad/trunk/src/librt/tests/CMakeLists.txt
___________________________________________________________________
Added: svn:mime-type
+ text/plain
Added: svn:eol-style
+ native
Added: brlcad/trunk/src/librt/tests/nurbs_surfaces.g
===================================================================
(Binary files differ)
Property changes on: brlcad/trunk/src/librt/tests/nurbs_surfaces.g
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: brlcad/trunk/src/librt/tests/nurbs_tests.cpp
===================================================================
--- brlcad/trunk/src/librt/tests/nurbs_tests.cpp
(rev 0)
+++ brlcad/trunk/src/librt/tests/nurbs_tests.cpp 2013-04-18 20:51:30 UTC
(rev 55183)
@@ -0,0 +1,139 @@
+#include "common.h"
+
+#include <cstdlib>
+#include <cstdio>
+#include <cstring>
+#include <limits>
+#include <iomanip>
+
+#include "brep.h"
+#include "raytrace.h"
+
+
+int
+get_surface(const char *name, struct db_i *dbip, struct rt_db_internal
*intern, struct rt_brep_internal **brep_ip) {
+
+ struct directory *dp;
+ (*brep_ip) = NULL;
+
+
+ dp = db_lookup(dbip, name, LOOKUP_QUIET);
+ if (dp == RT_DIR_NULL) {
+ bu_log("ERROR: Unable to look up object %s\n", name);
+ return -1;
+ }
+
+ if (rt_db_get_internal(intern, dp, dbip, NULL, &rt_uniresource) < 0) {
+ bu_log("ERROR: Unable to get internal representation of %s\n", name);
+ return -1;
+ }
+
+ if (intern->idb_minor_type != DB5_MINORTYPE_BRLCAD_BREP) {
+ bu_log("ERROR: object %s does not appear to be of type BRep\n", name);
+ return -1;
+ } else {
+ (*brep_ip) = (struct rt_brep_internal *)(intern->idb_ptr);
+ }
+
+ RT_BREP_CK_MAGIC(*brep_ip);
+ return 0;
+}
+
+int
+nurbs_test(long int test_number, struct db_i *dbip)
+{
+ /* Have various specific surface structures - there are potentially
+ * multiple surfaces in the .g, but multiple tests may re-use a single
+ * surface.*/
+
+ /* openNURBS pointers */
+ ON_Brep *brep = NULL;
+
+ switch (test_number)
+ {
+ case 1: /* 3d -> 2d pull-back function get_closest_point */
+ {
+ struct rt_db_internal case_1_intern;
+ struct rt_brep_internal *case_1_brep_ip = NULL;
+ RT_DB_INTERNAL_INIT(&case_1_intern);
+ if (get_surface("case_1_surface.s", dbip, &case_1_intern,
&case_1_brep_ip)) return -1;
+ brep = case_1_brep_ip->brep;
+ ON_BrepFace& c1_face = brep->m_F[0];
+ brlcad::SurfaceTree* c1_st = new brlcad::SurfaceTree(&c1_face,
false);
+
+ ON_2dPoint pt_2d_1, pt_2d_2;
+ ON_3dPoint pt_3d_1(11204.05366897583,-16726.489562988281,
3358.7263298034668);
+ ON_3dPoint pt_3d_2(11204.007682800293, -16726.479568481445,
3358.8327312469482);
+ (void)get_closest_point(pt_2d_1, &c1_face, pt_3d_1, c1_st);
+ (void)get_closest_point(pt_2d_2, &c1_face, pt_3d_2, c1_st);
+ rt_db_free_internal(&case_1_intern);
+ delete c1_st;
+ if (pt_2d_1 == pt_2d_2) {
+ std::cout << "NURBS test case 1 failure
(get_closest_point): Unexpectedly identical 2D pullbacks from different 3D
points\n";
+ std::cout <<
std::setprecision(std::numeric_limits<double>::digits10) << "Inputs: 3d Point
1: " << pt_3d_1.x << "," << pt_3d_1.y << "," << pt_3d_1.z << "\n";
+ std::cout <<
std::setprecision(std::numeric_limits<double>::digits10) << " 3d Point
2: " << pt_3d_2.x << "," << pt_3d_2.y << "," << pt_3d_2.z << "\n";
+ std::cout <<
std::setprecision(std::numeric_limits<double>::digits10) << "Output: 2d Point
1: " << pt_2d_1.x << "," << pt_2d_1.y << "\n";
+ std::cout <<
std::setprecision(std::numeric_limits<double>::digits10) << " 2d Point
2: " << pt_2d_2.x << "," << pt_2d_2.y << "\n";
+ return -1;
+ }
+ return 0;
+ }
+ default:
+ return -1;
+ }
+
+}
+
+
+int
+main(int argc, char **argv)
+{
+ long int test_number = 0;
+ int retval = 0;
+ long int all_tests = 1;
+ struct db_i *dbip;
+ char *endptr;
+
+ /* Set up the .g file */
+ if (argc > 3 || argc < 2) {
+ bu_log("Usage: %s file.g [test_number]", argv[0]);
+ return -1;
+ }
+
+ dbip = db_open(argv[1], "r+w");
+ if (dbip == DBI_NULL) {
+ bu_log("ERROR: Unable to read from %s\n", argv[1]);
+ return -1;
+ }
+
+ if (db_dirbuild(dbip) < 0){
+ bu_log("ERROR: Unable to read from %s\n", argv[1]);
+ return -1;
+ }
+
+ if (argc == 3) {
+ test_number = strtol(argv[2], &endptr, 10);
+ if (*endptr) return -1;
+ }
+
+ if (!test_number) {
+ int current_test = 1;
+ while (current_test <= all_tests && retval != -1) {
+ retval = nurbs_test(current_test, dbip);
+ current_test++;
+ }
+ } else {
+ retval = nurbs_test(test_number, dbip);
+ }
+ return retval;
+}
+
+/*
+ * Local Variables:
+ * mode: C
+ * tab-width: 8
+ * indent-tabs-mode: t
+ * c-file-style: "stroustrup"
+ * End:
+ * ex: shiftwidth=4 tabstop=8
+ */
Property changes on: brlcad/trunk/src/librt/tests/nurbs_tests.cpp
___________________________________________________________________
Added: svn:mime-type
+ text/plain
Added: svn:eol-style
+ native
This was sent by the SourceForge.net collaborative development platform, the
world's largest Open Source development site.
------------------------------------------------------------------------------
Precog is a next-generation analytics platform capable of advanced
analytics on semi-structured data. The platform includes APIs for building
apps and a phenomenal toolset for data science. Developers can use
our toolset for easy data analysis & visualization. Get a free account!
http://www2.precog.com/precogplatform/slashdotnewsletter
_______________________________________________
BRL-CAD Source Commits mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/brlcad-commits
