• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1#! /usr/bin/env perl
2# Copyright 2012-2016 The OpenSSL Project Authors. All Rights Reserved.
3#
4# Licensed under the OpenSSL license (the "License").  You may not use
5# this file except in compliance with the License.  You can obtain a copy
6# in the file LICENSE in the source distribution or at
7# https://www.openssl.org/source/license.html
8
9
10# ====================================================================
11# Written by Andy Polyakov <appro@openssl.org> for the OpenSSL
12# project. The module is, however, dual licensed under OpenSSL and
13# CRYPTOGAMS licenses depending on where you obtain it. For further
14# details see http://www.openssl.org/~appro/cryptogams/.
15#
16# Specific modes and adaptation for Linux kernel by Ard Biesheuvel
17# of Linaro. Permission to use under GPL terms is granted.
18# ====================================================================
19
20# Bit-sliced AES for ARM NEON
21#
22# February 2012.
23#
24# This implementation is direct adaptation of bsaes-x86_64 module for
25# ARM NEON. Except that this module is endian-neutral [in sense that
26# it can be compiled for either endianness] by courtesy of vld1.8's
27# neutrality. Initial version doesn't implement interface to OpenSSL,
28# only low-level primitives and unsupported entry points, just enough
29# to collect performance results, which for Cortex-A8 core are:
30#
31# encrypt	19.5 cycles per byte processed with 128-bit key
32# decrypt	22.1 cycles per byte processed with 128-bit key
33# key conv.	440  cycles per 128-bit key/0.18 of 8x block
34#
35# Snapdragon S4 encrypts byte in 17.6 cycles and decrypts in 19.7,
36# which is [much] worse than anticipated (for further details see
37# http://www.openssl.org/~appro/Snapdragon-S4.html).
38#
39# Cortex-A15 manages in 14.2/16.1 cycles [when integer-only code
40# manages in 20.0 cycles].
41#
42# When comparing to x86_64 results keep in mind that NEON unit is
43# [mostly] single-issue and thus can't [fully] benefit from
44# instruction-level parallelism. And when comparing to aes-armv4
45# results keep in mind key schedule conversion overhead (see
46# bsaes-x86_64.pl for further details)...
47#
48#						<appro@openssl.org>
49
50# April-August 2013
51# Add CBC, CTR and XTS subroutines and adapt for kernel use; courtesy of Ard.
52
53$flavour = shift;
54if ($flavour=~/\w[\w\-]*\.\w+$/) { $output=$flavour; undef $flavour; }
55else { while (($output=shift) && ($output!~/\w[\w\-]*\.\w+$/)) {} }
56
57if ($flavour && $flavour ne "void") {
58    $0 =~ m/(.*[\/\\])[^\/\\]+$/; $dir=$1;
59    ( $xlate="${dir}arm-xlate.pl" and -f $xlate ) or
60    ( $xlate="${dir}../../../perlasm/arm-xlate.pl" and -f $xlate) or
61    die "can't locate arm-xlate.pl";
62
63    open OUT,"| \"$^X\" $xlate $flavour $output";
64    *STDOUT=*OUT;
65} else {
66    open OUT,">$output";
67    *STDOUT=*OUT;
68}
69
70my ($inp,$out,$len,$key)=("r0","r1","r2","r3");
71my @XMM=map("q$_",(0..15));
72
73{
74my ($key,$rounds,$const)=("r4","r5","r6");
75
76sub Dlo()   { shift=~m|q([1]?[0-9])|?"d".($1*2):"";     }
77sub Dhi()   { shift=~m|q([1]?[0-9])|?"d".($1*2+1):"";   }
78
79sub Sbox {
80# input in  lsb > [b0, b1, b2, b3, b4, b5, b6, b7] < msb
81# output in lsb > [b0, b1, b4, b6, b3, b7, b2, b5] < msb
82my @b=@_[0..7];
83my @t=@_[8..11];
84my @s=@_[12..15];
85	&InBasisChange	(@b);
86	&Inv_GF256	(@b[6,5,0,3,7,1,4,2],@t,@s);
87	&OutBasisChange	(@b[7,1,4,2,6,5,0,3]);
88}
89
90sub InBasisChange {
91# input in  lsb > [b0, b1, b2, b3, b4, b5, b6, b7] < msb
92# output in lsb > [b6, b5, b0, b3, b7, b1, b4, b2] < msb
93my @b=@_[0..7];
94$code.=<<___;
95	veor	@b[2], @b[2], @b[1]
96	veor	@b[5], @b[5], @b[6]
97	veor	@b[3], @b[3], @b[0]
98	veor	@b[6], @b[6], @b[2]
99	veor	@b[5], @b[5], @b[0]
100
101	veor	@b[6], @b[6], @b[3]
102	veor	@b[3], @b[3], @b[7]
103	veor	@b[7], @b[7], @b[5]
104	veor	@b[3], @b[3], @b[4]
105	veor	@b[4], @b[4], @b[5]
106
107	veor	@b[2], @b[2], @b[7]
108	veor	@b[3], @b[3], @b[1]
109	veor	@b[1], @b[1], @b[5]
110___
111}
112
113sub OutBasisChange {
114# input in  lsb > [b0, b1, b2, b3, b4, b5, b6, b7] < msb
115# output in lsb > [b6, b1, b2, b4, b7, b0, b3, b5] < msb
116my @b=@_[0..7];
117$code.=<<___;
118	veor	@b[0], @b[0], @b[6]
119	veor	@b[1], @b[1], @b[4]
120	veor	@b[4], @b[4], @b[6]
121	veor	@b[2], @b[2], @b[0]
122	veor	@b[6], @b[6], @b[1]
123
124	veor	@b[1], @b[1], @b[5]
125	veor	@b[5], @b[5], @b[3]
126	veor	@b[3], @b[3], @b[7]
127	veor	@b[7], @b[7], @b[5]
128	veor	@b[2], @b[2], @b[5]
129
130	veor	@b[4], @b[4], @b[7]
131___
132}
133
134sub InvSbox {
135# input in lsb 	> [b0, b1, b2, b3, b4, b5, b6, b7] < msb
136# output in lsb	> [b0, b1, b6, b4, b2, b7, b3, b5] < msb
137my @b=@_[0..7];
138my @t=@_[8..11];
139my @s=@_[12..15];
140	&InvInBasisChange	(@b);
141	&Inv_GF256		(@b[5,1,2,6,3,7,0,4],@t,@s);
142	&InvOutBasisChange	(@b[3,7,0,4,5,1,2,6]);
143}
144
145sub InvInBasisChange {		# OutBasisChange in reverse (with twist)
146my @b=@_[5,1,2,6,3,7,0,4];
147$code.=<<___
148	 veor	@b[1], @b[1], @b[7]
149	veor	@b[4], @b[4], @b[7]
150
151	veor	@b[7], @b[7], @b[5]
152	 veor	@b[1], @b[1], @b[3]
153	veor	@b[2], @b[2], @b[5]
154	veor	@b[3], @b[3], @b[7]
155
156	veor	@b[6], @b[6], @b[1]
157	veor	@b[2], @b[2], @b[0]
158	 veor	@b[5], @b[5], @b[3]
159	veor	@b[4], @b[4], @b[6]
160	veor	@b[0], @b[0], @b[6]
161	veor	@b[1], @b[1], @b[4]
162___
163}
164
165sub InvOutBasisChange {		# InBasisChange in reverse
166my @b=@_[2,5,7,3,6,1,0,4];
167$code.=<<___;
168	veor	@b[1], @b[1], @b[5]
169	veor	@b[2], @b[2], @b[7]
170
171	veor	@b[3], @b[3], @b[1]
172	veor	@b[4], @b[4], @b[5]
173	veor	@b[7], @b[7], @b[5]
174	veor	@b[3], @b[3], @b[4]
175	 veor 	@b[5], @b[5], @b[0]
176	veor	@b[3], @b[3], @b[7]
177	 veor	@b[6], @b[6], @b[2]
178	 veor	@b[2], @b[2], @b[1]
179	veor	@b[6], @b[6], @b[3]
180
181	veor	@b[3], @b[3], @b[0]
182	veor	@b[5], @b[5], @b[6]
183___
184}
185
186sub Mul_GF4 {
187#;*************************************************************
188#;* Mul_GF4: Input x0-x1,y0-y1 Output x0-x1 Temp t0 (8) *
189#;*************************************************************
190my ($x0,$x1,$y0,$y1,$t0,$t1)=@_;
191$code.=<<___;
192	veor 	$t0, $y0, $y1
193	vand	$t0, $t0, $x0
194	veor	$x0, $x0, $x1
195	vand	$t1, $x1, $y0
196	vand	$x0, $x0, $y1
197	veor	$x1, $t1, $t0
198	veor	$x0, $x0, $t1
199___
200}
201
202sub Mul_GF4_N {				# not used, see next subroutine
203# multiply and scale by N
204my ($x0,$x1,$y0,$y1,$t0)=@_;
205$code.=<<___;
206	veor	$t0, $y0, $y1
207	vand	$t0, $t0, $x0
208	veor	$x0, $x0, $x1
209	vand	$x1, $x1, $y0
210	vand	$x0, $x0, $y1
211	veor	$x1, $x1, $x0
212	veor	$x0, $x0, $t0
213___
214}
215
216sub Mul_GF4_N_GF4 {
217# interleaved Mul_GF4_N and Mul_GF4
218my ($x0,$x1,$y0,$y1,$t0,
219    $x2,$x3,$y2,$y3,$t1)=@_;
220$code.=<<___;
221	veor	$t0, $y0, $y1
222	 veor 	$t1, $y2, $y3
223	vand	$t0, $t0, $x0
224	 vand	$t1, $t1, $x2
225	veor	$x0, $x0, $x1
226	 veor	$x2, $x2, $x3
227	vand	$x1, $x1, $y0
228	 vand	$x3, $x3, $y2
229	vand	$x0, $x0, $y1
230	 vand	$x2, $x2, $y3
231	veor	$x1, $x1, $x0
232	 veor	$x2, $x2, $x3
233	veor	$x0, $x0, $t0
234	 veor	$x3, $x3, $t1
235___
236}
237sub Mul_GF16_2 {
238my @x=@_[0..7];
239my @y=@_[8..11];
240my @t=@_[12..15];
241$code.=<<___;
242	veor	@t[0], @x[0], @x[2]
243	veor	@t[1], @x[1], @x[3]
244___
245	&Mul_GF4  	(@x[0], @x[1], @y[0], @y[1], @t[2..3]);
246$code.=<<___;
247	veor	@y[0], @y[0], @y[2]
248	veor	@y[1], @y[1], @y[3]
249___
250	Mul_GF4_N_GF4	(@t[0], @t[1], @y[0], @y[1], @t[3],
251			 @x[2], @x[3], @y[2], @y[3], @t[2]);
252$code.=<<___;
253	veor	@x[0], @x[0], @t[0]
254	veor	@x[2], @x[2], @t[0]
255	veor	@x[1], @x[1], @t[1]
256	veor	@x[3], @x[3], @t[1]
257
258	veor	@t[0], @x[4], @x[6]
259	veor	@t[1], @x[5], @x[7]
260___
261	&Mul_GF4_N_GF4	(@t[0], @t[1], @y[0], @y[1], @t[3],
262			 @x[6], @x[7], @y[2], @y[3], @t[2]);
263$code.=<<___;
264	veor	@y[0], @y[0], @y[2]
265	veor	@y[1], @y[1], @y[3]
266___
267	&Mul_GF4  	(@x[4], @x[5], @y[0], @y[1], @t[2..3]);
268$code.=<<___;
269	veor	@x[4], @x[4], @t[0]
270	veor	@x[6], @x[6], @t[0]
271	veor	@x[5], @x[5], @t[1]
272	veor	@x[7], @x[7], @t[1]
273___
274}
275sub Inv_GF256 {
276#;********************************************************************
277#;* Inv_GF256: Input x0-x7 Output x0-x7 Temp t0-t3,s0-s3 (144)       *
278#;********************************************************************
279my @x=@_[0..7];
280my @t=@_[8..11];
281my @s=@_[12..15];
282# direct optimizations from hardware
283$code.=<<___;
284	veor	@t[3], @x[4], @x[6]
285	veor	@t[2], @x[5], @x[7]
286	veor	@t[1], @x[1], @x[3]
287	veor	@s[1], @x[7], @x[6]
288	 vmov	@t[0], @t[2]
289	veor	@s[0], @x[0], @x[2]
290
291	vorr	@t[2], @t[2], @t[1]
292	veor	@s[3], @t[3], @t[0]
293	vand	@s[2], @t[3], @s[0]
294	vorr	@t[3], @t[3], @s[0]
295	veor	@s[0], @s[0], @t[1]
296	vand	@t[0], @t[0], @t[1]
297	veor	@t[1], @x[3], @x[2]
298	vand	@s[3], @s[3], @s[0]
299	vand	@s[1], @s[1], @t[1]
300	veor	@t[1], @x[4], @x[5]
301	veor	@s[0], @x[1], @x[0]
302	veor	@t[3], @t[3], @s[1]
303	veor	@t[2], @t[2], @s[1]
304	vand	@s[1], @t[1], @s[0]
305	vorr	@t[1], @t[1], @s[0]
306	veor	@t[3], @t[3], @s[3]
307	veor	@t[0], @t[0], @s[1]
308	veor	@t[2], @t[2], @s[2]
309	veor	@t[1], @t[1], @s[3]
310	veor	@t[0], @t[0], @s[2]
311	vand	@s[0], @x[7], @x[3]
312	veor	@t[1], @t[1], @s[2]
313	vand	@s[1], @x[6], @x[2]
314	vand	@s[2], @x[5], @x[1]
315	vorr	@s[3], @x[4], @x[0]
316	veor	@t[3], @t[3], @s[0]
317	veor	@t[1], @t[1], @s[2]
318	veor	@t[0], @t[0], @s[3]
319	veor	@t[2], @t[2], @s[1]
320
321	@ Inv_GF16 \t0, \t1, \t2, \t3, \s0, \s1, \s2, \s3
322
323	@ new smaller inversion
324
325	vand	@s[2], @t[3], @t[1]
326	vmov	@s[0], @t[0]
327
328	veor	@s[1], @t[2], @s[2]
329	veor	@s[3], @t[0], @s[2]
330	veor	@s[2], @t[0], @s[2]	@ @s[2]=@s[3]
331
332	vbsl	@s[1], @t[1], @t[0]
333	vbsl	@s[3], @t[3], @t[2]
334	veor	@t[3], @t[3], @t[2]
335
336	vbsl	@s[0], @s[1], @s[2]
337	vbsl	@t[0], @s[2], @s[1]
338
339	vand	@s[2], @s[0], @s[3]
340	veor	@t[1], @t[1], @t[0]
341
342	veor	@s[2], @s[2], @t[3]
343___
344# output in s3, s2, s1, t1
345
346# Mul_GF16_2 \x0, \x1, \x2, \x3, \x4, \x5, \x6, \x7, \t2, \t3, \t0, \t1, \s0, \s1, \s2, \s3
347
348# Mul_GF16_2 \x0, \x1, \x2, \x3, \x4, \x5, \x6, \x7, \s3, \s2, \s1, \t1, \s0, \t0, \t2, \t3
349	&Mul_GF16_2(@x,@s[3,2,1],@t[1],@s[0],@t[0,2,3]);
350
351### output msb > [x3,x2,x1,x0,x7,x6,x5,x4] < lsb
352}
353
354# AES linear components
355
356sub ShiftRows {
357my @x=@_[0..7];
358my @t=@_[8..11];
359my $mask=pop;
360$code.=<<___;
361	vldmia	$key!, {@t[0]-@t[3]}
362	veor	@t[0], @t[0], @x[0]
363	veor	@t[1], @t[1], @x[1]
364	vtbl.8	`&Dlo(@x[0])`, {@t[0]}, `&Dlo($mask)`
365	vtbl.8	`&Dhi(@x[0])`, {@t[0]}, `&Dhi($mask)`
366	vldmia	$key!, {@t[0]}
367	veor	@t[2], @t[2], @x[2]
368	vtbl.8	`&Dlo(@x[1])`, {@t[1]}, `&Dlo($mask)`
369	vtbl.8	`&Dhi(@x[1])`, {@t[1]}, `&Dhi($mask)`
370	vldmia	$key!, {@t[1]}
371	veor	@t[3], @t[3], @x[3]
372	vtbl.8	`&Dlo(@x[2])`, {@t[2]}, `&Dlo($mask)`
373	vtbl.8	`&Dhi(@x[2])`, {@t[2]}, `&Dhi($mask)`
374	vldmia	$key!, {@t[2]}
375	vtbl.8	`&Dlo(@x[3])`, {@t[3]}, `&Dlo($mask)`
376	vtbl.8	`&Dhi(@x[3])`, {@t[3]}, `&Dhi($mask)`
377	vldmia	$key!, {@t[3]}
378	veor	@t[0], @t[0], @x[4]
379	veor	@t[1], @t[1], @x[5]
380	vtbl.8	`&Dlo(@x[4])`, {@t[0]}, `&Dlo($mask)`
381	vtbl.8	`&Dhi(@x[4])`, {@t[0]}, `&Dhi($mask)`
382	veor	@t[2], @t[2], @x[6]
383	vtbl.8	`&Dlo(@x[5])`, {@t[1]}, `&Dlo($mask)`
384	vtbl.8	`&Dhi(@x[5])`, {@t[1]}, `&Dhi($mask)`
385	veor	@t[3], @t[3], @x[7]
386	vtbl.8	`&Dlo(@x[6])`, {@t[2]}, `&Dlo($mask)`
387	vtbl.8	`&Dhi(@x[6])`, {@t[2]}, `&Dhi($mask)`
388	vtbl.8	`&Dlo(@x[7])`, {@t[3]}, `&Dlo($mask)`
389	vtbl.8	`&Dhi(@x[7])`, {@t[3]}, `&Dhi($mask)`
390___
391}
392
393sub MixColumns {
394# modified to emit output in order suitable for feeding back to aesenc[last]
395my @x=@_[0..7];
396my @t=@_[8..15];
397my $inv=@_[16];	# optional
398$code.=<<___;
399	vext.8	@t[0], @x[0], @x[0], #12	@ x0 <<< 32
400	vext.8	@t[1], @x[1], @x[1], #12
401	 veor	@x[0], @x[0], @t[0]		@ x0 ^ (x0 <<< 32)
402	vext.8	@t[2], @x[2], @x[2], #12
403	 veor	@x[1], @x[1], @t[1]
404	vext.8	@t[3], @x[3], @x[3], #12
405	 veor	@x[2], @x[2], @t[2]
406	vext.8	@t[4], @x[4], @x[4], #12
407	 veor	@x[3], @x[3], @t[3]
408	vext.8	@t[5], @x[5], @x[5], #12
409	 veor	@x[4], @x[4], @t[4]
410	vext.8	@t[6], @x[6], @x[6], #12
411	 veor	@x[5], @x[5], @t[5]
412	vext.8	@t[7], @x[7], @x[7], #12
413	 veor	@x[6], @x[6], @t[6]
414
415	veor	@t[1], @t[1], @x[0]
416	 veor	@x[7], @x[7], @t[7]
417	 vext.8	@x[0], @x[0], @x[0], #8		@ (x0 ^ (x0 <<< 32)) <<< 64)
418	veor	@t[2], @t[2], @x[1]
419	veor	@t[0], @t[0], @x[7]
420	veor	@t[1], @t[1], @x[7]
421	 vext.8	@x[1], @x[1], @x[1], #8
422	veor	@t[5], @t[5], @x[4]
423	 veor	@x[0], @x[0], @t[0]
424	veor	@t[6], @t[6], @x[5]
425	 veor	@x[1], @x[1], @t[1]
426	 vext.8	@t[0], @x[4], @x[4], #8
427	veor	@t[4], @t[4], @x[3]
428	 vext.8	@t[1], @x[5], @x[5], #8
429	veor	@t[7], @t[7], @x[6]
430	 vext.8	@x[4], @x[3], @x[3], #8
431	veor	@t[3], @t[3], @x[2]
432	 vext.8	@x[5], @x[7], @x[7], #8
433	veor	@t[4], @t[4], @x[7]
434	 vext.8	@x[3], @x[6], @x[6], #8
435	veor	@t[3], @t[3], @x[7]
436	 vext.8	@x[6], @x[2], @x[2], #8
437	veor	@x[7], @t[1], @t[5]
438___
439$code.=<<___ if (!$inv);
440	veor	@x[2], @t[0], @t[4]
441	veor	@x[4], @x[4], @t[3]
442	veor	@x[5], @x[5], @t[7]
443	veor	@x[3], @x[3], @t[6]
444	 @ vmov	@x[2], @t[0]
445	veor	@x[6], @x[6], @t[2]
446	 @ vmov	@x[7], @t[1]
447___
448$code.=<<___ if ($inv);
449	veor	@t[3], @t[3], @x[4]
450	veor	@x[5], @x[5], @t[7]
451	veor	@x[2], @x[3], @t[6]
452	veor	@x[3], @t[0], @t[4]
453	veor	@x[4], @x[6], @t[2]
454	vmov	@x[6], @t[3]
455	 @ vmov	@x[7], @t[1]
456___
457}
458
459sub InvMixColumns_orig {
460my @x=@_[0..7];
461my @t=@_[8..15];
462
463$code.=<<___;
464	@ multiplication by 0x0e
465	vext.8	@t[7], @x[7], @x[7], #12
466	vmov	@t[2], @x[2]
467	veor	@x[2], @x[2], @x[5]		@ 2 5
468	veor	@x[7], @x[7], @x[5]		@ 7 5
469	vext.8	@t[0], @x[0], @x[0], #12
470	vmov	@t[5], @x[5]
471	veor	@x[5], @x[5], @x[0]		@ 5 0		[1]
472	veor	@x[0], @x[0], @x[1]		@ 0 1
473	vext.8	@t[1], @x[1], @x[1], #12
474	veor	@x[1], @x[1], @x[2]		@ 1 25
475	veor	@x[0], @x[0], @x[6]		@ 01 6		[2]
476	vext.8	@t[3], @x[3], @x[3], #12
477	veor	@x[1], @x[1], @x[3]		@ 125 3		[4]
478	veor	@x[2], @x[2], @x[0]		@ 25 016	[3]
479	veor	@x[3], @x[3], @x[7]		@ 3 75
480	veor	@x[7], @x[7], @x[6]		@ 75 6		[0]
481	vext.8	@t[6], @x[6], @x[6], #12
482	vmov	@t[4], @x[4]
483	veor	@x[6], @x[6], @x[4]		@ 6 4
484	veor	@x[4], @x[4], @x[3]		@ 4 375		[6]
485	veor	@x[3], @x[3], @x[7]		@ 375 756=36
486	veor	@x[6], @x[6], @t[5]		@ 64 5		[7]
487	veor	@x[3], @x[3], @t[2]		@ 36 2
488	vext.8	@t[5], @t[5], @t[5], #12
489	veor	@x[3], @x[3], @t[4]		@ 362 4		[5]
490___
491					my @y = @x[7,5,0,2,1,3,4,6];
492$code.=<<___;
493	@ multiplication by 0x0b
494	veor	@y[1], @y[1], @y[0]
495	veor	@y[0], @y[0], @t[0]
496	vext.8	@t[2], @t[2], @t[2], #12
497	veor	@y[1], @y[1], @t[1]
498	veor	@y[0], @y[0], @t[5]
499	vext.8	@t[4], @t[4], @t[4], #12
500	veor	@y[1], @y[1], @t[6]
501	veor	@y[0], @y[0], @t[7]
502	veor	@t[7], @t[7], @t[6]		@ clobber t[7]
503
504	veor	@y[3], @y[3], @t[0]
505	 veor	@y[1], @y[1], @y[0]
506	vext.8	@t[0], @t[0], @t[0], #12
507	veor	@y[2], @y[2], @t[1]
508	veor	@y[4], @y[4], @t[1]
509	vext.8	@t[1], @t[1], @t[1], #12
510	veor	@y[2], @y[2], @t[2]
511	veor	@y[3], @y[3], @t[2]
512	veor	@y[5], @y[5], @t[2]
513	veor	@y[2], @y[2], @t[7]
514	vext.8	@t[2], @t[2], @t[2], #12
515	veor	@y[3], @y[3], @t[3]
516	veor	@y[6], @y[6], @t[3]
517	veor	@y[4], @y[4], @t[3]
518	veor	@y[7], @y[7], @t[4]
519	vext.8	@t[3], @t[3], @t[3], #12
520	veor	@y[5], @y[5], @t[4]
521	veor	@y[7], @y[7], @t[7]
522	veor	@t[7], @t[7], @t[5]		@ clobber t[7] even more
523	veor	@y[3], @y[3], @t[5]
524	veor	@y[4], @y[4], @t[4]
525
526	veor	@y[5], @y[5], @t[7]
527	vext.8	@t[4], @t[4], @t[4], #12
528	veor	@y[6], @y[6], @t[7]
529	veor	@y[4], @y[4], @t[7]
530
531	veor	@t[7], @t[7], @t[5]
532	vext.8	@t[5], @t[5], @t[5], #12
533
534	@ multiplication by 0x0d
535	veor	@y[4], @y[4], @y[7]
536	 veor	@t[7], @t[7], @t[6]		@ restore t[7]
537	veor	@y[7], @y[7], @t[4]
538	vext.8	@t[6], @t[6], @t[6], #12
539	veor	@y[2], @y[2], @t[0]
540	veor	@y[7], @y[7], @t[5]
541	vext.8	@t[7], @t[7], @t[7], #12
542	veor	@y[2], @y[2], @t[2]
543
544	veor	@y[3], @y[3], @y[1]
545	veor	@y[1], @y[1], @t[1]
546	veor	@y[0], @y[0], @t[0]
547	veor	@y[3], @y[3], @t[0]
548	veor	@y[1], @y[1], @t[5]
549	veor	@y[0], @y[0], @t[5]
550	vext.8	@t[0], @t[0], @t[0], #12
551	veor	@y[1], @y[1], @t[7]
552	veor	@y[0], @y[0], @t[6]
553	veor	@y[3], @y[3], @y[1]
554	veor	@y[4], @y[4], @t[1]
555	vext.8	@t[1], @t[1], @t[1], #12
556
557	veor	@y[7], @y[7], @t[7]
558	veor	@y[4], @y[4], @t[2]
559	veor	@y[5], @y[5], @t[2]
560	veor	@y[2], @y[2], @t[6]
561	veor	@t[6], @t[6], @t[3]		@ clobber t[6]
562	vext.8	@t[2], @t[2], @t[2], #12
563	veor	@y[4], @y[4], @y[7]
564	veor	@y[3], @y[3], @t[6]
565
566	veor	@y[6], @y[6], @t[6]
567	veor	@y[5], @y[5], @t[5]
568	vext.8	@t[5], @t[5], @t[5], #12
569	veor	@y[6], @y[6], @t[4]
570	vext.8	@t[4], @t[4], @t[4], #12
571	veor	@y[5], @y[5], @t[6]
572	veor	@y[6], @y[6], @t[7]
573	vext.8	@t[7], @t[7], @t[7], #12
574	veor	@t[6], @t[6], @t[3]		@ restore t[6]
575	vext.8	@t[3], @t[3], @t[3], #12
576
577	@ multiplication by 0x09
578	veor	@y[4], @y[4], @y[1]
579	veor	@t[1], @t[1], @y[1]		@ t[1]=y[1]
580	veor	@t[0], @t[0], @t[5]		@ clobber t[0]
581	vext.8	@t[6], @t[6], @t[6], #12
582	veor	@t[1], @t[1], @t[5]
583	veor	@y[3], @y[3], @t[0]
584	veor	@t[0], @t[0], @y[0]		@ t[0]=y[0]
585	veor	@t[1], @t[1], @t[6]
586	veor	@t[6], @t[6], @t[7]		@ clobber t[6]
587	veor	@y[4], @y[4], @t[1]
588	veor	@y[7], @y[7], @t[4]
589	veor	@y[6], @y[6], @t[3]
590	veor	@y[5], @y[5], @t[2]
591	veor	@t[4], @t[4], @y[4]		@ t[4]=y[4]
592	veor	@t[3], @t[3], @y[3]		@ t[3]=y[3]
593	veor	@t[5], @t[5], @y[5]		@ t[5]=y[5]
594	veor	@t[2], @t[2], @y[2]		@ t[2]=y[2]
595	veor	@t[3], @t[3], @t[7]
596	veor	@XMM[5], @t[5], @t[6]
597	veor	@XMM[6], @t[6], @y[6]		@ t[6]=y[6]
598	veor	@XMM[2], @t[2], @t[6]
599	veor	@XMM[7], @t[7], @y[7]		@ t[7]=y[7]
600
601	vmov	@XMM[0], @t[0]
602	vmov	@XMM[1], @t[1]
603	@ vmov	@XMM[2], @t[2]
604	vmov	@XMM[3], @t[3]
605	vmov	@XMM[4], @t[4]
606	@ vmov	@XMM[5], @t[5]
607	@ vmov	@XMM[6], @t[6]
608	@ vmov	@XMM[7], @t[7]
609___
610}
611
612sub InvMixColumns {
613my @x=@_[0..7];
614my @t=@_[8..15];
615
616# Thanks to Jussi Kivilinna for providing pointer to
617#
618# | 0e 0b 0d 09 |   | 02 03 01 01 |   | 05 00 04 00 |
619# | 09 0e 0b 0d | = | 01 02 03 01 | x | 00 05 00 04 |
620# | 0d 09 0e 0b |   | 01 01 02 03 |   | 04 00 05 00 |
621# | 0b 0d 09 0e |   | 03 01 01 02 |   | 00 04 00 05 |
622
623$code.=<<___;
624	@ multiplication by 0x05-0x00-0x04-0x00
625	vext.8	@t[0], @x[0], @x[0], #8
626	vext.8	@t[6], @x[6], @x[6], #8
627	vext.8	@t[7], @x[7], @x[7], #8
628	veor	@t[0], @t[0], @x[0]
629	vext.8	@t[1], @x[1], @x[1], #8
630	veor	@t[6], @t[6], @x[6]
631	vext.8	@t[2], @x[2], @x[2], #8
632	veor	@t[7], @t[7], @x[7]
633	vext.8	@t[3], @x[3], @x[3], #8
634	veor	@t[1], @t[1], @x[1]
635	vext.8	@t[4], @x[4], @x[4], #8
636	veor	@t[2], @t[2], @x[2]
637	vext.8	@t[5], @x[5], @x[5], #8
638	veor	@t[3], @t[3], @x[3]
639	veor	@t[4], @t[4], @x[4]
640	veor	@t[5], @t[5], @x[5]
641
642	 veor	@x[0], @x[0], @t[6]
643	 veor	@x[1], @x[1], @t[6]
644	 veor	@x[2], @x[2], @t[0]
645	 veor	@x[4], @x[4], @t[2]
646	 veor	@x[3], @x[3], @t[1]
647	 veor	@x[1], @x[1], @t[7]
648	 veor	@x[2], @x[2], @t[7]
649	 veor	@x[4], @x[4], @t[6]
650	 veor	@x[5], @x[5], @t[3]
651	 veor	@x[3], @x[3], @t[6]
652	 veor	@x[6], @x[6], @t[4]
653	 veor	@x[4], @x[4], @t[7]
654	 veor	@x[5], @x[5], @t[7]
655	 veor	@x[7], @x[7], @t[5]
656___
657	&MixColumns	(@x,@t,1);	# flipped 2<->3 and 4<->6
658}
659
660sub swapmove {
661my ($a,$b,$n,$mask,$t)=@_;
662$code.=<<___;
663	vshr.u64	$t, $b, #$n
664	veor		$t, $t, $a
665	vand		$t, $t, $mask
666	veor		$a, $a, $t
667	vshl.u64	$t, $t, #$n
668	veor		$b, $b, $t
669___
670}
671sub swapmove2x {
672my ($a0,$b0,$a1,$b1,$n,$mask,$t0,$t1)=@_;
673$code.=<<___;
674	vshr.u64	$t0, $b0, #$n
675	 vshr.u64	$t1, $b1, #$n
676	veor		$t0, $t0, $a0
677	 veor		$t1, $t1, $a1
678	vand		$t0, $t0, $mask
679	 vand		$t1, $t1, $mask
680	veor		$a0, $a0, $t0
681	vshl.u64	$t0, $t0, #$n
682	 veor		$a1, $a1, $t1
683	 vshl.u64	$t1, $t1, #$n
684	veor		$b0, $b0, $t0
685	 veor		$b1, $b1, $t1
686___
687}
688
689sub bitslice {
690my @x=reverse(@_[0..7]);
691my ($t0,$t1,$t2,$t3)=@_[8..11];
692$code.=<<___;
693	vmov.i8	$t0,#0x55			@ compose .LBS0
694	vmov.i8	$t1,#0x33			@ compose .LBS1
695___
696	&swapmove2x(@x[0,1,2,3],1,$t0,$t2,$t3);
697	&swapmove2x(@x[4,5,6,7],1,$t0,$t2,$t3);
698$code.=<<___;
699	vmov.i8	$t0,#0x0f			@ compose .LBS2
700___
701	&swapmove2x(@x[0,2,1,3],2,$t1,$t2,$t3);
702	&swapmove2x(@x[4,6,5,7],2,$t1,$t2,$t3);
703
704	&swapmove2x(@x[0,4,1,5],4,$t0,$t2,$t3);
705	&swapmove2x(@x[2,6,3,7],4,$t0,$t2,$t3);
706}
707
708$code.=<<___;
709#ifndef __KERNEL__
710# include <GFp/arm_arch.h>
711
712# define VFP_ABI_PUSH	vstmdb	sp!,{d8-d15}
713# define VFP_ABI_POP	vldmia	sp!,{d8-d15}
714# define VFP_ABI_FRAME	0x40
715#else
716# define VFP_ABI_PUSH
717# define VFP_ABI_POP
718# define VFP_ABI_FRAME	0
719# define BSAES_ASM_EXTENDED_KEY
720# define __ARM_ARCH__ __LINUX_ARM_ARCH__
721# define __ARM_MAX_ARCH__ 7
722#endif
723
724#ifdef __thumb__
725# define adrl adr
726#endif
727
728#if __ARM_MAX_ARCH__>=7
729.arch	armv7-a
730.fpu	neon
731
732.text
733.syntax	unified 	@ ARMv7-capable assembler is expected to handle this
734#if defined(__thumb2__) && !defined(__APPLE__)
735.thumb
736#else
737.code   32
738# undef __thumb2__
739#endif
740
741.type	_bsaes_const,%object
742.align	6
743_bsaes_const:
744.LM0ISR:	@ InvShiftRows constants
745	.quad	0x0a0e0206070b0f03, 0x0004080c0d010509
746.LISR:
747	.quad	0x0504070602010003, 0x0f0e0d0c080b0a09
748.LISRM0:
749	.quad	0x01040b0e0205080f, 0x0306090c00070a0d
750.LM0SR:		@ ShiftRows constants
751	.quad	0x0a0e02060f03070b, 0x0004080c05090d01
752.LSR:
753	.quad	0x0504070600030201, 0x0f0e0d0c0a09080b
754.LSRM0:
755	.quad	0x0304090e00050a0f, 0x01060b0c0207080d
756.LM0:
757	.quad	0x02060a0e03070b0f, 0x0004080c0105090d
758.LREVM0SR:
759	.quad	0x090d01050c000408, 0x03070b0f060a0e02
760.asciz	"Bit-sliced AES for NEON, CRYPTOGAMS by <appro\@openssl.org>"
761.align	6
762.size	_bsaes_const,.-_bsaes_const
763
764.type	_bsaes_encrypt8,%function
765.align	4
766_bsaes_encrypt8:
767	adr	$const,.
768	vldmia	$key!, {@XMM[9]}		@ round 0 key
769#if defined(__thumb2__) || defined(__APPLE__)
770	adr	$const,.LM0SR
771#else
772	sub	$const,$const,#_bsaes_encrypt8-.LM0SR
773#endif
774
775	vldmia	$const!, {@XMM[8]}		@ .LM0SR
776_bsaes_encrypt8_alt:
777	veor	@XMM[10], @XMM[0], @XMM[9]	@ xor with round0 key
778	veor	@XMM[11], @XMM[1], @XMM[9]
779	 vtbl.8	`&Dlo(@XMM[0])`, {@XMM[10]}, `&Dlo(@XMM[8])`
780	 vtbl.8	`&Dhi(@XMM[0])`, {@XMM[10]}, `&Dhi(@XMM[8])`
781	veor	@XMM[12], @XMM[2], @XMM[9]
782	 vtbl.8	`&Dlo(@XMM[1])`, {@XMM[11]}, `&Dlo(@XMM[8])`
783	 vtbl.8	`&Dhi(@XMM[1])`, {@XMM[11]}, `&Dhi(@XMM[8])`
784	veor	@XMM[13], @XMM[3], @XMM[9]
785	 vtbl.8	`&Dlo(@XMM[2])`, {@XMM[12]}, `&Dlo(@XMM[8])`
786	 vtbl.8	`&Dhi(@XMM[2])`, {@XMM[12]}, `&Dhi(@XMM[8])`
787	veor	@XMM[14], @XMM[4], @XMM[9]
788	 vtbl.8	`&Dlo(@XMM[3])`, {@XMM[13]}, `&Dlo(@XMM[8])`
789	 vtbl.8	`&Dhi(@XMM[3])`, {@XMM[13]}, `&Dhi(@XMM[8])`
790	veor	@XMM[15], @XMM[5], @XMM[9]
791	 vtbl.8	`&Dlo(@XMM[4])`, {@XMM[14]}, `&Dlo(@XMM[8])`
792	 vtbl.8	`&Dhi(@XMM[4])`, {@XMM[14]}, `&Dhi(@XMM[8])`
793	veor	@XMM[10], @XMM[6], @XMM[9]
794	 vtbl.8	`&Dlo(@XMM[5])`, {@XMM[15]}, `&Dlo(@XMM[8])`
795	 vtbl.8	`&Dhi(@XMM[5])`, {@XMM[15]}, `&Dhi(@XMM[8])`
796	veor	@XMM[11], @XMM[7], @XMM[9]
797	 vtbl.8	`&Dlo(@XMM[6])`, {@XMM[10]}, `&Dlo(@XMM[8])`
798	 vtbl.8	`&Dhi(@XMM[6])`, {@XMM[10]}, `&Dhi(@XMM[8])`
799	 vtbl.8	`&Dlo(@XMM[7])`, {@XMM[11]}, `&Dlo(@XMM[8])`
800	 vtbl.8	`&Dhi(@XMM[7])`, {@XMM[11]}, `&Dhi(@XMM[8])`
801_bsaes_encrypt8_bitslice:
802___
803	&bitslice	(@XMM[0..7, 8..11]);
804$code.=<<___;
805	sub	$rounds,$rounds,#1
806	b	.Lenc_sbox
807.align	4
808.Lenc_loop:
809___
810	&ShiftRows	(@XMM[0..7, 8..12]);
811$code.=".Lenc_sbox:\n";
812	&Sbox		(@XMM[0..7, 8..15]);
813$code.=<<___;
814	subs	$rounds,$rounds,#1
815	bcc	.Lenc_done
816___
817	&MixColumns	(@XMM[0,1,4,6,3,7,2,5, 8..15]);
818$code.=<<___;
819	vldmia	$const, {@XMM[12]}		@ .LSR
820	ite	eq				@ Thumb2 thing, samity check in ARM
821	addeq	$const,$const,#0x10
822	bne	.Lenc_loop
823	vldmia	$const, {@XMM[12]}		@ .LSRM0
824	b	.Lenc_loop
825.align	4
826.Lenc_done:
827___
828	# output in lsb > [t0, t1, t4, t6, t3, t7, t2, t5] < msb
829	&bitslice	(@XMM[0,1,4,6,3,7,2,5, 8..11]);
830$code.=<<___;
831	vldmia	$key, {@XMM[8]}			@ last round key
832	veor	@XMM[4], @XMM[4], @XMM[8]
833	veor	@XMM[6], @XMM[6], @XMM[8]
834	veor	@XMM[3], @XMM[3], @XMM[8]
835	veor	@XMM[7], @XMM[7], @XMM[8]
836	veor	@XMM[2], @XMM[2], @XMM[8]
837	veor	@XMM[5], @XMM[5], @XMM[8]
838	veor	@XMM[0], @XMM[0], @XMM[8]
839	veor	@XMM[1], @XMM[1], @XMM[8]
840	bx	lr
841.size	_bsaes_encrypt8,.-_bsaes_encrypt8
842___
843}
844{
845my ($out,$inp,$rounds,$const)=("r12","r4","r5","r6");
846
847sub bitslice_key {
848my @x=reverse(@_[0..7]);
849my ($bs0,$bs1,$bs2,$t2,$t3)=@_[8..12];
850
851	&swapmove	(@x[0,1],1,$bs0,$t2,$t3);
852$code.=<<___;
853	@ &swapmove(@x[2,3],1,$t0,$t2,$t3);
854	vmov	@x[2], @x[0]
855	vmov	@x[3], @x[1]
856___
857	#&swapmove2x(@x[4,5,6,7],1,$t0,$t2,$t3);
858
859	&swapmove2x	(@x[0,2,1,3],2,$bs1,$t2,$t3);
860$code.=<<___;
861	@ &swapmove2x(@x[4,6,5,7],2,$t1,$t2,$t3);
862	vmov	@x[4], @x[0]
863	vmov	@x[6], @x[2]
864	vmov	@x[5], @x[1]
865	vmov	@x[7], @x[3]
866___
867	&swapmove2x	(@x[0,4,1,5],4,$bs2,$t2,$t3);
868	&swapmove2x	(@x[2,6,3,7],4,$bs2,$t2,$t3);
869}
870
871$code.=<<___;
872.type	_bsaes_key_convert,%function
873.align	4
874_bsaes_key_convert:
875	adr	$const,.
876	vld1.8	{@XMM[7]},  [$inp]!		@ load round 0 key
877#if defined(__thumb2__) || defined(__APPLE__)
878	adr	$const,.LM0
879#else
880	sub	$const,$const,#_bsaes_key_convert-.LM0
881#endif
882	vld1.8	{@XMM[15]}, [$inp]!		@ load round 1 key
883
884	vmov.i8	@XMM[8],  #0x01			@ bit masks
885	vmov.i8	@XMM[9],  #0x02
886	vmov.i8	@XMM[10], #0x04
887	vmov.i8	@XMM[11], #0x08
888	vmov.i8	@XMM[12], #0x10
889	vmov.i8	@XMM[13], #0x20
890	vldmia	$const, {@XMM[14]}		@ .LM0
891
892#ifdef __ARMEL__
893	vrev32.8	@XMM[7],  @XMM[7]
894	vrev32.8	@XMM[15], @XMM[15]
895#endif
896	sub	$rounds,$rounds,#1
897	vstmia	$out!, {@XMM[7]}		@ save round 0 key
898	b	.Lkey_loop
899
900.align	4
901.Lkey_loop:
902	vtbl.8	`&Dlo(@XMM[7])`,{@XMM[15]},`&Dlo(@XMM[14])`
903	vtbl.8	`&Dhi(@XMM[7])`,{@XMM[15]},`&Dhi(@XMM[14])`
904	vmov.i8	@XMM[6],  #0x40
905	vmov.i8	@XMM[15], #0x80
906
907	vtst.8	@XMM[0], @XMM[7], @XMM[8]
908	vtst.8	@XMM[1], @XMM[7], @XMM[9]
909	vtst.8	@XMM[2], @XMM[7], @XMM[10]
910	vtst.8	@XMM[3], @XMM[7], @XMM[11]
911	vtst.8	@XMM[4], @XMM[7], @XMM[12]
912	vtst.8	@XMM[5], @XMM[7], @XMM[13]
913	vtst.8	@XMM[6], @XMM[7], @XMM[6]
914	vtst.8	@XMM[7], @XMM[7], @XMM[15]
915	vld1.8	{@XMM[15]}, [$inp]!		@ load next round key
916	vmvn	@XMM[0], @XMM[0]		@ "pnot"
917	vmvn	@XMM[1], @XMM[1]
918	vmvn	@XMM[5], @XMM[5]
919	vmvn	@XMM[6], @XMM[6]
920#ifdef __ARMEL__
921	vrev32.8	@XMM[15], @XMM[15]
922#endif
923	subs	$rounds,$rounds,#1
924	vstmia	$out!,{@XMM[0]-@XMM[7]}		@ write bit-sliced round key
925	bne	.Lkey_loop
926
927	vmov.i8	@XMM[7],#0x63			@ compose .L63
928	@ don't save last round key
929	bx	lr
930.size	_bsaes_key_convert,.-_bsaes_key_convert
931___
932}
933
934{
935my ($inp,$out,$len,$key, $ctr,$fp,$rounds)=(map("r$_",(0..3,8..10)));
936my $const = "r6";	# shared with _bsaes_encrypt8_alt
937my $keysched = "sp";
938
939$code.=<<___;
940.global	GFp_bsaes_ctr32_encrypt_blocks
941.type	GFp_bsaes_ctr32_encrypt_blocks,%function
942.align	5
943GFp_bsaes_ctr32_encrypt_blocks:
944	@ In OpenSSL, short inputs fall back to aes_nohw_* here. We patch this
945	@ out to retain a constant-time implementation.
946	mov	ip, sp
947	stmdb	sp!, {r4-r10, lr}
948	VFP_ABI_PUSH
949	ldr	$ctr, [ip]			@ ctr is 1st arg on the stack
950	sub	sp, sp, #0x10			@ scratch space to carry over the ctr
951	mov	$fp, sp				@ save sp
952
953	ldr	$rounds, [$key, #240]		@ get # of rounds
954#ifndef	BSAES_ASM_EXTENDED_KEY
955	@ allocate the key schedule on the stack
956	sub	r12, sp, $rounds, lsl#7		@ 128 bytes per inner round key
957	add	r12, #`128-32`			@ size of bit-sliced key schedule
958
959	@ populate the key schedule
960	mov	r4, $key			@ pass key
961	mov	r5, $rounds			@ pass # of rounds
962	mov	sp, r12				@ sp is $keysched
963	bl	_bsaes_key_convert
964	veor	@XMM[7],@XMM[7],@XMM[15]	@ fix up last round key
965	vstmia	r12, {@XMM[7]}			@ save last round key
966
967	vld1.8	{@XMM[0]}, [$ctr]		@ load counter
968#ifdef	__APPLE__
969	mov	$ctr, #:lower16:(.LREVM0SR-.LM0)
970	add	$ctr, $const, $ctr
971#else
972	add	$ctr, $const, #.LREVM0SR-.LM0	@ borrow $ctr
973#endif
974	vldmia	$keysched, {@XMM[4]}		@ load round0 key
975#else
976	ldr	r12, [$key, #244]
977	eors	r12, #1
978	beq	0f
979
980	@ populate the key schedule
981	str	r12, [$key, #244]
982	mov	r4, $key			@ pass key
983	mov	r5, $rounds			@ pass # of rounds
984	add	r12, $key, #248			@ pass key schedule
985	bl	_bsaes_key_convert
986	veor	@XMM[7],@XMM[7],@XMM[15]	@ fix up last round key
987	vstmia	r12, {@XMM[7]}			@ save last round key
988
989.align	2
9900:	add	r12, $key, #248
991	vld1.8	{@XMM[0]}, [$ctr]		@ load counter
992	adrl	$ctr, .LREVM0SR			@ borrow $ctr
993	vldmia	r12, {@XMM[4]}			@ load round0 key
994	sub	sp, #0x10			@ place for adjusted round0 key
995#endif
996
997	vmov.i32	@XMM[8],#1		@ compose 1<<96
998	veor		@XMM[9],@XMM[9],@XMM[9]
999	vrev32.8	@XMM[0],@XMM[0]
1000	vext.8		@XMM[8],@XMM[9],@XMM[8],#4
1001	vrev32.8	@XMM[4],@XMM[4]
1002	vadd.u32	@XMM[9],@XMM[8],@XMM[8]	@ compose 2<<96
1003	vstmia	$keysched, {@XMM[4]}		@ save adjusted round0 key
1004	b	.Lctr_enc_loop
1005
1006.align	4
1007.Lctr_enc_loop:
1008	vadd.u32	@XMM[10], @XMM[8], @XMM[9]	@ compose 3<<96
1009	vadd.u32	@XMM[1], @XMM[0], @XMM[8]	@ +1
1010	vadd.u32	@XMM[2], @XMM[0], @XMM[9]	@ +2
1011	vadd.u32	@XMM[3], @XMM[0], @XMM[10]	@ +3
1012	vadd.u32	@XMM[4], @XMM[1], @XMM[10]
1013	vadd.u32	@XMM[5], @XMM[2], @XMM[10]
1014	vadd.u32	@XMM[6], @XMM[3], @XMM[10]
1015	vadd.u32	@XMM[7], @XMM[4], @XMM[10]
1016	vadd.u32	@XMM[10], @XMM[5], @XMM[10]	@ next counter
1017
1018	@ Borrow prologue from _bsaes_encrypt8 to use the opportunity
1019	@ to flip byte order in 32-bit counter
1020
1021	vldmia		$keysched, {@XMM[9]}		@ load round0 key
1022#ifndef	BSAES_ASM_EXTENDED_KEY
1023	add		r4, $keysched, #0x10		@ pass next round key
1024#else
1025	add		r4, $key, #`248+16`
1026#endif
1027	vldmia		$ctr, {@XMM[8]}			@ .LREVM0SR
1028	mov		r5, $rounds			@ pass rounds
1029	vstmia		$fp, {@XMM[10]}			@ save next counter
1030#ifdef	__APPLE__
1031	mov		$const, #:lower16:(.LREVM0SR-.LSR)
1032	sub		$const, $ctr, $const
1033#else
1034	sub		$const, $ctr, #.LREVM0SR-.LSR	@ pass constants
1035#endif
1036
1037	bl		_bsaes_encrypt8_alt
1038
1039	subs		$len, $len, #8
1040	blo		.Lctr_enc_loop_done
1041
1042	vld1.8		{@XMM[8]-@XMM[9]}, [$inp]!	@ load input
1043	vld1.8		{@XMM[10]-@XMM[11]}, [$inp]!
1044	veor		@XMM[0], @XMM[8]
1045	veor		@XMM[1], @XMM[9]
1046	vld1.8		{@XMM[12]-@XMM[13]}, [$inp]!
1047	veor		@XMM[4], @XMM[10]
1048	veor		@XMM[6], @XMM[11]
1049	vld1.8		{@XMM[14]-@XMM[15]}, [$inp]!
1050	veor		@XMM[3], @XMM[12]
1051	vst1.8		{@XMM[0]-@XMM[1]}, [$out]!	@ write output
1052	veor		@XMM[7], @XMM[13]
1053	veor		@XMM[2], @XMM[14]
1054	vst1.8		{@XMM[4]}, [$out]!
1055	veor		@XMM[5], @XMM[15]
1056	vst1.8		{@XMM[6]}, [$out]!
1057	vmov.i32	@XMM[8], #1			@ compose 1<<96
1058	vst1.8		{@XMM[3]}, [$out]!
1059	veor		@XMM[9], @XMM[9], @XMM[9]
1060	vst1.8		{@XMM[7]}, [$out]!
1061	vext.8		@XMM[8], @XMM[9], @XMM[8], #4
1062	vst1.8		{@XMM[2]}, [$out]!
1063	vadd.u32	@XMM[9],@XMM[8],@XMM[8]		@ compose 2<<96
1064	vst1.8		{@XMM[5]}, [$out]!
1065	vldmia		$fp, {@XMM[0]}			@ load counter
1066
1067	bne		.Lctr_enc_loop
1068	b		.Lctr_enc_done
1069
1070.align	4
1071.Lctr_enc_loop_done:
1072	add		$len, $len, #8
1073	vld1.8		{@XMM[8]}, [$inp]!	@ load input
1074	veor		@XMM[0], @XMM[8]
1075	vst1.8		{@XMM[0]}, [$out]!	@ write output
1076	cmp		$len, #2
1077	blo		.Lctr_enc_done
1078	vld1.8		{@XMM[9]}, [$inp]!
1079	veor		@XMM[1], @XMM[9]
1080	vst1.8		{@XMM[1]}, [$out]!
1081	beq		.Lctr_enc_done
1082	vld1.8		{@XMM[10]}, [$inp]!
1083	veor		@XMM[4], @XMM[10]
1084	vst1.8		{@XMM[4]}, [$out]!
1085	cmp		$len, #4
1086	blo		.Lctr_enc_done
1087	vld1.8		{@XMM[11]}, [$inp]!
1088	veor		@XMM[6], @XMM[11]
1089	vst1.8		{@XMM[6]}, [$out]!
1090	beq		.Lctr_enc_done
1091	vld1.8		{@XMM[12]}, [$inp]!
1092	veor		@XMM[3], @XMM[12]
1093	vst1.8		{@XMM[3]}, [$out]!
1094	cmp		$len, #6
1095	blo		.Lctr_enc_done
1096	vld1.8		{@XMM[13]}, [$inp]!
1097	veor		@XMM[7], @XMM[13]
1098	vst1.8		{@XMM[7]}, [$out]!
1099	beq		.Lctr_enc_done
1100	vld1.8		{@XMM[14]}, [$inp]
1101	veor		@XMM[2], @XMM[14]
1102	vst1.8		{@XMM[2]}, [$out]!
1103
1104.Lctr_enc_done:
1105	vmov.i32	q0, #0
1106	vmov.i32	q1, #0
1107#ifndef	BSAES_ASM_EXTENDED_KEY
1108.Lctr_enc_bzero:			@ wipe key schedule [if any]
1109	vstmia		$keysched!, {q0-q1}
1110	cmp		$keysched, $fp
1111	bne		.Lctr_enc_bzero
1112#else
1113	vstmia		$keysched, {q0-q1}
1114#endif
1115
1116	mov	sp, $fp
1117	add	sp, #0x10		@ add sp,$fp,#0x10 is no good for thumb
1118	VFP_ABI_POP
1119	ldmia	sp!, {r4-r10, pc}	@ return
1120
1121	@ OpenSSL contains aes_nohw_* fallback code here. We patch this
1122	@ out to retain a constant-time implementation.
1123.size	GFp_bsaes_ctr32_encrypt_blocks,.-GFp_bsaes_ctr32_encrypt_blocks
1124___
1125}
1126$code.=<<___;
1127#endif
1128___
1129
1130$code =~ s/\`([^\`]*)\`/eval($1)/gem;
1131
1132open SELF,$0;
1133while(<SELF>) {
1134	next if (/^#!/);
1135        last if (!s/^#/@/ and !/^$/);
1136        print;
1137}
1138close SELF;
1139
1140print $code;
1141
1142close STDOUT or die "error closing STDOUT";
1143