nag/parse_test.go at master · fumin/nag · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
package nag

import (
	"fmt"
	"testing"
)

func TestParse(t *testing.T) {
	tests := []struct {
		variables map[string]Symbol
		order     Order
		input     string
		p         *Polynomial[*Rat]
	}{
		{
			variables: map[string]Symbol{"a": 1, "b": 2},
			order:     Deglex,
			input:     "ba^3",
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(1, 1), Monomial: Monomial{2, 1, 1, 1}},
			),
		},
		{
			variables: map[string]Symbol{"a": 1, `{b^{\dagger}}`: 2},
			order:     Deglex,
			input:     `-{b^{\dagger}}^1a^3`,
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{2, 1, 1, 1}},
			),
		},
		{
			variables: map[string]Symbol{"a": 1, "b": 2},
			order:     Deglex,
			input:     "(-ba)^3a-a(3b-5b^2)",
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{2, 1, 2, 1, 2, 1, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-3, 1), Monomial: Monomial{1, 2}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(5, 1), Monomial: Monomial{1, 2, 2}},
			),
		},
		{
			variables: map[string]Symbol{"a": 1, "b": 2},
			order:     Deglex,
			input:     "(a-b)^3",
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(1, 1), Monomial: Monomial{1, 1, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{1, 1, 2}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{1, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(1, 1), Monomial: Monomial{1, 2, 2}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{2, 1, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(1, 1), Monomial: Monomial{2, 1, 2}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(1, 1), Monomial: Monomial{2, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-1, 1), Monomial: Monomial{2, 2, 2}},
			),
		},
		{
			variables: map[string]Symbol{"a": 1, "b": 2, "c": 3},
			order:     Deglex,
			input:     "-12/5a^3((a+cc)b)^2a+7/3ca-3/2b",
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(-12, 5), Monomial: Monomial{1, 1, 1, 3, 3, 2, 3, 3, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-12, 5), Monomial: Monomial{1, 1, 1, 1, 2, 3, 3, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-12, 5), Monomial: Monomial{1, 1, 1, 3, 3, 2, 1, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-12, 5), Monomial: Monomial{1, 1, 1, 1, 2, 1, 2, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(7, 3), Monomial: Monomial{3, 1}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(-3, 2), Monomial: Monomial{2}},
			),
		},
		{
			variables: map[string]Symbol{"a": 1, "b": 2, "c": 3},
			order:     Deglex,
			input:     "5/3b(a+b)^2c+9a",
			p: NewPolynomial(
				NewRat(0, 1), Deglex,
				PolynomialTerm[*Rat]{Coefficient: NewRat(5, 3), Monomial: Monomial{2, 1, 1, 3}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(5, 3), Monomial: Monomial{2, 1, 2, 3}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(5, 3), Monomial: Monomial{2, 2, 1, 3}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(5, 3), Monomial: Monomial{2, 2, 2, 3}},
				PolynomialTerm[*Rat]{Coefficient: NewRat(9, 1), Monomial: Monomial{1}},
			),
		},
	}

	for i, test := range tests {
		t.Run(fmt.Sprintf("%d", i), func(t *testing.T) {
			t.Parallel()
			p, err := Parse(test.variables, test.order, test.input)
			if err != nil {
				t.Fatalf("%+v", err)
			}
			if !p.Equal(test.p) {
				t.Errorf("%v %v", p, test.p)
			}
		})
	}
}