URL https://opencores.org/ocsvn/eco32/eco32/trunk

Subversion Repositories eco32

[/] [eco32/] [trunk/] [fp/] [implementation/] [mmix/] [saddle2.mms] - Blame information for rev 133

Go to most recent revision | Details | Compare with Previous | View Log


* Exercise 1.3.2'--18, Solution 2
      LOC #100
t IS $255
a00 GREG Data_Segment
a10 GREG Data_Segment+8
a20 GREG Data_Segment+8*2
ij  GREG  % element index
ii  GREG  % row index times 8
j   GREG  % column index
x   GREG  % current maximum
y   GREG  % current element
z   GREG  % current min max
ans IS $0 % return register
Phase1 SET j,8 Start at column 8.
       SET z,1000  $\.z\gets\infty$ (more or less).
3H     ADD ij,j,9*8-2*8
       LDB x,a20,ij
1H     LDB y,a10,ij
       CMP t,x,y    Is x
       CSN x,t,y   If so, update the maximum.
2H     SUB ij,ij,8  Move up one.
       PBP ij,1B
       STB x,a10,ij Store column maximum.
       CMP t,x,z    Is x
       CSN z,t,x   If so, update the min max.
       SUB  j,j,1   Move left a column.
       PBP  j,3B
Phase2 SET  ii,9*8-8 At this point $\.z=\min_jC(j)$
3H     ADD  ij,ii,8  Prepare to search a row.
       SET  j,8
1H     LDB  x,a10,ij
       SUB  t,z,x    Is $\.z>a_{ij}$?
       PBP  t,No     No saddle in this row
       PBN  t,2F
       LDB  x,a00,j  Is $a_{ij}=C(j)$?
       CMP  t,x,z
       CSZ  ans,t,ij If so, remember a possible saddle point.
2H     SUB  j,j,1    Move left in row.
       SUB  ij,ij,1
       PBP  j,1B
       LDA  ans,a10,ans  A saddle point was found here.
       POP  1,0
No     SUB  ii,ii,8
       PBP  ii,3B     Try another row.
       SET  ans,0
       POP  1,0     $\.{ans} = 0$; no saddle.\quad\slug
 
aaaa   GREG 6364136223846793005  C E Haynes's multiplier
Main   SET  ij,9*8      assume that $1 = seed
1H     MULU $1,$1,aaaa
       INCL $1,1
       MULU x,$1,5
       GET  x,rH
       SUB  x,x,2
       STB  x,a10,ij
       SUB  ij,ij,1
       PBP  ij,1B
       PUSHJ 2,Phase1
       JMP  Main

Line No.	Rev	Author	Line
1	15	hellwig	`* Exercise 1.3.2'--18, Solution 2`
2			`LOC #100`
3			`t IS $255`
4			`a00 GREG Data_Segment`
5			`a10 GREG Data_Segment+8`
6			`a20 GREG Data_Segment+8*2`
7			`ij GREG % element index`
8			`ii GREG % row index times 8`
9			`j GREG % column index`
10			`x GREG % current maximum`
11			`y GREG % current element`
12			`z GREG % current min max`
13			`ans IS $0 % return register`
14			`Phase1 SET j,8 Start at column 8.`
15			`SET z,1000 $\.z\gets\infty$ (more or less).`
16			`3H ADD ij,j,98-28`
17			`LDB x,a20,ij`
18			`1H LDB y,a10,ij`
19			`CMP t,x,y Is x`
20			`CSN x,t,y If so, update the maximum.`
21			`2H SUB ij,ij,8 Move up one.`
22			`PBP ij,1B`
23			`STB x,a10,ij Store column maximum.`
24			`CMP t,x,z Is x`
25			`CSN z,t,x If so, update the min max.`
26			`SUB j,j,1 Move left a column.`
27			`PBP j,3B`
28			`Phase2 SET ii,9*8-8 At this point $\.z=\min_jC(j)$`
29			`3H ADD ij,ii,8 Prepare to search a row.`
30			`SET j,8`
31			`1H LDB x,a10,ij`
32			`SUB t,z,x Is $\.z>a_{ij}$?`
33			`PBP t,No No saddle in this row`
34			`PBN t,2F`
35			`LDB x,a00,j Is $a_{ij}=C(j)$?`
36			`CMP t,x,z`
37			`CSZ ans,t,ij If so, remember a possible saddle point.`
38			`2H SUB j,j,1 Move left in row.`
39			`SUB ij,ij,1`
40			`PBP j,1B`
41			`LDA ans,a10,ans A saddle point was found here.`
42			`POP 1,0`
43			`No SUB ii,ii,8`
44			`PBP ii,3B Try another row.`
45			`SET ans,0`
46			`POP 1,0 $\.{ans} = 0$; no saddle.\quad\slug`
47
48			`aaaa GREG 6364136223846793005 C E Haynes's multiplier`
49			`Main SET ij,9*8 assume that $1 = seed`
50			`1H MULU $1,$1,aaaa`
51			`INCL $1,1`
52			`MULU x,$1,5`
53			`GET x,rH`
54			`SUB x,x,2`
55			`STB x,a10,ij`
56			`SUB ij,ij,1`
57			`PBP ij,1B`
58			`PUSHJ 2,Phase1`
59			`JMP Main`

Browse

Tools

Subversion Repositories eco32

[/] [eco32/] [trunk/] [fp/] [implementation/] [mmix/] [saddle2.mms] - Blame information for rev 133