jengelhdatetimespielwiese
Last change
on this file since 8930f2c was
8930f2c,
checked in by Hans Schoenemann <hannes@…>, 13 years ago

surfacesignature tests
gitsvnid: file:///usr/local/Singular/svn/trunk@13375 2c84dea37e6841379b89c4e89433aadc


Property mode set to
100755

File size:
1.3 KB

Line  

1  LIB "tst.lib"; 

2  tst_init(); 

3  

4  LIB "surfacesignature.lib"; 

5  

6  example brieskornSign; 

7  

8  ring R = 0, (x,y), dp; 

9  poly f, g, g1, g2, g3; 

10  int N; 

11  number a,b,c; 

12  

13  N = 5; 

14  f = x1521x14+8x13y6x1316x12y+20x11y2x12+8x11y36x10y2 

15  +24x9y3+4x9y216x8y3+26x7y46x6y4+8x5y5+4x3y6y8; 

16  a = signature(N,f,1); 

17  b = signature(N,f,2); 

18  a; 

19  a == b; 

20  

21  g1 = f^3 + x9y8; 

22  signature(N,g1,2); 

23  

24  g3 = f^5 + x20y20; 

25  signature(N,g3,2); 

26  

27  

28  N = 6; 

29  f = y4+2x3y2+x6+x5y; 

30  a = signature(N,f,1); 

31  b = signature(N,f,2); 

32  c = signature(N,f,3); 

33  a; 

34  a == b; 

35  b == c; 

36  

37  g1 = f^2 + x5y5; 

38  a = signature(N,g1,1); 

39  b = signature(N,g1,2); 

40  a; 

41  a == b; 

42  

43  g2 = f^3 + x11y11; 

44  a = signature(N,g2,1); 

45  b = signature(N,g2,2); 

46  a; 

47  a == b; 

48  

49  

50  N = 7; 

51  f = x5+y11; 

52  a = signature(N,f,1); 

53  b = signature(N,f,2); 

54  c = signature(N,f,3); 

55  a; 

56  a == b; 

57  b == c; 

58  

59  g1 = f^2 + x5y5; 

60  b = signature(N,g1,2); 

61  c = signature(N,g1,3); 

62  b; 

63  b == c; 

64  

65  g2 = f^3 + x11y11; 

66  a = signature(N,g2,1); 

67  b = signature(N,g2,2); 

68  a; 

69  a == b; 

70  

71  

72  N = 6; 

73  f = x71+6x65+15x59630x52y6+20x53+6230x46y6+910x39y12+15x47 

74  7530x40y6+14955x33y12285x26y18+6x41+1230x34y6+4680x27y12 

75  +1830x20y18+30x13y24+x355x28y6+10x21y1210x14y18+5x7y24y30; 

76  a = signature(N,f,1); 

77  b = signature(N,f,2); 

78  a; 

79  a == b; 

80  

81  

82  N = 12; 

83  f = 2x3y32y5+x4xy2; 

84  b = signature(N,f,2); 

85  c = signature(N,f,3); 

86  b; 

87  b == c; 

88  

89  f = x3y3+x6y+xy6x4y4; 

90  b = signature(N,f,2); 

91  c = signature(N,f,3); 

92  b; 

93  b == c; 

94  

95  tst_status(1);$ 

Note: See
TracBrowser
for help on using the repository browser.