Rank/Atop: Difference between revisions

From NARS2000
Jump to navigationJump to search
No edit summary
No edit summary
 
(4 intermediate revisions by 2 users not shown)
Line 1: Line 1:
==Monadic Derived Function==
==Rank Monadic Derived Function==


<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
Line 6: Line 6:
     <table border="0" cellpadding="5" cellspacing="0" summary="">
     <table border="0" cellpadding="5" cellspacing="0" summary="">
     <tr>
     <tr>
       <td valign="top"><apll>Z←(<i>f</i>⍤[X] Y) R</apll></td>
       <td valign="top"><apll>Z←<i>f</i>⍤[X] Y R</apll></td>
       <td></td>
       <td></td>
       <td></td>
       <td></td>
Line 36: Line 36:
</tr>
</tr>
<tr>
<tr>
   <td>The result <apll>Z</apll> is the <b>conforming</b> disclose of the above (see below).</td>
   <td>The result <apll>Z</apll> is the <b>conforming</b> disclose of the above, using <apll>[X]</apll>, if present.</td>
</tr>
</tr>
</table>
</table>
<br />
<br />
<p>The monadic derived function of this dyadic operator is partially implemented by calling the following internal magic function:</p>
<p>The monadic derived function of this dyadic operator is implemented by calling an internal magic function similar to the following:</p>


<apll>&nbsp;&nbsp;&nbsp;&nbsp;∇ Z←(LO #MonRank Y) R;YR<br />
<apll><pre>
[1]&nbsp;&nbsp;&nbsp;YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R<br />
    ∇ Z←(LO #MonRank[X] Y) R;YR
[2]&nbsp;&nbsp;&nbsp;:if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end<br />
[1]   YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R
[3]&nbsp;&nbsp;&nbsp;Z←LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →0<br />
[2]   if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[4]&nbsp;&nbsp;&nbsp;⎕PRO:YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R<br />
[3]   Z←LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →L1
[5]&nbsp;&nbsp;&nbsp;:if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end<br />
[4]   ⎕PRO:YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R
[6]&nbsp;&nbsp;&nbsp;Z←⊃LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R<br />
[5]   :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
&nbsp;&nbsp;&nbsp;&nbsp;∇</apll>
[6]   Z←⊃LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R
 
[7]  L1::if 0=⎕NC 'X' ⋄ Z←⊃Z ⋄ :else ⋄ Z←⊃[X] Z ⋄ :end
See the discussion below for details on the final processing of the result of this magic function.
    ∇</pre></apll>


<br />
<br />
<br />


==Dyadic Derived Function==
==Rank Dyadic Derived Function==


<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
Line 63: Line 62:
     <table border="0" cellpadding="5" cellspacing="0" summary="">
     <table border="0" cellpadding="5" cellspacing="0" summary="">
     <tr>
     <tr>
       <td valign="top"><apll>Z←L (<i>f</i>⍤[X] Y) R</apll></td>
       <td valign="top"><apll>Z←L <i>f</i>⍤[X] Y R</apll></td>
       <td></td>
       <td></td>
       <td></td>
       <td></td>
Line 101: Line 100:
</tr>
</tr>
<tr>
<tr>
   <td>The result <apll>Z</apll> is the <b>conforming</b> disclose of the above (see below).</td>
   <td>The result <apll>Z</apll> is the <b>conforming</b> disclose of the above, using <apll>[X]</apll>, if present.</td>
</tr>
</tr>
</table>
</table>
<br />
<br />
<p>The dyadic derived function of this dyadic operator is partially implemented by calling the following internal magic function:</p>
<p>The dyadic derived function of this dyadic operator is implemented by calling an internal magic function similar to the following:</p>
 
<apll><pre>
    ∇ Z←L (LO #DydRank[X] Y) R;YL;YR
[1]  (YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R
[2]  :if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end
[3]  :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[4]  Z←(⊂[(-YL)↑⍳⍴⍴L] L) LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →L1
[5]  ⎕PRO:(YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R
[6]  :if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end
[7]  :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[8]  Z←⊃(⊂[(-YL)↑⍳⍴⍴L]¨0⍴⊂L) LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R
[9]  L1::if 0=⎕NC 'X' ⋄ Z←⊃Z ⋄ :else ⋄ Z←⊃[X] Z ⋄ :end
    ∇</pre></apll>
 
For example, the rank operator can emulate laminate between the following two vectors to produce the first result, but not the second without help from the axis operator:


<apll>&nbsp;&nbsp;&nbsp;&nbsp;∇ Z←L (LO #DydRank Y) R;YL;YR<br />
<apll><pre>
[1]&nbsp;&nbsp;&nbsp;(YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R<br />
      L←'abcdef' {diamond} R←⍳⍴L
[2]&nbsp;&nbsp;&nbsp;:if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end<br />
      L ,⍤0 R                  a.k.a. L,[1.5] R
[3]&nbsp;&nbsp;&nbsp;:if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end<br />
a 1
[4]&nbsp;&nbsp;&nbsp;Z←(⊂[(-YL)↑⍳⍴⍴L] L) LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →0<br />
b 2
[5]&nbsp;&nbsp;&nbsp;⎕PRO:(YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R<br />
c 3
[6]&nbsp;&nbsp;&nbsp;:if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end<br />
d 4
[7]&nbsp;&nbsp;&nbsp;:if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end<br />
e 5
[8]&nbsp;&nbsp;&nbsp;Z←⊃(⊂[(-YL)↑⍳⍴⍴L]¨0⍴⊂L) LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R<br />
f 6
&nbsp;&nbsp;&nbsp;&nbsp;∇</apll>
      L ,⍤[1] 0 R              a.k.a. L,[0.5] R
a b c d e f
1 2 3 4 5 6</pre></apll>
 
==Atop Monadic Derived Function==
 
<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
<tr>
  <td>
    <table border="0" cellpadding="5" cellspacing="0" summary="">
    <tr>
      <td valign="top"><apll>Z←<i>f</i>⍤<i>g</i> R</apll></td>
      <td></td>
      <td></td>
      <td>Applies the monadic functions <apll><i>f</i></apll> and <apll><i>g</i></apll> to <apll>R</apll> as in <apll> f g R</apll>.</td>
    </tr>
    </table>
  </td>
</tr>
<tr>
  <td><apll>R</apll> is an arbitrary array, <apll><i>f</i></apll> and <apll><i>g</i></apll> are both arbitrary monadic functions.</td>
</tr>
</table>


== Conforming Disclose ==
<br />
<p>The monadic and dyadic derived functions of this dyadic operator are implemented by calling an internal magic function similar to the following:</p>


Both of the above magic functions implement their respective derived function except for some final processing which is essentially a disclose but one which allows for mismatched ranks.  That part is implemented by the following function where <apll>L</apll> has already been calculated as the maximum rank across all items:
<apll><pre>
#Atop←{⍺←⊢ ⋄ ⍺⍺ ⍺ ⍵⍵ ⍵}
</pre></apll>


<apll>&nbsp;&nbsp;&nbsp;&nbsp;∇ Z←L #Conform R</apll><br />
For example,
<apll>[1]&nbsp;&nbsp;&nbsp;Z←(((L-∊⍴∘⍴¨R)⍴¨1),¨⍴¨R)⍴¨R</apll><br />
<apll>&nbsp;&nbsp;&nbsp;&nbsp;∇</apll>


If the axis operator <apll>[X]</apll> is present on the rank operator, it is used in the final stage to disclose with axis (<apll>⊃[X]</apll>) the <apll>#Conform</apll> result to produce the final result <apll>Z</apll>. If the axis operator is not present, the final result <apll>Z</apll> is the disclose without axis (<apll>⊃</apll>) of the <apll>#Conform</apll> result.
<apll><pre>     ⌈⍤-¯1.1      ⍝ where f⍤g R ←→ f g R
2
</pre></apll>


For example, the rank operator can emulate laminate between the following two vectors to produce the first result, but not the second without help from the axis operator:
==Atop Dyadic Derived Function==
 
<table border="1" cellpadding="5" cellspacing="0" rules="none" summary="">
<tr>
  <td>
    <table border="0" cellpadding="5" cellspacing="0" summary="">
    <tr>
      <td valign="top"><apll>Z←L <i>f</i>⍤<i>g</i> R</apll></td>
      <td></td>
      <td></td>
      <td>Applies the monadic function <apll><i>f</i></apll> to the result of applying the dyadic function <apll><i>g</i></apll> between <apll>L</apll> and the <apll>R</apll>.  That is, it is equivalent to <apll><i>f</i> L <i>g</i> R</apll>.</td>
    </tr>
    </table>
  </td>
</tr>
<tr>
  <td><apll>L</apll> and <apll>R</apll> are arbitrary arrays, <apll><i>f</i></apll> is an arbitrary monadic function, <apll><i>g</i></apll> is an arbitrary dydadic function.</td>
</tr>
</table>
 
<br />
<p>The monadic and dyadic derived functions of this dyadic operator are implemented by calling an internal magic function similar to the following:</p>
 
<apll><pre>
#Atop←{⍺←⊢ ⋄ ⍺⍺ ⍺ ⍵⍵ ⍵}
</pre></apll>
 
For example,


<apll>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;L←'abcdef' {diamond} R←⍳⍴L</apll><br />
<apll><pre>     2⌈⍤-¯1.1      ⍝ where L f⍤g R ←→ f L g R
<apll>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;L ,⍤0 R</apll>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;a.k.a. <apll>L,[1.5] R</apll><br />
4
<apll>&nbsp;a 1</apll><br />
</pre></apll>
<apll>&nbsp;b 2</apll><br />
<apll>&nbsp;c 3</apll><br />
<apll>&nbsp;d 4</apll><br />
<apll>&nbsp;e 5</apll><br />
<apll>&nbsp;f 6</apll><br />
<apll>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;L ,⍤[1] 0 R</apll>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;a.k.a. <apll>L,[0.5] R</apll><br />
<apll>&nbsp;a b c d e f</apll><br />
<apll>&nbsp;1 2 3 4 5 6</apll><br />


== Acknowledgments ==
== Acknowledgments ==


The APL functions shown above were originally based upon those found in the paper J. Philip Benkhard, "[http://portal.acm.org/citation.cfm?id=114058&dl=GUIDE&coll=GUIDE&CFID=98347615&CFTOKEN=27795431 Extending structure, type, and expression in APL2]", ACM SIGAPL APL Quote Quad, v.21 n.4, p.20-38, Aug. 1991, but later was modified to reflect the definition of the Rank operator in the ISO-IEC 13751 Extended APL Standard.
The APL functions for the Rank operator shown above were originally based upon those found in the paper J. Philip Benkhard, "[http://portal.acm.org/citation.cfm?id=114058&dl=GUIDE&coll=GUIDE&CFID=98347615&CFTOKEN=27795431 Extending structure, type, and expression in APL2]", ACM SIGAPL APL Quote Quad, v.21 n.4, p.20-38, Aug. 1991, but later was modified to reflect the definition of the Rank operator in the ISO-IEC 13751 Extended APL Standard.
 
The Atop operator was defined by the developers at Dyalog APL.

Latest revision as of 11:11, 14 October 2023

Rank Monadic Derived Function

Z←f⍤[X] Y R Applies the monadic function f to the rank-r cells of R, where r is defined by Y and restores the cells to the result as per X.
R is an arbitrary array, f is an arbitrary monadic function, X is an origin-sensitive integer scalar or vector, and Y is an integer scalar or vector.
If 1<⍴⍴Y, signal a RANK ERROR.
If 1=⍴⍴Y and ~(⍴Y)∊1 2 3, signal a LENGTH ERROR.
Y is a scalar or one-, two-, or three-element integer vector that describes the cell size of the arguments. When it is a three-element vector, the first element specifies the cell size of the right argument when the derived function is called monadically; the second and third elements specify the left and right argument cell sizes when the derived function is called dyadically. When it is a two-element vector, the two values are used for the left and right arguments cell sizes when the derived function is called dyadically; if the derived function is called monadically, the cell size is the second value of the two. When it is a scalar or one-element vector, the single value is used for all three cell sizes. Moreover, the value of Y for the monadic derived function case may be positive or negative with an absolute value no larger than ⍴⍴R. In this case (first element only), the above rules may be summarized by setting Y to the one-element integer vector (-⍴⍴R)⌈(⍴⍴R)⌊1↑⌽3⍴⌽Y.
The cells from R are of rank |Y.
If Y is positive, the cells from R are taken from the right end of the shape vector; if negative, the cells are taken from the left end of the shape vector. In particular, CR←(-Y)↑⍴R is the shape of the righthand cells, and FR←(-Y)↓⍴R is the shape of the righthand frame.
Loop through the frame applying f to the shape CR cells from R.
The result Z is the conforming disclose of the above, using [X], if present.


The monadic derived function of this dyadic operator is implemented by calling an internal magic function similar to the following:

    ∇ Z←(LO #MonRank[X] Y) R;YR
[1]   YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R
[2]   if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[3]   Z←LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →L1
[4]   ⎕PRO:YR←(1↑⌽3⍴⌽Y)⌊⍴⍴R
[5]   :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[6]   Z←⊃LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R
[7]   L1::if 0=⎕NC 'X' ⋄ Z←⊃Z ⋄ :else ⋄ Z←⊃[X] Z ⋄ :end
    ∇


Rank Dyadic Derived Function

Z←L f⍤[X] Y R Applies the dyadic function f between the rank-l cells of L and the rank-r cells of R and restores the cells to the result as per X, where l and r are defined by Y.
L and R are arbitrary arrays, f is an arbitrary dyadic function, X is an origin-sensitive integer scalar or vector, and Y is an integer scalar or vector.
If 1<⍴⍴Y, signal a RANK ERROR.
If 1=⍴⍴Y and ~(⍴Y)∊1 2 3, signal a LENGTH ERROR.
Y is a scalar or one-, two-, or three-element integer vector that describes the cell size of the arguments. When it is a three-element vector, the first element specifies the cell size of the right argument when the derived function is called monadically; the second and third elements specify the left and right argument cell sizes when the derived function is called dyadically. When it is a two-element vector, the two values are used for the left and right arguments cell sizes when the derived function is called dyadically; if the derived function is called monadically, the cell size is the second value of the two. When it is a scalar or one-element vector, the single value is used for all three cell sizes. Moreover, the values of Y for the dyadic derived function case may be positive or negative with an absolute value no larger than (⍴⍴L),⍴⍴R. In this case (second and third elements only), the above rules may be summarized by setting Y to the two-element integer vector (-(⍴⍴L),⍴⍴R)⌈((⍴⍴L),⍴⍴R)⌊1↓⌽3⍴⌽Y.
The cells from L are of rank |1↑Y; the cells from R are of rank |1↓Y.
If 1↑Y is positive, the cells from L are taken from the right end of the shape vector; if negative, the cells are taken from the left end of the shape vector. In particular, CL←(-1↑Y)↑⍴L is the shape of the lefthand cells, and FL←(-1↑Y)↓⍴L is the shape of the lefthand frame.
If 1↓Y is positive, the cells from R are taken from the right end of the shape vector; if negative, the cells are taken from the left end of the shape vector. In particular, CR←(-1↓Y)↑⍴R is the shape of the righthand cells, and FR←(-1↓Y)↓⍴R is the shape of the righthand frame.
If FL and FR are both non-empty and (⍴FL)≠⍴FR, signal a RANK ERROR; if the shapes of FL and FR are the same, but their values differ, signal a LENGTH ERROR.
Loop through the frames (scalar extending as necessary) applying f between the shape CL cells from L and the shape CR cells from R.
The result Z is the conforming disclose of the above, using [X], if present.


The dyadic derived function of this dyadic operator is implemented by calling an internal magic function similar to the following:

    ∇ Z←L (LO #DydRank[X] Y) R;YL;YR
[1]   (YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R
[2]   :if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end
[3]   :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[4]   Z←(⊂[(-YL)↑⍳⍴⍴L] L) LO¨⊂[(-YR)↑⍳⍴⍴R] R ⋄ →L1
[5]   ⎕PRO:(YL YR)←(1↓⌽3⍴⌽Y)⌊(⍴⍴L),⍴⍴R
[6]   :if 0>YL ⋄ YL←0⌈YL+⍴⍴L ⋄ :end
[7]   :if 0>YR ⋄ YR←0⌈YR+⍴⍴R ⋄ :end
[8]   Z←⊃(⊂[(-YL)↑⍳⍴⍴L]¨0⍴⊂L) LO¨¨⊂[(-YR)↑⍳⍴⍴R]¨0⍴⊂R
[9]   L1::if 0=⎕NC 'X' ⋄ Z←⊃Z ⋄ :else ⋄ Z←⊃[X] Z ⋄ :end
    ∇

For example, the rank operator can emulate laminate between the following two vectors to produce the first result, but not the second without help from the axis operator:

      L←'abcdef' ⋄ R←⍳⍴L
      L ,⍤0 R                   a.k.a. L,[1.5] R
 a 1
 b 2
 c 3
 d 4
 e 5
 f 6
      L ,⍤[1] 0 R               a.k.a. L,[0.5] R
 a b c d e f
 1 2 3 4 5 6

Atop Monadic Derived Function

Z←fg R Applies the monadic functions f and g to R as in f g R.
R is an arbitrary array, f and g are both arbitrary monadic functions.


The monadic and dyadic derived functions of this dyadic operator are implemented by calling an internal magic function similar to the following:

#Atop←{⍺←⊢ ⋄ ⍺⍺ ⍺ ⍵⍵ ⍵}

For example,

      ⌈⍤-¯1.1      ⍝ where f⍤g R ←→ f g R
2

Atop Dyadic Derived Function

Z←L fg R Applies the monadic function f to the result of applying the dyadic function g between L and the R. That is, it is equivalent to f L g R.
L and R are arbitrary arrays, f is an arbitrary monadic function, g is an arbitrary dydadic function.


The monadic and dyadic derived functions of this dyadic operator are implemented by calling an internal magic function similar to the following:

#Atop←{⍺←⊢ ⋄ ⍺⍺ ⍺ ⍵⍵ ⍵}

For example,

      2⌈⍤-¯1.1      ⍝ where L f⍤g R ←→ f L g R
4

Acknowledgments

The APL functions for the Rank operator shown above were originally based upon those found in the paper J. Philip Benkhard, "Extending structure, type, and expression in APL2", ACM SIGAPL APL Quote Quad, v.21 n.4, p.20-38, Aug. 1991, but later was modified to reflect the definition of the Rank operator in the ISO-IEC 13751 Extended APL Standard.

The Atop operator was defined by the developers at Dyalog APL.