1 | /* -*- C++ -*- |
---|
2 | * |
---|
3 | * This file is a part of LEMON, a generic C++ optimization library |
---|
4 | * |
---|
5 | * Copyright (C) 2003-2006 |
---|
6 | * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
---|
7 | * (Egervary Research Group on Combinatorial Optimization, EGRES). |
---|
8 | * |
---|
9 | * Permission to use, modify and distribute this software is granted |
---|
10 | * provided that this copyright notice appears in all copies. For |
---|
11 | * precise terms see the accompanying LICENSE file. |
---|
12 | * |
---|
13 | * This software is provided "AS IS" with no warranty of any kind, |
---|
14 | * express or implied, and with no claim as to its suitability for any |
---|
15 | * purpose. |
---|
16 | * |
---|
17 | */ |
---|
18 | |
---|
19 | #include <lemon/polynomial.h> |
---|
20 | #include <lemon/xy.h> |
---|
21 | #include <iostream> |
---|
22 | #include "test_tools.h" |
---|
23 | |
---|
24 | using namespace std; |
---|
25 | using namespace lemon; |
---|
26 | int main() |
---|
27 | { |
---|
28 | Polynomial<int> pi(5); |
---|
29 | check(pi.deg()==5,"Something is wrong here."); |
---|
30 | pi[0]=12; |
---|
31 | pi[5]=3; |
---|
32 | pi[2]=7; |
---|
33 | check(pi[1]==0,"Uninitialized elements should be(?) zero."); |
---|
34 | { |
---|
35 | Polynomial<double> pd=pi; |
---|
36 | check(pd.deg()==5,"Something is wrong here."); |
---|
37 | check(pd[0]==12,"Something is wrong here."); |
---|
38 | check(pd[5]==3,"Something is wrong here."); |
---|
39 | |
---|
40 | } |
---|
41 | |
---|
42 | Polynomial<double> pd; |
---|
43 | pd=pi; |
---|
44 | check(pd.deg()==5,"Something is wrong here."); |
---|
45 | check(pd[0]==12,"Something is wrong here."); |
---|
46 | check(pd[5]==3,"Something is wrong here."); |
---|
47 | |
---|
48 | check(pd(0)==12,"Something is wrong here."); |
---|
49 | check(pd(1)==22,"Something is wrong here."); |
---|
50 | check(pd(2)==136,"Something is wrong here."); |
---|
51 | |
---|
52 | check((pd*pi).deg()==10,"Something is wrong here."); |
---|
53 | check((pd*pi)[10]==9,"Something is wrong here."); |
---|
54 | check((pd*pi)[7]==42,"Something is wrong here."); |
---|
55 | |
---|
56 | Polynomial<double> pd2=pd+pi; |
---|
57 | check(pd2[5]==6,"Something is wrong here."); |
---|
58 | pd2+=pd; |
---|
59 | pd2+=pi; |
---|
60 | check(pd2[5]==12,"Something is wrong here."); |
---|
61 | pd2-=pd; |
---|
62 | pd2-=pi; |
---|
63 | check(pd2[5]==6,"Something is wrong here."); |
---|
64 | check((pd-pi)[5]==0,"Something is wrong here."); |
---|
65 | |
---|
66 | Polynomial<double> pdd=pd.derivate(); |
---|
67 | pd.derivateMyself(); |
---|
68 | check(pdd==pd,"Something is wrong here."); |
---|
69 | check(pd.deg()==4,"Something is wrong here."); |
---|
70 | check(pd[4]==15,"Something is wrong here."); |
---|
71 | |
---|
72 | |
---|
73 | Polynomial<double> pdi=pd.integrate(); |
---|
74 | pd.integrateMyself(); |
---|
75 | check(pdi==pd,"Something is wrong here."); |
---|
76 | check(pd.deg()==5,"Something is wrong here."); |
---|
77 | check(pd[5]==3,"Something is wrong here."); |
---|
78 | check(pd[0]==0,"Something is wrong here."); |
---|
79 | |
---|
80 | |
---|
81 | } |
---|