|
| 1 | +# RV64IDV system |
| 2 | +# |
| 3 | +# void |
| 4 | +# sgemm_nn(size_t n, |
| 5 | +# size_t m, |
| 6 | +# size_t k, |
| 7 | +# const float*a, // m * k matrix |
| 8 | +# size_t lda, |
| 9 | +# const float*b, // k * n matrix |
| 10 | +# size_t ldb, |
| 11 | +# float*c, // m * n matrix |
| 12 | +# size_t ldc) |
| 13 | +# |
| 14 | +# c += a*b (alpha=1, no transpose on input matrices) |
| 15 | +# matrices stored in C row-major order |
| 16 | + |
| 17 | +#define n a0 |
| 18 | +#define m a1 |
| 19 | +#define k a2 |
| 20 | +#define ap a3 |
| 21 | +#define astride a4 |
| 22 | +#define bp a5 |
| 23 | +#define bstride a6 |
| 24 | +#define cp a7 |
| 25 | +#define cstride t0 |
| 26 | +#define kt t1 |
| 27 | +#define nt t2 |
| 28 | +#define bnp t3 |
| 29 | +#define cnp t4 |
| 30 | +#define akp t5 |
| 31 | +#define bkp s0 |
| 32 | +#define nvl s1 |
| 33 | +#define ccp s2 |
| 34 | +#define amp s3 |
| 35 | + |
| 36 | +# Use args as additional temporaries |
| 37 | +#define ft12 fa0 |
| 38 | +#define ft13 fa1 |
| 39 | +#define ft14 fa2 |
| 40 | +#define ft15 fa3 |
| 41 | + |
| 42 | +# This version holds a 16*VLMAX block of C matrix in vector registers |
| 43 | +# in inner loop, but otherwise does not cache or TLB tiling. |
| 44 | + |
| 45 | +sgemm_nn: |
| 46 | + addi sp, sp, -FRAMESIZE |
| 47 | + sd s0, OFFSET(sp) |
| 48 | + sd s1, OFFSET(sp) |
| 49 | + sd s2, OFFSET(sp) |
| 50 | + |
| 51 | + # Check for zero size matrices |
| 52 | + beqz n, exit |
| 53 | + beqz m, exit |
| 54 | + beqz k, exit |
| 55 | + |
| 56 | + # Convert elements strides to byte strides. |
| 57 | + ld cstride, OFFSET(sp) # Get arg from stack frame |
| 58 | + slli astride, astride, 2 |
| 59 | + slli bstride, bstride, 2 |
| 60 | + slli cstride, cstride, 2 |
| 61 | + |
| 62 | + slti t6, m, 16 |
| 63 | + bnez t6, end_rows |
| 64 | + |
| 65 | +c_row_loop: # Loop across rows of C blocks |
| 66 | + |
| 67 | + mv nt, n # Initialize n counter for next row of C blocks |
| 68 | + |
| 69 | + mv bnp, bp # Initialize B n-loop pointer to start |
| 70 | + mv cnp, cp # Initialize C n-loop pointer |
| 71 | + |
| 72 | +c_col_loop: # Loop across one row of C blocks |
| 73 | + vsetvli nvl, nt, e32, ta, ma # 32-bit vectors, LMUL=1 |
| 74 | + |
| 75 | + mv akp, ap # reset pointer into A to beginning |
| 76 | + mv bkp, bnp # step to next column in B matrix |
| 77 | + |
| 78 | + # Initalize current C submatrix block from memory. |
| 79 | + vle32.v v0, (cnp); add ccp, cnp, cstride; |
| 80 | + vle32.v v1, (ccp); add ccp, ccp, cstride; |
| 81 | + vle32.v v2, (ccp); add ccp, ccp, cstride; |
| 82 | + vle32.v v3, (ccp); add ccp, ccp, cstride; |
| 83 | + vle32.v v4, (ccp); add ccp, ccp, cstride; |
| 84 | + vle32.v v5, (ccp); add ccp, ccp, cstride; |
| 85 | + vle32.v v6, (ccp); add ccp, ccp, cstride; |
| 86 | + vle32.v v7, (ccp); add ccp, ccp, cstride; |
| 87 | + vle32.v v8, (ccp); add ccp, ccp, cstride; |
| 88 | + vle32.v v9, (ccp); add ccp, ccp, cstride; |
| 89 | + vle32.v v10, (ccp); add ccp, ccp, cstride; |
| 90 | + vle32.v v11, (ccp); add ccp, ccp, cstride; |
| 91 | + vle32.v v12, (ccp); add ccp, ccp, cstride; |
| 92 | + vle32.v v13, (ccp); add ccp, ccp, cstride; |
| 93 | + vle32.v v14, (ccp); add ccp, ccp, cstride; |
| 94 | + vle32.v v15, (ccp) |
| 95 | + |
| 96 | + |
| 97 | + mv kt, k # Initialize inner loop counter |
| 98 | + |
| 99 | + # Inner loop scheduled assuming 4-clock occupancy of vfmacc instruction and single-issue pipeline |
| 100 | + # Software pipeline loads |
| 101 | + flw ft0, (akp); add amp, akp, astride; |
| 102 | + flw ft1, (amp); add amp, amp, astride; |
| 103 | + flw ft2, (amp); add amp, amp, astride; |
| 104 | + flw ft3, (amp); add amp, amp, astride; |
| 105 | + # Get vector from B matrix |
| 106 | + vle32.v v16, (bkp) |
| 107 | + |
| 108 | + # Loop on inner dimension for current C block |
| 109 | + k_loop: |
| 110 | + vfmacc.vf v0, ft0, v16 |
| 111 | + add bkp, bkp, bstride |
| 112 | + flw ft4, (amp) |
| 113 | + add amp, amp, astride |
| 114 | + vfmacc.vf v1, ft1, v16 |
| 115 | + addi kt, kt, -1 # Decrement k counter |
| 116 | + flw ft5, (amp) |
| 117 | + add amp, amp, astride |
| 118 | + vfmacc.vf v2, ft2, v16 |
| 119 | + flw ft6, (amp) |
| 120 | + add amp, amp, astride |
| 121 | + flw ft7, (amp) |
| 122 | + vfmacc.vf v3, ft3, v16 |
| 123 | + add amp, amp, astride |
| 124 | + flw ft8, (amp) |
| 125 | + add amp, amp, astride |
| 126 | + vfmacc.vf v4, ft4, v16 |
| 127 | + flw ft9, (amp) |
| 128 | + add amp, amp, astride |
| 129 | + vfmacc.vf v5, ft5, v16 |
| 130 | + flw ft10, (amp) |
| 131 | + add amp, amp, astride |
| 132 | + vfmacc.vf v6, ft6, v16 |
| 133 | + flw ft11, (amp) |
| 134 | + add amp, amp, astride |
| 135 | + vfmacc.vf v7, ft7, v16 |
| 136 | + flw ft12, (amp) |
| 137 | + add amp, amp, astride |
| 138 | + vfmacc.vf v8, ft8, v16 |
| 139 | + flw ft13, (amp) |
| 140 | + add amp, amp, astride |
| 141 | + vfmacc.vf v9, ft9, v16 |
| 142 | + flw ft14, (amp) |
| 143 | + add amp, amp, astride |
| 144 | + vfmacc.vf v10, ft10, v16 |
| 145 | + flw ft15, (amp) |
| 146 | + add amp, amp, astride |
| 147 | + addi akp, akp, 4 # Move to next column of a |
| 148 | + vfmacc.vf v11, ft11, v16 |
| 149 | + beqz kt, 1f # Don't load past end of matrix |
| 150 | + flw ft0, (akp) |
| 151 | + add amp, akp, astride |
| 152 | +1: vfmacc.vf v12, ft12, v16 |
| 153 | + beqz kt, 1f |
| 154 | + flw ft1, (amp) |
| 155 | + add amp, amp, astride |
| 156 | +1: vfmacc.vf v13, ft13, v16 |
| 157 | + beqz kt, 1f |
| 158 | + flw ft2, (amp) |
| 159 | + add amp, amp, astride |
| 160 | +1: vfmacc.vf v14, ft14, v16 |
| 161 | + beqz kt, 1f # Exit out of loop |
| 162 | + flw ft3, (amp) |
| 163 | + add amp, amp, astride |
| 164 | + vfmacc.vf v15, ft15, v16 |
| 165 | + vle32.v v16, (bkp) # Get next vector from B matrix, overlap loads with jump stalls |
| 166 | + j k_loop |
| 167 | + |
| 168 | +1: vfmacc.vf v15, ft15, v16 |
| 169 | + |
| 170 | + # Save C matrix block back to memory |
| 171 | + vse32.v v0, (cnp); add ccp, cnp, cstride; |
| 172 | + vse32.v v1, (ccp); add ccp, ccp, cstride; |
| 173 | + vse32.v v2, (ccp); add ccp, ccp, cstride; |
| 174 | + vse32.v v3, (ccp); add ccp, ccp, cstride; |
| 175 | + vse32.v v4, (ccp); add ccp, ccp, cstride; |
| 176 | + vse32.v v5, (ccp); add ccp, ccp, cstride; |
| 177 | + vse32.v v6, (ccp); add ccp, ccp, cstride; |
| 178 | + vse32.v v7, (ccp); add ccp, ccp, cstride; |
| 179 | + vse32.v v8, (ccp); add ccp, ccp, cstride; |
| 180 | + vse32.v v9, (ccp); add ccp, ccp, cstride; |
| 181 | + vse32.v v10, (ccp); add ccp, ccp, cstride; |
| 182 | + vse32.v v11, (ccp); add ccp, ccp, cstride; |
| 183 | + vse32.v v12, (ccp); add ccp, ccp, cstride; |
| 184 | + vse32.v v13, (ccp); add ccp, ccp, cstride; |
| 185 | + vse32.v v14, (ccp); add ccp, ccp, cstride; |
| 186 | + vse32.v v15, (ccp) |
| 187 | + |
| 188 | + # Following tail instructions should be scheduled earlier in free slots during C block save. |
| 189 | + # Leaving here for clarity. |
| 190 | + |
| 191 | + # Bump pointers for loop across blocks in one row |
| 192 | + slli t6, nvl, 2 |
| 193 | + add cnp, cnp, t6 # Move C block pointer over |
| 194 | + add bnp, bnp, t6 # Move B block pointer over |
| 195 | + sub nt, nt, nvl # Decrement element count in n dimension |
| 196 | + bnez nt, c_col_loop # Any more to do? |
| 197 | + |
| 198 | + # Move to next set of rows |
| 199 | + addi m, m, -16 # Did 16 rows above |
| 200 | + slli t6, astride, 4 # Multiply astride by 16 |
| 201 | + add ap, ap, t6 # Move A matrix pointer down 16 rows |
| 202 | + slli t6, cstride, 4 # Multiply cstride by 16 |
| 203 | + add cp, cp, t6 # Move C matrix pointer down 16 rows |
| 204 | + |
| 205 | + slti t6, m, 16 |
| 206 | + beqz t6, c_row_loop |
| 207 | + |
| 208 | + # Handle end of matrix with fewer than 16 rows. |
| 209 | + # Can use smaller versions of above decreasing in powers-of-2 depending on code-size concerns. |
| 210 | +end_rows: |
| 211 | + # Not done. |
| 212 | + |
| 213 | +exit: |
| 214 | + ld s0, OFFSET(sp) |
| 215 | + ld s1, OFFSET(sp) |
| 216 | + ld s2, OFFSET(sp) |
| 217 | + addi sp, sp, FRAMESIZE |
| 218 | + ret |
0 commit comments