ntt_n256_l8_s32_mont.s

///
/// Copyright (c) 2021 Arm Limited
/// Copyright (c) 2022 Hanno Becker
/// Copyright (c) 2023 Amin Abdulrahman, Matthias Kannwischer
/// SPDX-License-Identifier: MIT
///
/// Permission is hereby granted, free of charge, to any person obtaining a copy
/// of this software and associated documentation files (the "Software"), to deal
/// in the Software without restriction, including without limitation the rights
/// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
/// copies of the Software, and to permit persons to whom the Software is
/// furnished to do so, subject to the following conditions:
///
/// The above copyright notice and this permission notice shall be included in all
/// copies or substantial portions of the Software.
///
/// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
/// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
/// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
/// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
/// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
/// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
/// SOFTWARE.
///



///
/// This assembly code has been auto-generated.
/// Don't modify it directly.
///

.data
roots:
.word   29095681 // zeta^128 * 2^31 = 28678040^128 * 2^31 = 17702291 * 2^31
.word 3280343807 // zeta^128 * f(q^(-1) mod 2^32) * 2^31 = 28678040^128 * 375649793 * 2^31
.word   14476917 // zeta^ 64 * 2^31 = 28678040^ 64 * 2^31 = 3260327 * 2^31
.word 2356128651 // zeta^ 64 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 64 * 375649793 * 2^31
.word   43317805 // zeta^192 * 2^31 = 28678040^192 * 2^31 = 14579576 * 2^31
.word  933021651 // zeta^192 * f(q^(-1) mod 2^32) * 2^31 = 28678040^192 * 375649793 * 2^31
.word   18598075 // zeta^ 32 * 2^31 = 28678040^ 32 * 2^31 = 6733847 * 2^31
.word 2578416965 // zeta^ 32 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 32 * 375649793 * 2^31
.word   39999747 // zeta^ 16 * 2^31 = 28678040^ 16 * 2^31 = 20428075 * 2^31
.word 3454780669 // zeta^ 16 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 16 * 375649793 * 2^31
.word   45317587 // zeta^144 * 2^31 = 28678040^144 * 2^31 = 14626653 * 2^31
.word 3083517997 // zeta^144 * f(q^(-1) mod 2^32) * 2^31 = 28678040^144 * 375649793 * 2^31
.word    4885007 // zeta^160 * 2^31 = 28678040^160 * 2^31 = 12909577 * 2^31
.word 2973633521 // zeta^160 * f(q^(-1) mod 2^32) * 2^31 = 28678040^160 * 375649793 * 2^31
.word   48811299 // zeta^ 80 * 2^31 = 28678040^ 80 * 2^31 = 29737761 * 2^31
.word 4050555101 // zeta^ 80 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 80 * 375649793 * 2^31
.word   54571669 // zeta^208 * 2^31 = 28678040^208 * 2^31 = 30285189 * 2^31
.word 4085587819 // zeta^208 * f(q^(-1) mod 2^32) * 2^31 = 28678040^208 * 375649793 * 2^31
.word   64683161 // zeta^ 96 * 2^31 = 28678040^ 96 * 2^31 = 14745691 * 2^31
.word 3091135847 // zeta^ 96 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 96 * 375649793 * 2^31
.word   59281651 // zeta^ 48 * 2^31 = 28678040^ 48 * 2^31 = 21289485 * 2^31
.word 3509906701 // zeta^ 48 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 48 * 375649793 * 2^31
.word   40500013 // zeta^176 * 2^31 = 28678040^176 * 2^31 = 9914896 * 2^31
.word  634504915 // zeta^176 * f(q^(-1) mod 2^32) * 2^31 = 28678040^176 * 375649793 * 2^31
.word   34427601 // zeta^224 * 2^31 = 28678040^224 * 2^31 = 13512548 * 2^31
.word  864737071 // zeta^224 * f(q^(-1) mod 2^32) * 2^31 = 28678040^224 * 375649793 * 2^31
.word   25917637 // zeta^112 * 2^31 = 28678040^112 * 2^31 = 22603682 * 2^31
.word 1446525243 // zeta^112 * f(q^(-1) mod 2^32) * 2^31 = 28678040^112 * 375649793 * 2^31
.word    8356523 // zeta^240 * 2^31 = 28678040^240 * 2^31 = 16204162 * 2^31
.word 1036987221 // zeta^240 * f(q^(-1) mod 2^32) * 2^31 = 28678040^240 * 375649793 * 2^31
.word   31719253 // zeta^  8 * 2^31 = 28678040^  8 * 2^31 = 23825509 * 2^31
.word 3672199851 // zeta^  8 * f(q^(-1) mod 2^32) * 2^31 = 28678040^  8 * 375649793 * 2^31
.word    5075563 // zeta^  4 * 2^31 = 28678040^  4 * 2^31 = 9010590 * 2^31
.word  576633749 // zeta^  4 * f(q^(-1) mod 2^32) * 2^31 = 28678040^  4 * 375649793 * 2^31
.word   43115375 // zeta^132 * 2^31 = 28678040^132 * 2^31 = 20699126 * 2^31
.word 1324642961 // zeta^132 * f(q^(-1) mod 2^32) * 2^31 = 28678040^132 * 375649793 * 2^31
.word   54842419 // zeta^136 * 2^31 = 28678040^136 * 2^31 = 27028662 * 2^31
.word 1729702349 // zeta^136 * f(q^(-1) mod 2^32) * 2^31 = 28678040^136 * 375649793 * 2^31
.word   35131011 // zeta^ 68 * 2^31 = 28678040^ 68 * 2^31 = 341080 * 2^31
.word   21827453 // zeta^ 68 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 68 * 375649793 * 2^31
.word   44664611 // zeta^196 * 2^31 = 28678040^196 * 2^31 = 21220783 * 2^31
.word 3505510109 // zeta^196 * f(q^(-1) mod 2^32) * 2^31 = 28678040^196 * 375649793 * 2^31
.word    1316163 // zeta^ 72 * 2^31 = 28678040^ 72 * 2^31 = 14833295 * 2^31
.word 3096742077 // zeta^ 72 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 72 * 375649793 * 2^31
.word   65968403 // zeta^ 36 * 2^31 = 28678040^ 36 * 2^31 = 25331745 * 2^31
.word 3768591597 // zeta^ 36 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 36 * 375649793 * 2^31
.word   53949037 // zeta^164 * 2^31 = 28678040^164 * 2^31 = 5289426 * 2^31
.word  338497427 // zeta^164 * f(q^(-1) mod 2^32) * 2^31 = 28678040^164 * 375649793 * 2^31
.word   10391631 // zeta^200 * 2^31 = 28678040^200 * 2^31 = 2138810 * 2^31
.word  136873393 // zeta^200 * f(q^(-1) mod 2^32) * 2^31 = 28678040^200 * 375649793 * 2^31
.word   52363231 // zeta^100 * 2^31 = 28678040^100 * 2^31 = 5705868 * 2^31
.word  365147681 // zeta^100 * f(q^(-1) mod 2^32) * 2^31 = 28678040^100 * 375649793 * 2^31
.word   39928117 // zeta^228 * 2^31 = 28678040^228 * 2^31 = 17686665 * 2^31
.word 3279343819 // zeta^228 * f(q^(-1) mod 2^32) * 2^31 = 28678040^228 * 375649793 * 2^31
.word   54335767 // zeta^ 40 * 2^31 = 28678040^ 40 * 2^31 = 6490403 * 2^31
.word 2562837737 // zeta^ 40 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 40 * 375649793 * 2^31
.word   54457727 // zeta^ 20 * 2^31 = 28678040^ 20 * 2^31 = 9106105 * 2^31
.word 2730229889 // zeta^ 20 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 20 * 375649793 * 2^31
.word   27596809 // zeta^148 * 2^31 = 28678040^148 * 2^31 = 18817700 * 2^31
.word 1204240887 // zeta^148 * f(q^(-1) mod 2^32) * 2^31 = 28678040^148 * 375649793 * 2^31
.word   46002083 // zeta^168 * 2^31 = 28678040^168 * 2^31 = 19648405 * 2^31
.word 3404885597 // zeta^168 * f(q^(-1) mod 2^32) * 2^31 = 28678040^168 * 375649793 * 2^31
.word   14847715 // zeta^ 84 * 2^31 = 28678040^ 84 * 2^31 = 1579445 * 2^31
.word 2248560413 // zeta^ 84 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 84 * 375649793 * 2^31
.word    1129279 // zeta^212 * 2^31 = 28678040^212 * 2^31 = 7769916 * 2^31
.word  497236673 // zeta^212 * f(q^(-1) mod 2^32) * 2^31 = 28678040^212 * 375649793 * 2^31
.word   35733845 // zeta^104 * 2^31 = 28678040^104 * 2^31 = 31254932 * 2^31
.word 2000162987 // zeta^104 * f(q^(-1) mod 2^32) * 2^31 = 28678040^104 * 375649793 * 2^31
.word   54563587 // zeta^ 52 * 2^31 = 28678040^ 52 * 2^31 = 21843119 * 2^31
.word 3545336573 // zeta^ 52 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 52 * 375649793 * 2^31
.word   35404977 // zeta^180 * 2^31 = 28678040^180 * 2^31 = 11828796 * 2^31
.word  756985167 // zeta^180 * f(q^(-1) mod 2^32) * 2^31 = 28678040^180 * 375649793 * 2^31
.word   61099389 // zeta^232 * 2^31 = 28678040^232 * 2^31 = 26362414 * 2^31
.word 1687065731 // zeta^232 * f(q^(-1) mod 2^32) * 2^31 = 28678040^232 * 375649793 * 2^31
.word   52947923 // zeta^116 * 2^31 = 28678040^116 * 2^31 = 19828530 * 2^31
.word 1268929069 // zeta^116 * f(q^(-1) mod 2^32) * 2^31 = 28678040^116 * 375649793 * 2^31
.word   41822583 // zeta^244 * 2^31 = 28678040^244 * 2^31 = 33201112 * 2^31
.word 2124709001 // zeta^244 * f(q^(-1) mod 2^32) * 2^31 = 28678040^244 * 375649793 * 2^31
.word   26241327 // zeta^ 24 * 2^31 = 28678040^ 24 * 2^31 = 572895 * 2^31
.word 2184146129 // zeta^ 24 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 24 * 375649793 * 2^31
.word   12770159 // zeta^ 12 * 2^31 = 28678040^ 12 * 2^31 = 23713020 * 2^31
.word 1517517457 // zeta^ 12 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 12 * 375649793 * 2^31
.word   24980679 // zeta^140 * 2^31 = 28678040^140 * 2^31 = 19537976 * 2^31
.word 1250335033 // zeta^140 * f(q^(-1) mod 2^32) * 2^31 = 28678040^140 * 375649793 * 2^31
.word    5033605 // zeta^152 * 2^31 = 28678040^152 * 2^31 = 26691971 * 2^31
.word 3855639419 // zeta^152 * f(q^(-1) mod 2^32) * 2^31 = 28678040^152 * 375649793 * 2^31
.word   61827033 // zeta^ 76 * 2^31 = 28678040^ 76 * 2^31 = 8285889 * 2^31
.word 2677740071 // zeta^ 76 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 76 * 375649793 * 2^31
.word   11221523 // zeta^204 * 2^31 = 28678040^204 * 2^31 = 24690028 * 2^31
.word 1580041197 // zeta^204 * f(q^(-1) mod 2^32) * 2^31 = 28678040^204 * 375649793 * 2^31
.word    8316793 // zeta^ 88 * 2^31 = 28678040^ 88 * 2^31 = 9249292 * 2^31
.word  591909511 // zeta^ 88 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 88 * 375649793 * 2^31
.word   19091691 // zeta^ 44 * 2^31 = 28678040^ 44 * 2^31 = 4778209 * 2^31
.word 2453265685 // zeta^ 44 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 44 * 375649793 * 2^31
.word   32210035 // zeta^172 * 2^31 = 28678040^172 * 2^31 = 13113327 * 2^31
.word 2986672525 // zeta^172 * f(q^(-1) mod 2^32) * 2^31 = 28678040^172 * 375649793 * 2^31
.word   16634213 // zeta^216 * 2^31 = 28678040^216 * 2^31 = 29292862 * 2^31
.word 1874600091 // zeta^216 * f(q^(-1) mod 2^32) * 2^31 = 28678040^216 * 375649793 * 2^31
.word   20871313 // zeta^108 * 2^31 = 28678040^108 * 2^31 = 25384023 * 2^31
.word 3771937135 // zeta^108 * f(q^(-1) mod 2^32) * 2^31 = 28678040^108 * 375649793 * 2^31
.word   46581651 // zeta^236 * 2^31 = 28678040^236 * 2^31 = 10905370 * 2^31
.word  697890413 // zeta^236 * f(q^(-1) mod 2^32) * 2^31 = 28678040^236 * 375649793 * 2^31
.word   63329695 // zeta^ 56 * 2^31 = 28678040^ 56 * 2^31 = 8247799 * 2^31
.word 2675302497 // zeta^ 56 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 56 * 375649793 * 2^31
.word   51221435 // zeta^ 28 * 2^31 = 28678040^ 28 * 2^31 = 16167867 * 2^31
.word 3182148165 // zeta^ 28 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 28 * 375649793 * 2^31
.word   18467171 // zeta^156 * 2^31 = 28678040^156 * 2^31 = 22046437 * 2^31
.word 3558347933 // zeta^156 * f(q^(-1) mod 2^32) * 2^31 = 28678040^156 * 375649793 * 2^31
.word    9983051 // zeta^184 * 2^31 = 28678040^184 * 2^31 = 5086187 * 2^31
.word 2472974773 // zeta^184 * f(q^(-1) mod 2^32) * 2^31 = 28678040^184 * 375649793 * 2^31
.word   37083207 // zeta^ 92 * 2^31 = 28678040^ 92 * 2^31 = 656361 * 2^31
.word 2189487545 // zeta^ 92 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 92 * 375649793 * 2^31
.word   52674527 // zeta^220 * 2^31 = 28678040^220 * 2^31 = 18153794 * 2^31
.word 1161754145 // zeta^220 * f(q^(-1) mod 2^32) * 2^31 = 28678040^220 * 375649793 * 2^31
.word    7721125 // zeta^120 * 2^31 = 28678040^120 * 2^31 = 28113639 * 2^31
.word 3946619227 // zeta^120 * f(q^(-1) mod 2^32) * 2^31 = 28678040^120 * 375649793 * 2^31
.word    8896309 // zeta^ 60 * 2^31 = 28678040^ 60 * 2^31 = 3732072 * 2^31
.word  238834379 // zeta^ 60 * f(q^(-1) mod 2^32) * 2^31 = 28678040^ 60 * 375649793 * 2^31
.word    2061353 // zeta^188 * 2^31 = 28678040^188 * 2^31 = 22126384 * 2^31
.word 1415980503 // zeta^188 * f(q^(-1) mod 2^32) * 2^31 = 28678040^188 * 375649793 * 2^31
.word    9383201 // zeta^248 * 2^31 = 28678040^248 * 2^31 = 8471290 * 2^31
.word  542121183 // zeta^248 * f(q^(-1) mod 2^32) * 2^31 = 28678040^248 * 375649793 * 2^31
.word   23761465 // zeta^124 * 2^31 = 28678040^124 * 2^31 = 9445744 * 2^31
.word  604481479 // zeta^124 * f(q^(-1) mod 2^32) * 2^31 = 28678040^124 * 375649793 * 2^31
.word   24512363 // zeta^252 * 2^31 = 28678040^252 * 2^31 = 794839 * 2^31
.word 2198349461 // zeta^252 * f(q^(-1) mod 2^32) * 2^31 = 28678040^252 * 375649793 * 2^31
.word   13704133 // zeta^  2 * 2^31 = 28678040^  2 * 2^31
.word   41177999 // zeta^130 * 2^31 = 28678040^130 * 2^31
.word   26703739 // zeta^ 66 * 2^31 = 28678040^ 66 * 2^31
.word   65289035 // zeta^194 * 2^31 = 28678040^194 * 2^31
.word 1666225723 // zeta^  2 * (q^(-1) mod 2^32) * 2^31 = 28678040^  2 * 375649793 * 2^31
.word 2599633521 // zeta^130 * (q^(-1) mod 2^32) * 2^31 = 28678040^130 * 375649793 * 2^31
.word 2869384837 // zeta^ 66 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 66 * 375649793 * 2^31
.word 1260434101 // zeta^194 * (q^(-1) mod 2^32) * 2^31 = 28678040^194 * 375649793 * 2^31
.word   50326315 // zeta^  1 * 2^31 = 28678040^  1 * 2^31
.word   37746191 // zeta^ 65 * 2^31 = 28678040^ 65 * 2^31
.word   49080301 // zeta^ 33 * 2^31 = 28678040^ 33 * 2^31
.word   34232193 // zeta^ 97 * 2^31 = 28678040^ 97 * 2^31
.word 1835254485 // zeta^  1 * (q^(-1) mod 2^32) * 2^31 = 28678040^  1 * 375649793 * 2^31
.word  360751089 // zeta^ 65 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 65 * 375649793 * 2^31
.word 1200511507 // zeta^ 33 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 33 * 375649793 * 2^31
.word  553431679 // zeta^ 97 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 97 * 375649793 * 2^31
.word   22955837 // zeta^129 * 2^31 = 28678040^129 * 2^31
.word   31411079 // zeta^193 * 2^31 = 28678040^193 * 2^31
.word     492607 // zeta^161 * 2^31 = 28678040^161 * 2^31
.word   22217509 // zeta^225 * 2^31 = 28678040^225 * 2^31
.word 2610305731 // zeta^129 * (q^(-1) mod 2^32) * 2^31 = 28678040^129 * 375649793 * 2^31
.word 2623011449 // zeta^193 * (q^(-1) mod 2^32) * 2^31 = 28678040^193 * 375649793 * 2^31
.word  948367809 // zeta^161 * (q^(-1) mod 2^32) * 2^31 = 28678040^161 * 375649793 * 2^31
.word 3562565339 // zeta^225 * (q^(-1) mod 2^32) * 2^31 = 28678040^225 * 375649793 * 2^31
.word    5481609 // zeta^ 34 * 2^31 = 28678040^ 34 * 2^31
.word   12552175 // zeta^162 * 2^31 = 28678040^162 * 2^31
.word   54494203 // zeta^ 98 * 2^31 = 28678040^ 98 * 2^31
.word   32704019 // zeta^226 * 2^31 = 28678040^226 * 2^31
.word  949335415 // zeta^ 34 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 34 * 375649793 * 2^31
.word 3610496529 // zeta^162 * (q^(-1) mod 2^32) * 2^31 = 28678040^162 * 375649793 * 2^31
.word 1474054661 // zeta^ 98 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 98 * 375649793 * 2^31
.word 2061350893 // zeta^226 * (q^(-1) mod 2^32) * 2^31 = 28678040^226 * 375649793 * 2^31
.word   48767307 // zeta^ 17 * 2^31 = 28678040^ 17 * 2^31
.word   39600285 // zeta^ 81 * 2^31 = 28678040^ 81 * 2^31
.word   31654617 // zeta^ 49 * 2^31 = 28678040^ 49 * 2^31
.word    4736231 // zeta^113 * 2^31 = 28678040^113 * 2^31
.word 2602093749 // zeta^ 17 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 17 * 375649793 * 2^31
.word 3705004387 // zeta^ 81 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 81 * 375649793 * 2^31
.word  427128615 // zeta^ 49 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 49 * 375649793 * 2^31
.word  237814041 // zeta^113 * (q^(-1) mod 2^32) * 2^31 = 28678040^113 * 375649793 * 2^31
.word   18965555 // zeta^145 * 2^31 = 28678040^145 * 2^31
.word   50771049 // zeta^209 * 2^31 = 28678040^209 * 2^31
.word    8794671 // zeta^177 * 2^31 = 28678040^177 * 2^31
.word   59508707 // zeta^241 * 2^31 = 28678040^241 * 2^31
.word 3336346061 // zeta^145 * (q^(-1) mod 2^32) * 2^31 = 28678040^145 * 375649793 * 2^31
.word 3238094231 // zeta^209 * (q^(-1) mod 2^32) * 2^31 = 28678040^209 * 375649793 * 2^31
.word 2568201169 // zeta^177 * (q^(-1) mod 2^32) * 2^31 = 28678040^177 * 375649793 * 2^31
.word 3726038557 // zeta^241 * (q^(-1) mod 2^32) * 2^31 = 28678040^241 * 375649793 * 2^31
.word   43973433 // zeta^ 18 * 2^31 = 28678040^ 18 * 2^31
.word   14453865 // zeta^146 * 2^31 = 28678040^146 * 2^31
.word   14937153 // zeta^ 82 * 2^31 = 28678040^ 82 * 2^31
.word   39701997 // zeta^210 * 2^31 = 28678040^210 * 2^31
.word  720191175 // zeta^ 18 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 18 * 375649793 * 2^31
.word 3181088151 // zeta^146 * (q^(-1) mod 2^32) * 2^31 = 28678040^146 * 375649793 * 2^31
.word  116563391 // zeta^ 82 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 82 * 375649793 * 2^31
.word 3642323987 // zeta^210 * (q^(-1) mod 2^32) * 2^31 = 28678040^210 * 375649793 * 2^31
.word   53455571 // zeta^  9 * 2^31 = 28678040^  9 * 2^31
.word   35877127 // zeta^ 73 * 2^31 = 28678040^ 73 * 2^31
.word     681755 // zeta^ 41 * 2^31 = 28678040^ 41 * 2^31
.word   63245537 // zeta^105 * 2^31 = 28678040^105 * 2^31
.word 4245721901 // zeta^  9 * (q^(-1) mod 2^32) * 2^31 = 28678040^  9 * 375649793 * 2^31
.word 2676675833 // zeta^ 73 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 73 * 375649793 * 2^31
.word 3480266469 // zeta^ 41 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 41 * 375649793 * 2^31
.word 1356315935 // zeta^105 * (q^(-1) mod 2^32) * 2^31 = 28678040^105 * 375649793 * 2^31
.word   11718751 // zeta^137 * 2^31 = 28678040^137 * 2^31
.word   41885553 // zeta^201 * 2^31 = 28678040^201 * 2^31
.word   54210213 // zeta^169 * 2^31 = 28678040^169 * 2^31
.word   16838301 // zeta^233 * 2^31 = 28678040^233 * 2^31
.word 1817503137 // zeta^137 * (q^(-1) mod 2^32) * 2^31 = 28678040^137 * 375649793 * 2^31
.word 4137110159 // zeta^201 * (q^(-1) mod 2^32) * 2^31 = 28678040^201 * 375649793 * 2^31
.word 3861070683 // zeta^169 * (q^(-1) mod 2^32) * 2^31 = 28678040^169 * 375649793 * 2^31
.word 1425879907 // zeta^233 * (q^(-1) mod 2^32) * 2^31 = 28678040^233 * 375649793 * 2^31
.word   40841465 // zeta^ 50 * 2^31 = 28678040^ 50 * 2^31
.word    3577749 // zeta^178 * 2^31 = 28678040^178 * 2^31
.word   33845545 // zeta^114 * 2^31 = 28678040^114 * 2^31
.word   19555165 // zeta^242 * 2^31 = 28678040^242 * 2^31
.word 3459680519 // zeta^ 50 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 50 * 375649793 * 2^31
.word  495008363 // zeta^178 * (q^(-1) mod 2^32) * 2^31 = 28678040^178 * 375649793 * 2^31
.word 1885546711 // zeta^114 * (q^(-1) mod 2^32) * 2^31 = 28678040^114 * 375649793 * 2^31
.word 3630382755 // zeta^242 * (q^(-1) mod 2^32) * 2^31 = 28678040^242 * 375649793 * 2^31
.word   62758213 // zeta^ 25 * 2^31 = 28678040^ 25 * 2^31
.word    8005843 // zeta^ 89 * 2^31 = 28678040^ 89 * 2^31
.word   51922779 // zeta^ 57 * 2^31 = 28678040^ 57 * 2^31
.word    7245689 // zeta^121 * 2^31 = 28678040^121 * 2^31
.word  124982459 // zeta^ 25 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 25 * 375649793 * 2^31
.word 2964460845 // zeta^ 89 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 89 * 375649793 * 2^31
.word 1042630309 // zeta^ 57 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 57 * 375649793 * 2^31
.word 3756534407 // zeta^121 * (q^(-1) mod 2^32) * 2^31 = 28678040^121 * 375649793 * 2^31
.word   30225471 // zeta^153 * 2^31 = 28678040^153 * 2^31
.word   44151511 // zeta^217 * 2^31 = 28678040^217 * 2^31
.word   64890121 // zeta^185 * 2^31 = 28678040^185 * 2^31
.word   65259669 // zeta^249 * 2^31 = 28678040^249 * 2^31
.word  799097281 // zeta^153 * (q^(-1) mod 2^32) * 2^31 = 28678040^153 * 375649793 * 2^31
.word 4166955817 // zeta^217 * (q^(-1) mod 2^32) * 2^31 = 28678040^217 * 375649793 * 2^31
.word 1813001463 // zeta^185 * (q^(-1) mod 2^32) * 2^31 = 28678040^185 * 375649793 * 2^31
.word 3116239211 // zeta^249 * (q^(-1) mod 2^32) * 2^31 = 28678040^249 * 375649793 * 2^31
.word   12974361 // zeta^ 10 * 2^31 = 28678040^ 10 * 2^31
.word   41807515 // zeta^138 * 2^31 = 28678040^138 * 2^31
.word   56379967 // zeta^ 74 * 2^31 = 28678040^ 74 * 2^31
.word   13380915 // zeta^202 * 2^31 = 28678040^202 * 2^31
.word 1194393831 // zeta^ 10 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 10 * 375649793 * 2^31
.word 1648893797 // zeta^138 * (q^(-1) mod 2^32) * 2^31 = 28678040^138 * 375649793 * 2^31
.word  753806273 // zeta^ 74 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 74 * 375649793 * 2^31
.word 4010528973 // zeta^202 * (q^(-1) mod 2^32) * 2^31 = 28678040^202 * 375649793 * 2^31
.word   16772797 // zeta^  5 * 2^31 = 28678040^  5 * 2^31
.word   58675875 // zeta^ 69 * 2^31 = 28678040^ 69 * 2^31
.word   59974505 // zeta^ 37 * 2^31 = 28678040^ 37 * 2^31
.word   33980107 // zeta^101 * 2^31 = 28678040^101 * 2^31
.word 2122281795 // zeta^  5 * (q^(-1) mod 2^32) * 2^31 = 28678040^  5 * 375649793 * 2^31
.word 2886667101 // zeta^ 69 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 69 * 375649793 * 2^31
.word 3771397783 // zeta^ 37 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 37 * 375649793 * 2^31
.word 1168207669 // zeta^101 * (q^(-1) mod 2^32) * 2^31 = 28678040^101 * 375649793 * 2^31
.word   28448893 // zeta^133 * 2^31 = 28678040^133 * 2^31
.word   24378249 // zeta^197 * 2^31 = 28678040^197 * 2^31
.word   62687027 // zeta^165 * 2^31 = 28678040^165 * 2^31
.word   65645595 // zeta^229 * 2^31 = 28678040^229 * 2^31
.word 2758723971 // zeta^133 * (q^(-1) mod 2^32) * 2^31 = 28678040^133 * 375649793 * 2^31
.word 2128305783 // zeta^197 * (q^(-1) mod 2^32) * 2^31 = 28678040^197 * 375649793 * 2^31
.word  664762061 // zeta^165 * (q^(-1) mod 2^32) * 2^31 = 28678040^165 * 375649793 * 2^31
.word 2420335077 // zeta^229 * (q^(-1) mod 2^32) * 2^31 = 28678040^229 * 375649793 * 2^31
.word   52771617 // zeta^ 42 * 2^31 = 28678040^ 42 * 2^31
.word   23396495 // zeta^170 * 2^31 = 28678040^170 * 2^31
.word   51483005 // zeta^106 * 2^31 = 28678040^106 * 2^31
.word   11487943 // zeta^234 * 2^31 = 28678040^234 * 2^31
.word 2185629407 // zeta^ 42 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 42 * 375649793 * 2^31
.word 1858377073 // zeta^170 * (q^(-1) mod 2^32) * 2^31 = 28678040^170 * 375649793 * 2^31
.word  432623747 // zeta^106 * (q^(-1) mod 2^32) * 2^31 = 28678040^106 * 375649793 * 2^31
.word 2290121529 // zeta^234 * (q^(-1) mod 2^32) * 2^31 = 28678040^234 * 375649793 * 2^31
.word   63287737 // zeta^ 21 * 2^31 = 28678040^ 21 * 2^31
.word   56338313 // zeta^ 85 * 2^31 = 28678040^ 85 * 2^31
.word   19445427 // zeta^ 53 * 2^31 = 28678040^ 53 * 2^31
.word   29167561 // zeta^117 * 2^31 = 28678040^117 * 2^31
.word 1659340871 // zeta^ 21 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 21 * 375649793 * 2^31
.word 1504424567 // zeta^ 85 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 85 * 375649793 * 2^31
.word 3591259981 // zeta^ 53 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 53 * 375649793 * 2^31
.word 4032612919 // zeta^117 * (q^(-1) mod 2^32) * 2^31 = 28678040^117 * 375649793 * 2^31
.word    7740335 // zeta^149 * 2^31 = 28678040^149 * 2^31
.word   23515783 // zeta^213 * 2^31 = 28678040^213 * 2^31
.word   33583453 // zeta^181 * 2^31 = 28678040^181 * 2^31
.word   60337403 // zeta^245 * 2^31 = 28678040^245 * 2^31
.word 3259152977 // zeta^149 * (q^(-1) mod 2^32) * 2^31 = 28678040^149 * 375649793 * 2^31
.word  739668857 // zeta^213 * (q^(-1) mod 2^32) * 2^31 = 28678040^213 * 375649793 * 2^31
.word 3155767459 // zeta^181 * (q^(-1) mod 2^32) * 2^31 = 28678040^181 * 375649793 * 2^31
.word 3474792709 // zeta^245 * (q^(-1) mod 2^32) * 2^31 = 28678040^245 * 375649793 * 2^31
.word   35192755 // zeta^ 26 * 2^31 = 28678040^ 26 * 2^31
.word   36544119 // zeta^154 * 2^31 = 28678040^154 * 2^31
.word    6787663 // zeta^ 90 * 2^31 = 28678040^ 90 * 2^31
.word   63484749 // zeta^218 * 2^31 = 28678040^218 * 2^31
.word 3019374157 // zeta^ 26 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 26 * 375649793 * 2^31
.word 2777089929 // zeta^154 * (q^(-1) mod 2^32) * 2^31 = 28678040^154 * 375649793 * 2^31
.word  443777969 // zeta^ 90 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 90 * 375649793 * 2^31
.word  723799731 // zeta^218 * (q^(-1) mod 2^32) * 2^31 = 28678040^218 * 375649793 * 2^31
.word   61997615 // zeta^ 13 * 2^31 = 28678040^ 13 * 2^31
.word    4479011 // zeta^ 77 * 2^31 = 28678040^ 77 * 2^31
.word   38089877 // zeta^ 45 * 2^31 = 28678040^ 45 * 2^31
.word   16590903 // zeta^109 * 2^31 = 28678040^109 * 2^31
.word  201839569 // zeta^ 13 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 13 * 375649793 * 2^31
.word  998311389 // zeta^ 77 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 77 * 375649793 * 2^31
.word 1502911851 // zeta^ 45 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 45 * 375649793 * 2^31
.word 1931017673 // zeta^109 * (q^(-1) mod 2^32) * 2^31 = 28678040^109 * 375649793 * 2^31
.word   43852787 // zeta^141 * 2^31 = 28678040^141 * 2^31
.word   24597857 // zeta^205 * 2^31 = 28678040^205 * 2^31
.word   43936833 // zeta^173 * 2^31 = 28678040^173 * 2^31
.word   15636061 // zeta^237 * 2^31 = 28678040^237 * 2^31
.word  870210061 // zeta^141 * (q^(-1) mod 2^32) * 2^31 = 28678040^141 * 375649793 * 2^31
.word 4160386207 // zeta^205 * (q^(-1) mod 2^32) * 2^31 = 28678040^205 * 375649793 * 2^31
.word 1312300479 // zeta^173 * (q^(-1) mod 2^32) * 2^31 = 28678040^173 * 375649793 * 2^31
.word 2526874531 // zeta^237 * (q^(-1) mod 2^32) * 2^31 = 28678040^237 * 375649793 * 2^31
.word   55869129 // zeta^ 58 * 2^31 = 28678040^ 58 * 2^31
.word   16038683 // zeta^186 * 2^31 = 28678040^186 * 2^31
.word   43560065 // zeta^122 * 2^31 = 28678040^122 * 2^31
.word   25949329 // zeta^250 * 2^31 = 28678040^250 * 2^31
.word 2098944823 // zeta^ 58 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 58 * 375649793 * 2^31
.word  634278629 // zeta^186 * (q^(-1) mod 2^32) * 2^31 = 28678040^186 * 375649793 * 2^31
.word 2076204415 // zeta^122 * (q^(-1) mod 2^32) * 2^31 = 28678040^122 * 375649793 * 2^31
.word 2002629999 // zeta^250 * (q^(-1) mod 2^32) * 2^31 = 28678040^250 * 375649793 * 2^31
.word    6591765 // zeta^ 29 * 2^31 = 28678040^ 29 * 2^31
.word    1696249 // zeta^ 93 * 2^31 = 28678040^ 93 * 2^31
.word   21795289 // zeta^ 61 * 2^31 = 28678040^ 61 * 2^31
.word   17734591 // zeta^125 * 2^31 = 28678040^125 * 2^31
.word 3812244715 // zeta^ 29 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 29 * 375649793 * 2^31
.word 1467340807 // zeta^ 93 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 93 * 375649793 * 2^31
.word 1570891815 // zeta^ 61 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 61 * 375649793 * 2^31
.word 1349179969 // zeta^125 * (q^(-1) mod 2^32) * 2^31 = 28678040^125 * 375649793 * 2^31
.word   66853037 // zeta^157 * 2^31 = 28678040^157 * 2^31
.word   24930199 // zeta^221 * 2^31 = 28678040^221 * 2^31
.word   54854635 // zeta^189 * 2^31 = 28678040^189 * 2^31
.word   39952565 // zeta^253 * 2^31 = 28678040^253 * 2^31
.word 1399236947 // zeta^157 * (q^(-1) mod 2^32) * 2^31 = 28678040^157 * 375649793 * 2^31
.word 1771273833 // zeta^221 * (q^(-1) mod 2^32) * 2^31 = 28678040^221 * 375649793 * 2^31
.word 4111870485 // zeta^189 * (q^(-1) mod 2^32) * 2^31 = 28678040^189 * 375649793 * 2^31
.word 2033283403 // zeta^253 * (q^(-1) mod 2^32) * 2^31 = 28678040^253 * 375649793 * 2^31
.word    5623923 // zeta^  6 * 2^31 = 28678040^  6 * 2^31
.word   38701067 // zeta^134 * 2^31 = 28678040^134 * 2^31
.word   18571677 // zeta^ 70 * 2^31 = 28678040^ 70 * 2^31
.word   14491707 // zeta^198 * 2^31 = 28678040^198 * 2^31
.word  182627725 // zeta^  6 * (q^(-1) mod 2^32) * 2^31 = 28678040^  6 * 375649793 * 2^31
.word 4172670453 // zeta^134 * (q^(-1) mod 2^32) * 2^31 = 28678040^134 * 375649793 * 2^31
.word 1902166115 // zeta^ 70 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 70 * 375649793 * 2^31
.word 4183371205 // zeta^198 * (q^(-1) mod 2^32) * 2^31 = 28678040^198 * 375649793 * 2^31
.word   17941849 // zeta^  3 * 2^31 = 28678040^  3 * 2^31
.word   12982967 // zeta^ 67 * 2^31 = 28678040^ 67 * 2^31
.word    8061707 // zeta^ 35 * 2^31 = 28678040^ 35 * 2^31
.word   17774995 // zeta^ 99 * 2^31 = 28678040^ 99 * 2^31
.word 4091524263 // zeta^  3 * (q^(-1) mod 2^32) * 2^31 = 28678040^  3 * 375649793 * 2^31
.word 2462649161 // zeta^ 67 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 67 * 375649793 * 2^31
.word 2874632949 // zeta^ 35 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 35 * 375649793 * 2^31
.word 2009367661 // zeta^ 99 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 99 * 375649793 * 2^31
.word   61107981 // zeta^131 * 2^31 = 28678040^131 * 2^31
.word   38975641 // zeta^195 * 2^31 = 28678040^195 * 2^31
.word   40352225 // zeta^163 * 2^31 = 28678040^163 * 2^31
.word   49569327 // zeta^227 * 2^31 = 28678040^227 * 2^31
.word 3919450867 // zeta^131 * (q^(-1) mod 2^32) * 2^31 = 28678040^131 * 375649793 * 2^31
.word 4146020711 // zeta^195 * (q^(-1) mod 2^32) * 2^31 = 28678040^195 * 375649793 * 2^31
.word  418844703 // zeta^163 * (q^(-1) mod 2^32) * 2^31 = 28678040^163 * 375649793 * 2^31
.word 3026024401 // zeta^227 * (q^(-1) mod 2^32) * 2^31 = 28678040^227 * 375649793 * 2^31
.word   26799603 // zeta^ 38 * 2^31 = 28678040^ 38 * 2^31
.word   33463463 // zeta^166 * 2^31 = 28678040^166 * 2^31
.word   39332725 // zeta^102 * 2^31 = 28678040^102 * 2^31
.word   61125067 // zeta^230 * 2^31 = 28678040^230 * 2^31
.word  583438349 // zeta^ 38 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 38 * 375649793 * 2^31
.word 1692658009 // zeta^166 * (q^(-1) mod 2^32) * 2^31 = 28678040^166 * 375649793 * 2^31
.word 1738958475 // zeta^102 * (q^(-1) mod 2^32) * 2^31 = 28678040^102 * 375649793 * 2^31
.word 2248227893 // zeta^230 * (q^(-1) mod 2^32) * 2^31 = 28678040^230 * 375649793 * 2^31
.word   40014327 // zeta^ 19 * 2^31 = 28678040^ 19 * 2^31
.word     562885 // zeta^ 83 * 2^31 = 28678040^ 83 * 2^31
.word   51009393 // zeta^ 51 * 2^31 = 28678040^ 51 * 2^31
.word   51995259 // zeta^115 * 2^31 = 28678040^115 * 2^31
.word 2564101129 // zeta^ 19 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 19 * 375649793 * 2^31
.word 2196183867 // zeta^ 83 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 83 * 375649793 * 2^31
.word 2252083855 // zeta^ 51 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 51 * 375649793 * 2^31
.word 4038290309 // zeta^115 * (q^(-1) mod 2^32) * 2^31 = 28678040^115 * 375649793 * 2^31
.word   24330211 // zeta^147 * 2^31 = 28678040^147 * 2^31
.word    7682101 // zeta^211 * 2^31 = 28678040^211 * 2^31
.word    7401943 // zeta^179 * 2^31 = 28678040^179 * 2^31
.word   41757453 // zeta^243 * 2^31 = 28678040^243 * 2^31
.word  140484125 // zeta^147 * (q^(-1) mod 2^32) * 2^31 = 28678040^147 * 375649793 * 2^31
.word  285792715 // zeta^211 * (q^(-1) mod 2^32) * 2^31 = 28678040^211 * 375649793 * 2^31
.word 1996846121 // zeta^179 * (q^(-1) mod 2^32) * 2^31 = 28678040^179 * 375649793 * 2^31
.word 4062009075 // zeta^243 * (q^(-1) mod 2^32) * 2^31 = 28678040^243 * 375649793 * 2^31
.word   65375453 // zeta^ 22 * 2^31 = 28678040^ 22 * 2^31
.word   40797001 // zeta^150 * 2^31 = 28678040^150 * 2^31
.word   59835311 // zeta^ 86 * 2^31 = 28678040^ 86 * 2^31
.word   32875577 // zeta^214 * 2^31 = 28678040^214 * 2^31
.word 4014413091 // zeta^ 22 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 22 * 375649793 * 2^31
.word 3224262327 // zeta^150 * (q^(-1) mod 2^32) * 2^31 = 28678040^150 * 375649793 * 2^31
.word  741855825 // zeta^ 86 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 86 * 375649793 * 2^31
.word 2318439879 // zeta^214 * (q^(-1) mod 2^32) * 2^31 = 28678040^214 * 375649793 * 2^31
.word   10045293 // zeta^ 11 * 2^31 = 28678040^ 11 * 2^31
.word   53076657 // zeta^ 75 * 2^31 = 28678040^ 75 * 2^31
.word   17896617 // zeta^ 43 * 2^31 = 28678040^ 43 * 2^31
.word   58413331 // zeta^107 * 2^31 = 28678040^107 * 2^31
.word 3080518291 // zeta^ 11 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 11 * 375649793 * 2^31
.word 3700229967 // zeta^ 75 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 75 * 375649793 * 2^31
.word  297370967 // zeta^ 43 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 43 * 375649793 * 2^31
.word 2151902445 // zeta^107 * (q^(-1) mod 2^32) * 2^31 = 28678040^107 * 375649793 * 2^31
.word   19472551 // zeta^139 * 2^31 = 28678040^139 * 2^31
.word    6043561 // zeta^203 * 2^31 = 28678040^203 * 2^31
.word   20934449 // zeta^171 * 2^31 = 28678040^171 * 2^31
.word   37620445 // zeta^235 * 2^31 = 28678040^235 * 2^31
.word 2128700761 // zeta^139 * (q^(-1) mod 2^32) * 2^31 = 28678040^139 * 375649793 * 2^31
.word 1439457879 // zeta^203 * (q^(-1) mod 2^32) * 2^31 = 28678040^203 * 375649793 * 2^31
.word 3556014799 // zeta^171 * (q^(-1) mod 2^32) * 2^31 = 28678040^171 * 375649793 * 2^31
.word  769300259 // zeta^235 * (q^(-1) mod 2^32) * 2^31 = 28678040^235 * 375649793 * 2^31
.word   12921459 // zeta^ 54 * 2^31 = 28678040^ 54 * 2^31
.word   63769677 // zeta^182 * 2^31 = 28678040^182 * 2^31
.word   61505033 // zeta^118 * 2^31 = 28678040^118 * 2^31
.word   65692461 // zeta^246 * 2^31 = 28678040^246 * 2^31
.word 1006064525 // zeta^ 54 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 54 * 375649793 * 2^31
.word 2459563443 // zeta^182 * (q^(-1) mod 2^32) * 2^31 = 28678040^182 * 375649793 * 2^31
.word 2747128823 // zeta^118 * (q^(-1) mod 2^32) * 2^31 = 28678040^118 * 375649793 * 2^31
.word 2288082643 // zeta^246 * (q^(-1) mod 2^32) * 2^31 = 28678040^246 * 375649793 * 2^31
.word   20171011 // zeta^ 27 * 2^31 = 28678040^ 27 * 2^31
.word   36495001 // zeta^ 91 * 2^31 = 28678040^ 91 * 2^31
.word   62685175 // zeta^ 59 * 2^31 = 28678040^ 59 * 2^31
.word     664745 // zeta^123 * 2^31 = 28678040^123 * 2^31
.word 1031427325 // zeta^ 27 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 27 * 375649793 * 2^31
.word 2764118887 // zeta^ 91 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 91 * 375649793 * 2^31
.word  583476745 // zeta^ 59 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 59 * 375649793 * 2^31
.word 2371908951 // zeta^123 * (q^(-1) mod 2^32) * 2^31 = 28678040^123 * 375649793 * 2^31
.word   56713759 // zeta^155 * 2^31 = 28678040^155 * 2^31
.word   59594509 // zeta^219 * 2^31 = 28678040^219 * 2^31
.word   41235703 // zeta^187 * 2^31 = 28678040^187 * 2^31
.word   11581499 // zeta^251 * 2^31 = 28678040^251 * 2^31
.word 3428307937 // zeta^155 * (q^(-1) mod 2^32) * 2^31 = 28678040^155 * 375649793 * 2^31
.word 1657088755 // zeta^219 * (q^(-1) mod 2^32) * 2^31 = 28678040^219 * 375649793 * 2^31
.word 2803921161 // zeta^187 * (q^(-1) mod 2^32) * 2^31 = 28678040^187 * 375649793 * 2^31
.word 3715470789 // zeta^251 * (q^(-1) mod 2^32) * 2^31 = 28678040^251 * 375649793 * 2^31
.word   23458751 // zeta^ 14 * 2^31 = 28678040^ 14 * 2^31
.word    9406759 // zeta^142 * 2^31 = 28678040^142 * 2^31
.word   33711991 // zeta^ 78 * 2^31 = 28678040^ 78 * 2^31
.word   32167773 // zeta^206 * 2^31 = 28678040^206 * 2^31
.word 1501790785 // zeta^ 14 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 14 * 375649793 * 2^31
.word 2911894745 // zeta^142 * (q^(-1) mod 2^32) * 2^31 = 28678040^142 * 375649793 * 2^31
.word 1905016457 // zeta^ 78 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 78 * 375649793 * 2^31
.word  204130979 // zeta^206 * (q^(-1) mod 2^32) * 2^31 = 28678040^206 * 375649793 * 2^31
.word   26043621 // zeta^  7 * 2^31 = 28678040^  7 * 2^31
.word   51942461 // zeta^ 71 * 2^31 = 28678040^ 71 * 2^31
.word   14401009 // zeta^ 39 * 2^31 = 28678040^ 39 * 2^31
.word   60574133 // zeta^103 * 2^31 = 28678040^103 * 2^31
.word 1827638555 // zeta^  7 * (q^(-1) mod 2^32) * 2^31 = 28678040^  7 * 375649793 * 2^31
.word 3437088195 // zeta^ 71 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 71 * 375649793 * 2^31
.word 2892737551 // zeta^ 39 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 39 * 375649793 * 2^31
.word 3197159499 // zeta^103 * (q^(-1) mod 2^32) * 2^31 = 28678040^103 * 375649793 * 2^31
.word   16031087 // zeta^135 * 2^31 = 28678040^135 * 2^31
.word   25566271 // zeta^199 * 2^31 = 28678040^199 * 2^31
.word   54040269 // zeta^167 * 2^31 = 28678040^167 * 2^31
.word   36895029 // zeta^231 * 2^31 = 28678040^231 * 2^31
.word 2211821713 // zeta^135 * (q^(-1) mod 2^32) * 2^31 = 28678040^135 * 375649793 * 2^31
.word 3371719105 // zeta^199 * (q^(-1) mod 2^32) * 2^31 = 28678040^199 * 375649793 * 2^31
.word 2895604531 // zeta^167 * (q^(-1) mod 2^32) * 2^31 = 28678040^167 * 375649793 * 2^31
.word  349509835 // zeta^231 * (q^(-1) mod 2^32) * 2^31 = 28678040^231 * 375649793 * 2^31
.word   41803191 // zeta^ 46 * 2^31 = 28678040^ 46 * 2^31
.word   19377381 // zeta^174 * 2^31 = 28678040^174 * 2^31
.word    9664027 // zeta^110 * 2^31 = 28678040^110 * 2^31
.word   55794235 // zeta^238 * 2^31 = 28678040^238 * 2^31
.word 2460960841 // zeta^ 46 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 46 * 375649793 * 2^31
.word 1411728667 // zeta^174 * (q^(-1) mod 2^32) * 2^31 = 28678040^174 * 375649793 * 2^31
.word 1300076517 // zeta^110 * (q^(-1) mod 2^32) * 2^31 = 28678040^110 * 375649793 * 2^31
.word 3978752965 // zeta^238 * (q^(-1) mod 2^32) * 2^31 = 28678040^238 * 375649793 * 2^31
.word   19675339 // zeta^ 23 * 2^31 = 28678040^ 23 * 2^31
.word   21359151 // zeta^ 87 * 2^31 = 28678040^ 87 * 2^31
.word   63140729 // zeta^ 55 * 2^31 = 28678040^ 55 * 2^31
.word   23160723 // zeta^119 * 2^31 = 28678040^119 * 2^31
.word  398439733 // zeta^ 23 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 23 * 375649793 * 2^31
.word  897838033 // zeta^ 87 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 87 * 375649793 * 2^31
.word  494618247 // zeta^ 55 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 55 * 375649793 * 2^31
.word 3040761453 // zeta^119 * (q^(-1) mod 2^32) * 2^31 = 28678040^119 * 375649793 * 2^31
.word    9258847 // zeta^151 * 2^31 = 28678040^151 * 2^31
.word    4669959 // zeta^215 * 2^31 = 28678040^215 * 2^31
.word   41266143 // zeta^183 * 2^31 = 28678040^183 * 2^31
.word   61464071 // zeta^247 * 2^31 = 28678040^247 * 2^31
.word 2032168609 // zeta^151 * (q^(-1) mod 2^32) * 2^31 = 28678040^151 * 375649793 * 2^31
.word 1670448121 // zeta^215 * (q^(-1) mod 2^32) * 2^31 = 28678040^215 * 375649793 * 2^31
.word 1227164193 // zeta^183 * (q^(-1) mod 2^32) * 2^31 = 28678040^183 * 375649793 * 2^31
.word 1246128121 // zeta^247 * (q^(-1) mod 2^32) * 2^31 = 28678040^247 * 375649793 * 2^31
.word   43355169 // zeta^ 30 * 2^31 = 28678040^ 30 * 2^31
.word    5591977 // zeta^158 * 2^31 = 28678040^158 * 2^31
.word   40694335 // zeta^ 94 * 2^31 = 28678040^ 94 * 2^31
.word   25071607 // zeta^222 * 2^31 = 28678040^222 * 2^31
.word 1107279327 // zeta^ 30 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 30 * 375649793 * 2^31
.word  552289879 // zeta^158 * (q^(-1) mod 2^32) * 2^31 = 28678040^158 * 375649793 * 2^31
.word  879592385 // zeta^ 94 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 94 * 375649793 * 2^31
.word 2040862217 // zeta^222 * (q^(-1) mod 2^32) * 2^31 = 28678040^222 * 375649793 * 2^31
.word   34737117 // zeta^ 15 * 2^31 = 28678040^ 15 * 2^31
.word   45994147 // zeta^ 79 * 2^31 = 28678040^ 79 * 2^31
.word   42273719 // zeta^ 47 * 2^31 = 28678040^ 47 * 2^31
.word   60428681 // zeta^111 * 2^31 = 28678040^111 * 2^31
.word  303076899 // zeta^ 15 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 15 * 375649793 * 2^31
.word 3854339421 // zeta^ 79 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 79 * 375649793 * 2^31
.word 3799259721 // zeta^ 47 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 47 * 375649793 * 2^31
.word 1636911223 // zeta^111 * (q^(-1) mod 2^32) * 2^31 = 28678040^111 * 375649793 * 2^31
.word   26028927 // zeta^143 * 2^31 = 28678040^143 * 2^31
.word   64083527 // zeta^207 * 2^31 = 28678040^207 * 2^31
.word   60382541 // zeta^175 * 2^31 = 28678040^175 * 2^31
.word   31337387 // zeta^239 * 2^31 = 28678040^239 * 2^31
.word 2592721537 // zeta^143 * (q^(-1) mod 2^32) * 2^31 = 28678040^143 * 375649793 * 2^31
.word 1624305593 // zeta^207 * (q^(-1) mod 2^32) * 2^31 = 28678040^207 * 375649793 * 2^31
.word 3925320883 // zeta^175 * (q^(-1) mod 2^32) * 2^31 = 28678040^175 * 375649793 * 2^31
.word 3943334485 // zeta^239 * (q^(-1) mod 2^32) * 2^31 = 28678040^239 * 375649793 * 2^31
.word   27553395 // zeta^ 62 * 2^31 = 28678040^ 62 * 2^31
.word    7648471 // zeta^190 * 2^31 = 28678040^190 * 2^31
.word     689375 // zeta^126 * 2^31 = 28678040^126 * 2^31
.word   46555773 // zeta^254 * 2^31 = 28678040^254 * 2^31
.word 1673531277 // zeta^ 62 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 62 * 375649793 * 2^31
.word 1889513769 // zeta^190 * (q^(-1) mod 2^32) * 2^31 = 28678040^190 * 375649793 * 2^31
.word 1477062945 // zeta^126 * (q^(-1) mod 2^32) * 2^31 = 28678040^126 * 375649793 * 2^31
.word 2252242819 // zeta^254 * (q^(-1) mod 2^32) * 2^31 = 28678040^254 * 375649793 * 2^31
.word   15797163 // zeta^ 31 * 2^31 = 28678040^ 31 * 2^31
.word   40170027 // zeta^ 95 * 2^31 = 28678040^ 95 * 2^31
.word   10866061 // zeta^ 63 * 2^31 = 28678040^ 63 * 2^31
.word   56298001 // zeta^127 * 2^31 = 28678040^127 * 2^31
.word  683123285 // zeta^ 31 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 31 * 375649793 * 2^31
.word 2755967957 // zeta^ 95 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 95 * 375649793 * 2^31
.word  273527923 // zeta^ 63 * (q^(-1) mod 2^32) * 2^31 = 28678040^ 63 * 375649793 * 2^31
.word  644194287 // zeta^127 * (q^(-1) mod 2^32) * 2^31 = 28678040^127 * 375649793 * 2^31
.word   50400667 // zeta^159 * 2^31 = 28678040^159 * 2^31
.word   33861863 // zeta^223 * 2^31 = 28678040^223 * 2^31
.word   53736885 // zeta^191 * 2^31 = 28678040^191 * 2^31
.word   31774129 // zeta^255 * 2^31 = 28678040^255 * 2^31
.word 1694171237 // zeta^159 * (q^(-1) mod 2^32) * 2^31 = 28678040^159 * 375649793 * 2^31
.word  950555929 // zeta^223 * (q^(-1) mod 2^32) * 2^31 = 28678040^223 * 375649793 * 2^31
.word 2075204683 // zeta^191 * (q^(-1) mod 2^32) * 2^31 = 28678040^191 * 375649793 * 2^31
.word 1062212687 // zeta^255 * (q^(-1) mod 2^32) * 2^31 = 28678040^255 * 375649793 * 2^31
.text

// Montgomery multiplication via rounding
.macro mulmod dst, src, const, const_twisted
        vqrdmulh.s32   \dst,  \src, \const
        vmul.u32       \src,  \src, \const_twisted
        vqrdmlah.s32   \dst,  \src, modulus
.endm

.macro ct_butterfly a, b, root, root_twisted
        mulmod tmp, \b, \root, \root_twisted
        vsub.u32       \b,    \a, tmp
        vadd.u32       \a,    \a, tmp
.endm

.align 4
roots_addr: .word roots
.syntax unified
.type ntt_n256_u32_33556993_28678040_complete_manual, %function
.global ntt_n256_u32_33556993_28678040_complete_manual
ntt_n256_u32_33556993_28678040_complete_manual:

        push {r4-r11,lr}
        // Save MVE vector registers
        vpush {d8-d15}

        modulus  .req r12
        root_ptr .req r11

        .equ modulus_const, 33556993
        movw modulus, #:lower16:modulus_const
        movt modulus, #:upper16:modulus_const
        ldr  root_ptr, roots_addr

        in_low       .req r0
        in_high      .req r1

        add in_high, in_low, #(4*128)

        root0         .req r2
        root0_twisted .req r3
        root1         .req r4
        root1_twisted .req r5
        root2         .req r6
        root2_twisted .req r7

        data0 .req q0
        data1 .req q1
        data2 .req q2
        data3 .req q3

        tmp .req q4

        // Layers 1-2

        ldrd root0, root0_twisted, [root_ptr], #+8
        ldrd root1, root1_twisted, [root_ptr], #+8
        ldrd root2, root2_twisted, [root_ptr], #+8

        mov lr, #16
layer12_loop:
        vldrw.u32 data0, [in_low]
        vldrw.u32 data1, [in_low, #(4*64)]
        vldrw.u32 data2, [in_high]
        vldrw.u32 data3, [in_high, #(4*64)]

        ct_butterfly data0, data2, root0, root0_twisted
        ct_butterfly data1, data3, root0, root0_twisted
        ct_butterfly data0, data1, root1, root1_twisted
        ct_butterfly data2, data3, root2, root2_twisted

        vstrw.u32 data0, [in_low], #16
        vstrw.u32 data1, [in_low, #(4*64 - 16)]
        vstrw.u32 data2, [in_high], #16
        vstrw.u32 data3, [in_high, #(4*64-16)]

        le lr, layer12_loop
        .unreq in_high
        .unreq in_low

        in .req r0
        sub in, in, #(64*4)

        // Layers 3,4

        // 4 butterfly blocks per root config, 4 root configs
        // loop over root configs

        count .req r1
        mov count, #4

out_start:
        ldrd root0, root0_twisted, [root_ptr], #+8
        ldrd root1, root1_twisted, [root_ptr], #+8
        ldrd root2, root2_twisted, [root_ptr], #+8

        mov lr, #4
layer34_loop:
        vldrw.u32 data0, [in]
        vldrw.u32 data1, [in, #(4*1*16)]
        vldrw.u32 data2, [in, #(4*2*16)]
        vldrw.u32 data3, [in, #(4*3*16)]

        ct_butterfly data0, data2, root0, root0_twisted
        ct_butterfly data1, data3, root0, root0_twisted
        ct_butterfly data0, data1, root1, root1_twisted
        ct_butterfly data2, data3, root2, root2_twisted

        vstrw.u32 data0, [in], #16
        vstrw.u32 data1, [in, #(4*1*16 - 16)]
        vstrw.u32 data2, [in, #(4*2*16 - 16)]
        vstrw.u32 data3, [in, #(4*3*16 - 16)]

        le lr, layer34_loop

        add in, in, #(4*64 - 4*16)
        subs count, count, #1
        bne out_start

        sub in, in, #(4*256)

        // Layers 5,6

        // 1 butterfly blocks per root config, 16 root configs
        // loop over root configs

        mov lr, #16
layer56_loop:
        ldrd root0, root0_twisted, [root_ptr], #+24
        ldrd root1, root1_twisted, [root_ptr, #(-16)]
        ldrd root2, root2_twisted, [root_ptr, #(-8)]

        vldrw.u32 data0, [in]
        vldrw.u32 data1, [in, #(4*1*4)]
        vldrw.u32 data2, [in, #(4*2*4)]
        vldrw.u32 data3, [in, #(4*3*4)]

        ct_butterfly data0, data2, root0, root0_twisted
        ct_butterfly data1, data3, root0, root0_twisted
        ct_butterfly data0, data1, root1, root1_twisted
        ct_butterfly data2, data3, root2, root2_twisted

        vst40.u32 {data0, data1, data2, data3}, [in]
        vst41.u32 {data0, data1, data2, data3}, [in]
        vst42.u32 {data0, data1, data2, data3}, [in]
        vst43.u32 {data0, data1, data2, data3}, [in]!

        le lr, layer56_loop

        sub in, in, #(4*256)

        // Layers 7,8

        .unreq root0
        .unreq root0_twisted
        .unreq root1
        .unreq root1_twisted
        .unreq root2
        .unreq root2_twisted

        root0         .req q5
        root0_twisted .req q6
        root1         .req q5
        root1_twisted .req q6
        root2         .req q5
        root2_twisted .req q6

        mov lr, #16
layer78_loop:
        vldrw.u32 data0, [in], #64
        vldrw.u32 data1, [in, #-48]
        vldrw.u32 data2, [in, #-32]
        vldrw.u32 data3, [in, #-16]

        vldrw.u32 root0,         [root_ptr], #+96
        vldrw.u32 root0_twisted, [root_ptr, #(+16-96)]
        ct_butterfly data0, data2, root0, root0_twisted
        ct_butterfly data1, data3, root0, root0_twisted

        vldrw.u32 root1,         [root_ptr, #(32 - 96)]
        vldrw.u32 root1_twisted, [root_ptr, #(48 - 96)]
        ct_butterfly data0, data1, root1, root1_twisted

        vldrw.u32 root2,         [root_ptr, #(64-96)]
        vldrw.u32 root2_twisted, [root_ptr, #(80-96)]
        ct_butterfly data2, data3, root2, root2_twisted

        vstrw.u32 data0, [in, #(4*0*4 - 64)]
        vstrw.u32 data1, [in, #(4*1*4 - 64)]
        vstrw.u32 data2, [in, #(4*2*4 - 64)]
        vstrw.u32 data3, [in, #(4*3*4 - 64)]
        le lr, layer78_loop

        // Restore MVE vector registers
        vpop {d8-d15}
        // Restore GPRs
        pop {r4-r11,lr}
        bx lr