Annotation of OpenXM/src/kan96xx/Doc/factor-a.sm1, Revision 1.1.1.1
1.1 maekawa 1: %% This package requires kan/sm1 version 2.981101 or later.
2:
3: /factor-a.version (2.981101) def
4: (factor-a.sm1 : kan/sm1 package to factor polynomials by calling risa/asir.)
5: message
6: ( : kan/sm1 package to simplify rationals by calling risa/asir.)
7: message
8: ( Version : ) messagen factor-a.version message
9:
10: factor-a.version [(Version)] system_variable gt
11: { (This package requires the latest version of kan/sm1) message
12: (Please get it from http://www.math.kobe-u.ac.jp/KAN) message
13: error
14: } { } ifelse
15:
16:
17: [(factor)
18: [(polynomial factor list_of_strings)
19: (Example: (x^2-1). factor :: ---> [[$1$ $1$] [$x-1$ $1$] [$x+1$ $1$]])
20: (cf.: data_conversion, map, get, pushfile)
21: (Note: The function call creates work files asir-tmp.t, asir-tmp.tt,)
22: ( asir-tmp-out.t, asri-tmp-log.t and asir-tmp-out.tt )
23: ( in the current directory.)
24: ]
25: ] putUsages
26:
27: %% /f (Dx^10*d*a-d*a) def
28:
29: /factor-asir-1 {
30: /arg1 set
31: [/f /fd /fnewline] pushVariables
32: [
33: arg1 /f set
34: %% (factor-asir-1 is tested with Asir version 950831 on Linux.) message
35: (asir-tmp.t) (w) file /fd set
36: /fnewline { fd 10 (string) data_conversion writestring } def
37: fd $output("asir-tmp-out.t");$ writestring fnewline
38: fd $fctr($ writestring
39: fd f writestring
40: fd $); output(); quit(); $ writestring fnewline
41: fd closefile
42: (/bin/rm -f asir-tmp.tt) system
43: (sed "s/D/kanD/g" asir-tmp.t | sed "s/E/kanE/g" | sed "s/Q/kanQ/g" >asir-tmp.tt) system
44: (/bin/rm -f asir-tmp-out.t asir-tmp-out.tt asir-tmp-log.t) system
45: (asir <asir-tmp.tt >asir-tmp-log.t) system
46: (sed "s/\[1\]/ /g" asir-tmp-out.t | sed "s/\[2\]/ /g" | sed "1s/1/ /g"| sed "s/\[/{/g" | sed "s/\]/}/g" | sed "s/kanD/D/g" | sed "s/kanE/E/g" | sed "s/kanQ/Q/g" >asir-tmp-out.tt) system
47: ] pop
48: popVariables
49: } def
50:
51: /clean-workfiles {
52: (/bin/rm -f asir-tmp-out.t asir-tmp-out.tt asir-tmp.t asir-tmp.tt sm0-tmp.t sm0-tmp-out.t asir-tmp-log.t) system
53: } def
54:
55:
56: %% comment: there is not data conversion function from string --> array
57: %% e.g. (abc) ---> [0x61, 0x62, 0x63]
58: %% We can do (abc) 1 10 put, but "get" does not work for strings.
59:
60: %% f factor-asir-1
61:
62: %%/aaa
63: %% ({{1,1},{x-1,1},{x+1,1},{x^4+x^3+x^2+x+1,1},{x^4-x^3+x^2-x+1,1}})
64: %%def
65:
66: /asir-list-to-kan {
67: /arg1 set
68: [/aaa /ftmp /ftmp2] pushVariables
69: [
70: /aaa arg1 def
71: [ aaa to_records pop ] /ftmp set
72: ftmp { to_records pop [ 3 1 roll ] } map /ftmp2 set
73: /arg1 ftmp2 def
74: ] pop
75: popVariables
76: arg1
77: } def
78:
79: /foo {
80: (input string is in f) message
81: f ::
82: f factor-asir-1
83: %% (asir-tmp-out.tt) run
84: %% (answer in @asir.out) message
85: %% bug of run.
86: (asir-tmp-out.tt) pushfile /@asir.out set
87: @asir.out asir-list-to-kan /ff2 set
88: (answer in ff2) message
89: } def
90:
91: /factor {
92: (string) data_conversion
93: factor-asir-1
94: (asir-tmp-out.tt) pushfile asir-list-to-kan
95: } def
96:
97: %%%%%%%%%%%%%%%%% macros for simplification (reduction, cancel)
98: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
99: [(cancel)
100: [(polynomial cancel list_of_strings)
101: (This function simplifies rationals.)
102: (Example: $x^2-1$. $x+1$. div cancel :: ---> [[$x-1$ , $1$]])
103: (Note: The function call creates work files asir-tmp.t, asir-tmp.tt,)
104: ( asir-tmp-out.t, asri-tmp-log.t and asir-tmp-out.tt )
105: ( in the current directory.)
106: ]
107: ] putUsages
108:
109: /reduce-asir-1 {
110: /arg1 set
111: [/f /fd /fnewline] pushVariables
112: [
113: arg1 /f set
114: %% (reduce-asir-1 is tested with Asir version 950831 on Linux.) message
115: (asir-tmp.t) (w) file /fd set
116: /fnewline { fd 10 (string) data_conversion writestring } def
117: fd $output("asir-tmp-out.t");$ writestring fnewline
118: fd $AsirTmp012=red($ writestring
119: fd f writestring
120: fd $)$ writestring
121: fd ($ ) writestring fnewline
122: fd $AsirTmp013=ptozp(nm(AsirTmp012))$ writestring
123: fd ($ ) writestring fnewline
124: fd $AsirTmp014=red(nm(AsirTmp012)/AsirTmp013)$ writestring
125: fd ($ ) writestring fnewline
126: fd $[[nm(AsirTmp014)*AsirTmp013,dn(AsirTmp014)*dn(AsirTmp012)]];output();quit(); $ writestring fnewline
127: fd closefile
128: (/bin/rm -f asir-tmp.tt) system
129: (sed "s/D/kanD/g" asir-tmp.t | sed "s/E/kanE/g" | sed "s/Q/kanQ/g" >asir-tmp.tt) system
130: (/bin/rm -f asir-tmp-out.t asir-tmp-out.tt asir-tmp-log.t) system
131: (asir <asir-tmp.tt >asir-tmp-log.t) system
132: (sed "s/\[1\]/ /g" asir-tmp-out.t | sed "s/\[2\]/ /g" |sed "s/\[3\]/ /g" |sed "s/\[4\]/ /g" |sed "s/\[5\]/ /g" | sed "1s/1/ /g"| sed "s/\[/{/g" | sed "s/\]/}/g" | sed "s/kanD/D/g" | sed "s/kanE/E/g" | sed "s/kanQ/Q/g" | sed "s/kanLeftBrace/\[/g" | sed "s/kanRightBrace/\]/g" | sed "s/kanComma/\,/g" >asir-tmp-out.tt) system
133:
134:
135: ] pop
136: popVariables
137: } def
138:
139: /cancel {
140: (string) data_conversion
141: reduce-asir-1
142: (asir-tmp-out.tt) pushfile asir-list-to-kan
143: } def
144:
145:
146: (Loaded macros "factor", "cancel") message
147:
148:
149:
150:
151:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to factor-a.sm1 CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM / src / kan96xx / Doc