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
|
# AMAPR -- Linear transformation, type real. The range of pixel values
# A1 to A2 are mapped into the range B1 to B2 using a linear transformation.
# Values less than A1 or greater than A2 are mapped into the values B1 and
# B2 upon output.
.data 0
.set A, 4
.set B, 8
.set NPIX, 12
.set A1, 16
.set A2, 20
.set B1, 24
.set B2, 28
.align 2
.globl _amapr_
.text
# AMAPR (a, b, npix, a1, a2, b1, b2)
#
# scalar = real (b2 - b1) / real (a2 - a1)
# minout = min (b1, b2)
# maxout = max (b1, b2)
#
# do i = 1, npix
# b[i] = max(minout, min(maxout,
# PIXEL((a[i] - a1) * scalar) + b1))
#
# Registers:
# r0 last_a
# r1 a
# r2 b
# r3 scalar
# r4 minout
# r5 maxout
# r6 a1
# r7 b1
_amapr_:
.word 01774 # save r2-r9
movl A(ap), r1
movl B(ap), r2
mull3 $4, *NPIX(ap), r0
addl2 r1, r0
movf *A1(ap), r6
movf *B1(ap), r7
movf *A2(ap), r8
movf *B2(ap), r9
subf3 r7, r9, r3 # r3 = (b2 - b1) / (a2 - a1)
subf3 r6, r8, r4
divf2 r4, r3
cmpf r7, r9 # b1 <= b2
bleq L1
movf r9, r4 # no, min=b2, max=b1
movf r7, r5
jbr L2
L1: movf r7, r4 # yes, min=b1, max=b2
movf r9, r5
L2:
subf3 r6, (r1)+, r8 # r8 = a[i] - a1
mulf2 r3, r8 # (..) * scalar
addf2 r7, r8 # (..) + b1
cmpf r8, r4 # r8 < minout?
bgtr L3
movf r4, (r2)+
jbr L5
L3:
cmpf r8, r5 # r8 > maxout?
blss L4
movf r5, (r2)+
jbr L5
L4:
movf r8, (r2)+ # new value in range
L5:
cmpl r1, r0
blssu L2 # loop again
ret
|
