Revision: 56500
http://sourceforge.net/p/brlcad/code/56500
Author: starseeker
Date: 2013-08-03 13:03:08 +0000 (Sat, 03 Aug 2013)
Log Message:
-----------
Add patch #221 from Nyah Check (Ch3ck) implementing a libbn test for
bn_poly_multiply.
Modified Paths:
--------------
brlcad/trunk/src/libbn/tests/CMakeLists.txt
Added Paths:
-----------
brlcad/trunk/src/libbn/tests/bn_poly_multiply.c
Modified: brlcad/trunk/src/libbn/tests/CMakeLists.txt
===================================================================
--- brlcad/trunk/src/libbn/tests/CMakeLists.txt 2013-08-03 10:46:14 UTC (rev
56499)
+++ brlcad/trunk/src/libbn/tests/CMakeLists.txt 2013-08-03 13:03:08 UTC (rev
56500)
@@ -1,6 +1,7 @@
BRLCAD_ADDEXEC(tester_bn_tri_tri_isect_coplanar bn_tri_tri_isect_coplanar.c
"libbu;libbn" NO_INSTALL)
BRLCAD_ADDEXEC(tester_bn_tri_tri_isect bn_tri_tri_isect.c "libbu;libbn"
NO_INSTALL)
BRLCAD_ADDEXEC(tester_bn_list bn_list.c "libbu;libbn" NO_INSTALL)
+BRLCAD_ADDEXEC(tester_bn_poly_multiply bn_poly_multiply.c "libbu;libbn"
NO_INSTALL)
# For tester_bn_tri_tri_isect_coplanar, the input format is as follows:
#
@@ -58,6 +59,14 @@
add_test(bn_list_3d_3 tester_bn_list 3 3,7,8 1,-5,8)
add_test(bn_list_3d_4 tester_bn_list 3 3,7,8 1,-5.321,3
67,0,67 0,100,23)
+
+#
+# *********** poly.c tests *************
+#
+#Format is:
+# tester_bn_poly < void > void
+add_test(bn_poly_multiply tester_bn_poly_multiply)
+
# Local Variables:
# tab-width: 8
# mode: cmake
Added: brlcad/trunk/src/libbn/tests/bn_poly_multiply.c
===================================================================
--- brlcad/trunk/src/libbn/tests/bn_poly_multiply.c
(rev 0)
+++ brlcad/trunk/src/libbn/tests/bn_poly_multiply.c 2013-08-03 13:03:08 UTC
(rev 56500)
@@ -0,0 +1,165 @@
+/* T E S T _ B N _ P O L Y _ M U L T I P L Y . C
+ * BRL-CAD
+ *
+ * Copyright (c) 2004-2013 United States Government as represented by
+ * the U.S. Army Research Laboratory.
+ *
+ * This library is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public License
+ * version 2.1 as published by the Free Software Foundation.
+ *
+ * This library is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with this file; see the file named COPYING for more
+ * information.
+ */
+/**
+ *@Brief
+ *The known values used for these tests are generated from Octave
+ *GNU Octave, version 3.4.3
+ *Copyright (C) 2011 John W. Eaton and others.
+ *Octave was configured for "i386-redhat-linux-gnu".
+ */
+
+
+#include "common.h"
+
+#include <stdio.h>
+#include <math.h>
+#include <string.h>
+#include <signal.h>
+
+#include "bu.h"
+#include "vmath.h"
+#include "bn.h"
+#include "magic.h"
+
+
+struct bn_poly bn_Zero_poly = { BN_POLY_MAGIC, 0, {0.0} };
+
+/*holds three polynomials to be used in test.*/
+bn_poly_t input[3], output[3];
+
+
+/**
+ *Initialises polnomial storing a negative, positive and zero coefficients.
+ */
+void
+poly_init(void)
+{
+
+ /*stores o coefficients to polynomial for input and output.*/
+ output[0] = bn_Zero_poly;
+ input[0] = bn_Zero_poly;
+ input[0].dgr = 2;
+ input[0].cf[0] = input[0].cf[1] = input[0].cf[2] = input[0].cf[3] = 0.0;
+
+ output[0].dgr = 4;
+ output[0].cf[0] = output[0].cf[1] = output[0].cf[2] = output[0].cf[3] =
output[0].cf[4] = 0.0;
+
+ /*stores negative coefficients to polynomial.*/
+ output[1] = bn_Zero_poly;
+ input[1] = bn_Zero_poly;
+ input[1].dgr = 2;
+ output[1].dgr = 4;
+
+ input[1].cf[0] = -4;
+ input[1].cf[1] = -3;
+ input[1].cf[2] = -2;
+
+ output[1].cf[0] = 16;
+ output[1].cf[1] = 24;
+ output[1].cf[2] = 25;
+ output[1].cf[3] = 12;
+ output[1].cf[4] = 4;
+
+ /*stores positive coefficients to to polynomial input.*/
+ output[2] = bn_Zero_poly;
+ input[2] = bn_Zero_poly;
+ input[2].dgr = 2;
+ output[2].dgr = 4;
+
+ input[2].cf[0] = 7854;
+ input[2].cf[1] = 2136;
+ input[2].cf[2] = 1450;
+
+ output[2].cf[0] = 61685316;
+ output[2].cf[1] = 33552288;
+ output[2].cf[2] = 27339096;
+ output[2].cf[3] = 6194400;
+ output[2].cf[4] = 2102500;
+
+}
+
+
+/* compares the values of the array and returns 0. */
+int
+check_results(fastf_t a[], fastf_t b[], int n)
+{
+ int i;
+
+ for (i = 0; i < n; i++) {
+ if (!EQUAL(a[i],b[i]))
+ return -1;
+ }
+
+ return 0;
+}
+
+
+/*tests the polynomials to make sure bn_poly_mul() works properly.*/
+int test_bn_poly(void)
+{
+ int val, val1, val2;
+ bn_poly_t a, b, c;
+ a = bn_Zero_poly, b = bn_Zero_poly, c = bn_Zero_poly;
+
+ bn_poly_mul(&a, &input[0], &input[0]);
+ bn_poly_mul(&b, &input[1], &input[1]);
+ bn_poly_mul(&c, &input[2], &input[2]);
+
+ val = check_results(a.cf,output[0].cf, output[0].dgr + 1);
+ val1 = check_results( b.cf, output[1].cf, output[1].dgr + 1);
+ val2 = check_results(c.cf, output[2].cf, output[2].dgr + 1);
+
+ if (val == 0 && val1 == 0 && val2 == 0)
+ return val;
+
+ return -1;
+}
+
+
+int main(void)
+{
+ int ret;
+
+ poly_init();
+
+ ret = test_bn_poly();
+
+ if( ret == 0){
+ bu_log("\nFunction computes correctly\n");
+
+ return ret;
+ }
+
+ else {
+ bu_log("\nInvalid output.\n");
+ }
+ return -1;
+
+}
+
+/*
+ * 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/libbn/tests/bn_poly_multiply.c
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+text/plain
\ No newline at end of property
Added: svn:eol-style
## -0,0 +1 ##
+native
\ No newline at end of property
This was sent by the SourceForge.net collaborative development platform, the
world's largest Open Source development site.
------------------------------------------------------------------------------
Get your SQL database under version control now!
Version control is standard for application code, but databases havent
caught up. So what steps can you take to put your SQL databases under
version control? Why should you start doing it? Read more to find out.
http://pubads.g.doubleclick.net/gampad/clk?id=49501711&iu=/4140/ostg.clktrk
_______________________________________________
BRL-CAD Source Commits mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/brlcad-commits
