From cca66b9ec4e494c1d919bff0f71a820d8afab1fa Mon Sep 17 00:00:00 2001 From: Daniel Baumann Date: Sun, 7 Apr 2024 20:24:48 +0200 Subject: Adding upstream version 1.2.2. Signed-off-by: Daniel Baumann --- src/3rdparty/2geom/tests/implicitization-test.cpp | 130 ++++++++++++++++++++++ 1 file changed, 130 insertions(+) create mode 100644 src/3rdparty/2geom/tests/implicitization-test.cpp (limited to 'src/3rdparty/2geom/tests/implicitization-test.cpp') diff --git a/src/3rdparty/2geom/tests/implicitization-test.cpp b/src/3rdparty/2geom/tests/implicitization-test.cpp new file mode 100644 index 0000000..bfc4c58 --- /dev/null +++ b/src/3rdparty/2geom/tests/implicitization-test.cpp @@ -0,0 +1,130 @@ +/* + * Test program for implicitization routines + * + * Authors: + * Marco Cecchetti + * + * Copyright 2008 authors + * + * This library is free software; you can redistribute it and/or + * modify it either under the terms of the GNU Lesser General Public + * License version 2.1 as published by the Free Software Foundation + * (the "LGPL") or, at your option, under the terms of the Mozilla + * Public License Version 1.1 (the "MPL"). If you do not alter this + * notice, a recipient may use your version of this file under either + * the MPL or the LGPL. + * + * You should have received a copy of the LGPL along with this library + * in the file COPYING-LGPL-2.1; if not, write to the Free Software + * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * You should have received a copy of the MPL along with this library + * in the file COPYING-MPL-1.1 + * + * The contents of this file are subject to the Mozilla Public License + * Version 1.1 (the "License"); you may not use this file except in + * compliance with the License. You may obtain a copy of the License at + * http://www.mozilla.org/MPL/ + * + * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY + * OF ANY KIND, either express or implied. See the LGPL or the MPL for + * the specific language governing rights and limitations. + */ + + + +#include <2geom/symbolic/implicit.h> + +#include "pick.h" + +#include + + + + +void print_basis(Geom::SL::basis_type const& b) +{ + for (size_t i= 0; i < 2; ++i) + { + for (size_t j= 0; j < 3; ++j) + { + std::cout << "b[" << i << "][" << j << "] = " << b[i][j] << "\n"; + } + } +} + + + + +int main( int argc, char * argv[] ) +{ + // degree of polinomial parametrization + // warning: not set N to a value greater than 20! + // (10 in case you don't utilize the micro-basis) + // determinant computation becomes very expensive + unsigned int N = 4; + // max modulus of polynomial coefficients + unsigned int M = 1000; + + if (argc > 1) + N = std::atoi(argv[1]); + if (argc > 2) + M = std::atoi(argv[2]); + + Geom::SL::MVPoly1 f, g; + Geom::SL::basis_type b; + Geom::SL::MVPoly3 p, q; + Geom::SL::Matrix B; + Geom::SL::MVPoly2 r; + + // generate two univariate polynomial with degree N + // and coeffcient in the range [-M, M] + f = pick_multipoly_max<1>(N, M); + g = pick_multipoly_max<1>(N, M); + + std::cout << "parametrization: \n"; + std::cout << "f = " << f << std::endl; + std::cout << "g = " << g << "\n\n"; + + // computes the micro-basis + microbasis(b, f, g); + // in case you want utilize directly the initial basis + // you should uncomment the next row and comment + // the microbasis function call + //make_initial_basis(b, f, g); + + std::cout << "generators in vector form : \n"; + print_basis(b); + std::cout << std::endl; + + // micro-basis generators + basis_to_poly(p, b[0]); + basis_to_poly(q, b[1]); + + std::cout << "generators as polynomial in R[t,x,y] : \n"; + std::cout << "p = " << p << std::endl; + std::cout << "q = " << q << "\n\n"; + + // make up the Bezout matrix and compute the determinant + B = make_bezout_matrix(p, q); + r = determinant_minor(B); + r.normalize(); + + std::cout << "Bezout matrix: (entries are bivariate polynomials) \n"; + std::cout << "B = " << B << "\n\n"; + std::cout << "determinant: \n"; + std::cout << "r(x, y) = " << r << "\n\n"; + + return EXIT_SUCCESS; +} + + +/* + Local Variables: + mode:c++ + c-file-style:"stroustrup" + c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +)) + indent-tabs-mode:nil + fill-column:99 + End: +*/ +// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99 : -- cgit v1.2.3