crypto/openssl/crypto/bn/asm/x86/comba.pl

   1 #!/usr/local/bin/perl
   2 # x86 assember
   3
   4 sub mul_add_c
   5         {
   6         local($a,$ai,$b,$bi,$c0,$c1,$c2,$pos,$i,$na,$nb)=@_;
   7
   8         # pos == -1 if eax and edx are pre-loaded, 0 to load from next
   9         # words, and 1 if load return value
  10
  11         &comment("mul a[$ai]*b[$bi]");
  12
  13         # "eax" and "edx" will always be pre-loaded.
  14         # &mov("eax",&DWP($ai*4,$a,"",0)) ;
  15         # &mov("edx",&DWP($bi*4,$b,"",0));
  16
  17         &mul("edx");
  18         &add($c0,"eax");
  19          &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 0;        # laod next a
  20          &mov("eax",&wparam(0)) if $pos > 0;                    # load r[]
  21          ###
  22         &adc($c1,"edx");
  23          &mov("edx",&DWP(($nb)*4,$b,"",0)) if $pos == 0;        # laod next b
  24          &mov("edx",&DWP(($nb)*4,$b,"",0)) if $pos == 1;        # laod next b
  25          ###
  26         &adc($c2,0);
  27          # is pos > 1, it means it is the last loop
  28          &mov(&DWP($i*4,"eax","",0),$c0) if $pos > 0;           # save r[];
  29         &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 1;         # laod next a
  30         }
  31
  32 sub sqr_add_c
  33         {
  34         local($r,$a,$ai,$bi,$c0,$c1,$c2,$pos,$i,$na,$nb)=@_;
  35
  36         # pos == -1 if eax and edx are pre-loaded, 0 to load from next
  37         # words, and 1 if load return value
  38
  39         &comment("sqr a[$ai]*a[$bi]");
  40
  41         # "eax" and "edx" will always be pre-loaded.
  42         # &mov("eax",&DWP($ai*4,$a,"",0)) ;
  43         # &mov("edx",&DWP($bi*4,$b,"",0));
  44
  45         if ($ai == $bi)
  46                 { &mul("eax");}
  47         else
  48                 { &mul("edx");}
  49         &add($c0,"eax");
  50          &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 0;        # load next a
  51          ###
  52         &adc($c1,"edx");
  53          &mov("edx",&DWP(($nb)*4,$a,"",0)) if ($pos == 1) && ($na != $nb);
  54          ###
  55         &adc($c2,0);
  56          # is pos > 1, it means it is the last loop
  57          &mov(&DWP($i*4,$r,"",0),$c0) if $pos > 0;              # save r[];
  58         &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 1;         # load next b
  59         }
  60
  61 sub sqr_add_c2
  62         {
  63         local($r,$a,$ai,$bi,$c0,$c1,$c2,$pos,$i,$na,$nb)=@_;
  64
  65         # pos == -1 if eax and edx are pre-loaded, 0 to load from next
  66         # words, and 1 if load return value
  67
  68         &comment("sqr a[$ai]*a[$bi]");
  69
  70         # "eax" and "edx" will always be pre-loaded.
  71         # &mov("eax",&DWP($ai*4,$a,"",0)) ;
  72         # &mov("edx",&DWP($bi*4,$a,"",0));
  73
  74         if ($ai == $bi)
  75                 { &mul("eax");}
  76         else
  77                 { &mul("edx");}
  78         &add("eax","eax");
  79          ###
  80         &adc("edx","edx");
  81          ###
  82         &adc($c2,0);
  83          &add($c0,"eax");
  84         &adc($c1,"edx");
  85          &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 0;        # load next a
  86          &mov("eax",&DWP(($na)*4,$a,"",0)) if $pos == 1;        # load next b
  87         &adc($c2,0);
  88         &mov(&DWP($i*4,$r,"",0),$c0) if $pos > 0;               # save r[];
  89          &mov("edx",&DWP(($nb)*4,$a,"",0)) if ($pos <= 1) && ($na != $nb);
  90          ###
  91         }
  92
  93 sub bn_mul_comba
  94         {
  95         local($name,$num)=@_;
  96         local($a,$b,$c0,$c1,$c2);
  97         local($i,$as,$ae,$bs,$be,$ai,$bi);
  98         local($tot,$end);
  99
 100         &function_begin_B($name,"");
 101
 102         $c0="ebx";
 103         $c1="ecx";
 104         $c2="ebp";
 105         $a="esi";
 106         $b="edi";
 107
 108         $as=0;
 109         $ae=0;
 110         $bs=0;
 111         $be=0;
 112         $tot=$num+$num-1;
 113
 114         &push("esi");
 115          &mov($a,&wparam(1));
 116         &push("edi");
 117          &mov($b,&wparam(2));
 118         &push("ebp");
 119          &push("ebx");
 120
 121         &xor($c0,$c0);
 122          &mov("eax",&DWP(0,$a,"",0));   # load the first word
 123         &xor($c1,$c1);
 124          &mov("edx",&DWP(0,$b,"",0));   # load the first second
 125
 126         for ($i=0; $i<$tot; $i++)
 127                 {
 128                 $ai=$as;
 129                 $bi=$bs;
 130                 $end=$be+1;
 131
 132                 &comment("################## Calculate word $i");
 133
 134                 for ($j=$bs; $j<$end; $j++)
 135                         {
 136                         &xor($c2,$c2) if ($j == $bs);
 137                         if (($j+1) == $end)
 138                                 {
 139                                 $v=1;
 140                                 $v=2 if (($i+1) == $tot);
 141                                 }
 142                         else
 143                                 { $v=0; }
 144                         if (($j+1) != $end)
 145                                 {
 146                                 $na=($ai-1);
 147                                 $nb=($bi+1);
 148                                 }
 149                         else
 150                                 {
 151                                 $na=$as+($i < ($num-1));
 152                                 $nb=$bs+($i >= ($num-1));
 153                                 }
 154 #printf STDERR "[$ai,$bi] -> [$na,$nb]\n";
 155                         &mul_add_c($a,$ai,$b,$bi,$c0,$c1,$c2,$v,$i,$na,$nb);
 156                         if ($v)
 157                                 {
 158                                 &comment("saved r[$i]");
 159                                 # &mov("eax",&wparam(0));
 160                                 # &mov(&DWP($i*4,"eax","",0),$c0);
 161                                 ($c0,$c1,$c2)=($c1,$c2,$c0);
 162                                 }
 163                         $ai--;
 164                         $bi++;
 165                         }
 166                 $as++ if ($i < ($num-1));
 167                 $ae++ if ($i >= ($num-1));
 168
 169                 $bs++ if ($i >= ($num-1));
 170                 $be++ if ($i < ($num-1));
 171                 }
 172         &comment("save r[$i]");
 173         # &mov("eax",&wparam(0));
 174         &mov(&DWP($i*4,"eax","",0),$c0);
 175
 176         &pop("ebx");
 177         &pop("ebp");
 178         &pop("edi");
 179         &pop("esi");
 180         &ret();
 181         &function_end_B($name);
 182         }
 183
 184 sub bn_sqr_comba
 185         {
 186         local($name,$num)=@_;
 187         local($r,$a,$c0,$c1,$c2)=@_;
 188         local($i,$as,$ae,$bs,$be,$ai,$bi);
 189         local($b,$tot,$end,$half);
 190
 191         &function_begin_B($name,"");
 192
 193         $c0="ebx";
 194         $c1="ecx";
 195         $c2="ebp";
 196         $a="esi";
 197         $r="edi";
 198
 199         &push("esi");
 200          &push("edi");
 201         &push("ebp");
 202          &push("ebx");
 203         &mov($r,&wparam(0));
 204          &mov($a,&wparam(1));
 205         &xor($c0,$c0);
 206          &xor($c1,$c1);
 207         &mov("eax",&DWP(0,$a,"",0)); # load the first word
 208
 209         $as=0;
 210         $ae=0;
 211         $bs=0;
 212         $be=0;
 213         $tot=$num+$num-1;
 214
 215         for ($i=0; $i<$tot; $i++)
 216                 {
 217                 $ai=$as;
 218                 $bi=$bs;
 219                 $end=$be+1;
 220
 221                 &comment("############### Calculate word $i");
 222                 for ($j=$bs; $j<$end; $j++)
 223                         {
 224                         &xor($c2,$c2) if ($j == $bs);
 225                         if (($ai-1) < ($bi+1))
 226                                 {
 227                                 $v=1;
 228                                 $v=2 if ($i+1) == $tot;
 229                                 }
 230                         else
 231                                 { $v=0; }
 232                         if (!$v)
 233                                 {
 234                                 $na=$ai-1;
 235                                 $nb=$bi+1;
 236                                 }
 237                         else
 238                                 {
 239                                 $na=$as+($i < ($num-1));
 240                                 $nb=$bs+($i >= ($num-1));
 241                                 }
 242                         if ($ai == $bi)
 243                                 {
 244                                 &sqr_add_c($r,$a,$ai,$bi,
 245                                         $c0,$c1,$c2,$v,$i,$na,$nb);
 246                                 }
 247                         else
 248                                 {
 249                                 &sqr_add_c2($r,$a,$ai,$bi,
 250                                         $c0,$c1,$c2,$v,$i,$na,$nb);
 251                                 }
 252                         if ($v)
 253                                 {
 254                                 &comment("saved r[$i]");
 255                                 #&mov(&DWP($i*4,$r,"",0),$c0);
 256                                 ($c0,$c1,$c2)=($c1,$c2,$c0);
 257                                 last;
 258                                 }
 259                         $ai--;
 260                         $bi++;
 261                         }
 262                 $as++ if ($i < ($num-1));
 263                 $ae++ if ($i >= ($num-1));
 264
 265                 $bs++ if ($i >= ($num-1));
 266                 $be++ if ($i < ($num-1));
 267                 }
 268         &mov(&DWP($i*4,$r,"",0),$c0);
 269         &pop("ebx");
 270         &pop("ebp");
 271         &pop("edi");
 272         &pop("esi");
 273         &ret();
 274         &function_end_B($name);
 275         }
 276
 277 1;