#ifndef around namespace GiNaC { }
[ginac.git] / check / lsolve_onedim.cpp
1 /** @file lsolve_onedim.cpp
2  *
3  * This test routine does some simple checks on solving a polynomial for a
4  * variable. */
5
6 /*
7  *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
8  *
9  *  This program is free software; you can redistribute it and/or modify
10  *  it under the terms of the GNU General Public License as published by
11  *  the Free Software Foundation; either version 2 of the License, or
12  *  (at your option) any later version.
13  *
14  *  This program is distributed in the hope that it will be useful,
15  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
16  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
17  *  GNU General Public License for more details.
18  *
19  *  You should have received a copy of the GNU General Public License
20  *  along with this program; if not, write to the Free Software
21  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
22  */
23
24 #include <ginac/ginac.h>
25
26 #ifndef NO_GINAC_NAMESPACE
27 using namespace GiNaC;
28 #endif // ndef NO_GINAC_NAMESPACE
29
30 unsigned lsolve_onedim(void)
31 {
32     unsigned result = 0;
33     symbol x("x");
34     ex eq, aux;
35     
36     cout << "checking linear solve..." << flush;
37     clog << "---------linear solve:" << endl;
38     
39     eq = (3*x+5 == numeric(8));
40     aux = lsolve(eq,x);
41     if (aux != 1) {
42         result++;
43         clog << "solution of 3*x+5==8 erroneously returned "
44              << aux << endl;
45     }
46     
47     if (! result) {
48         cout << " passed ";
49         clog << "(no output)" << endl;
50     } else {
51         cout << " failed ";
52     }
53     
54     return result;
55 }