← Back to branch summary

~ubuntu-branches/ubuntu/saucy/fftw/saucy

« back to all changes in this revision

Viewing changes to fftw/ftw_32.c

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-29 01:48:33 UTC
  • mfrom: (5.1.1 sid)
  • Revision ID: package-import@ubuntu.com-20111129014833-3leawsns2nghyaoj
Tags: 2.1.5-1
http://bugs.debian.org/571446
http://bugs.debian.org/633175
* Team upload.
* New upstream release
* Package moved into the Debian Science team (no answer from the previous
  maintainer) and package not maintained.
* Standards-Version updated to version 3.9.2
* Vcs-Browser & Vcs-Svn updated
* Switch to mpi-default-dev (Closes: #571446)
* Get ride of .la files (Closes: #633175)
* Fix lintian warning debhelper-but-no-misc-depends
* Fix lintian warning patch-system-but-direct-changes-in-diff 
* Switch to dpkg-source 3.0 (quilt) format

Show diffs side-by-side

added added

removed removed

Lines of Context:
fftw/ftw_32.c
1
1
/*
2
 
 * Copyright (c) 1997-1999 Massachusetts Institute of Technology
 
2
 * Copyright (c) 1997-1999, 2003 Massachusetts Institute of Technology
3
3
 *
4
4
 * This program is free software; you can redistribute it and/or modify
5
5
 * it under the terms of the GNU General Public License as published by
18
18
 */
19
19
 
20
20
/* This file was automatically generated --- DO NOT EDIT */
21
 
/* Generated on Sun Nov  7 20:44:37 EST 1999 */
22
 
 
23
 
#include <fftw-int.h>
24
 
#include <fftw.h>
25
 
 
26
 
/* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -twiddle 32 */
 
21
/* Generated on Mon Mar 24 02:07:45 EST 2003 */
 
22
 
 
23
#include "fftw-int.h"
 
24
#include "fftw.h"
 
25
 
 
26
/* Generated by: /homee/stevenj/cvs/fftw/gensrc/genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -twiddle 32 */
27
27
 
28
28
/*
29
29
 * This function contains 434 FP additions, 208 FP multiplications,
30
30
 * (or, 340 additions, 114 multiplications, 94 fused multiply/add),
31
31
 * 90 stack variables, and 128 memory accesses
32
32
 */
33
 
static const fftw_real K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
34
 
static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
35
 
static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
36
 
static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
37
 
static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
38
 
static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
39
 
static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
 
33
static const fftw_real K195090322 =
 
34
FFTW_KONST(+0.195090322016128267848284868477022240927691618);
 
35
static const fftw_real K980785280 =
 
36
FFTW_KONST(+0.980785280403230449126182236134239036973933731);
 
37
static const fftw_real K831469612 =
 
38
FFTW_KONST(+0.831469612302545237078788377617905756738560812);
 
39
static const fftw_real K555570233 =
 
40
FFTW_KONST(+0.555570233019602224742830813948532874374937191);
 
41
static const fftw_real K382683432 =
 
42
FFTW_KONST(+0.382683432365089771728459984030398866761344562);
 
43
static const fftw_real K923879532 =
 
44
FFTW_KONST(+0.923879532511286756128183189396788286822416626);
 
45
static const fftw_real K707106781 =
 
46
FFTW_KONST(+0.707106781186547524400844362104849039284835938);
40
47
 
41
48
/*
42
49
 * Generator Id's : 
43
 
 * $Id: exprdag.ml,v 1.41 1999/05/26 15:44:14 fftw Exp $
44
 
 * $Id: fft.ml,v 1.43 1999/05/17 19:44:18 fftw Exp $
45
 
 * $Id: to_c.ml,v 1.25 1999/10/26 21:41:32 stevenj Exp $
 
50
 * $Id: exprdag.ml,v 1.43 2003/03/16 23:43:46 stevenj Exp $
 
51
 * $Id: fft.ml,v 1.44 2003/03/16 23:43:46 stevenj Exp $
 
52
 * $Id: to_c.ml,v 1.26 2003/03/16 23:43:46 stevenj Exp $
46
53
 */
47
54
 
48
 
void fftw_twiddle_32(fftw_complex *A, const fftw_complex *W, int iostride, int m, int dist)
 
55
void fftw_twiddle_32(fftw_complex *A, const fftw_complex *W, int iostride,
 
56
                     int m, int dist)
49
57
{
50
58
     int i;
51
59
     fftw_complex *inout;
1016
1024
                         tmp462 = K707106781 * (tmp397 - tmp396);
1017
1025
                         tmp372 = tmp370 - tmp371;
1018
1026
                         tmp377 = tmp373 - tmp376;
1019
 
                         tmp378 = (K923879532 * tmp372) + (K382683432 * tmp377);
1020
 
                         tmp392 = (K382683432 * tmp372) - (K923879532 * tmp377);
 
1027
                         tmp378 =
 
1028
                             (K923879532 * tmp372) + (K382683432 * tmp377);
 
1029
                         tmp392 =
 
1030
                             (K382683432 * tmp372) - (K923879532 * tmp377);
1021
1031
                    }
1022
1032
                    tmp457 = tmp455 + tmp456;
1023
1033
                    tmp463 = tmp456 - tmp455;
1033
1043
                         ASSERT_ALIGNED_DOUBLE;
1034
1044
                         tmp383 = tmp379 - tmp382;
1035
1045
                         tmp388 = tmp386 - tmp387;
1036
 
                         tmp389 = (K382683432 * tmp383) - (K923879532 * tmp388);
1037
 
                         tmp393 = (K382683432 * tmp388) + (K923879532 * tmp383);
 
1046
                         tmp389 =
 
1047
                             (K382683432 * tmp383) - (K923879532 * tmp388);
 
1048
                         tmp393 =
 
1049
                             (K382683432 * tmp388) + (K923879532 * tmp383);
1038
1050
                         tmp403 = tmp379 + tmp382;
1039
1051
                         tmp404 = tmp386 + tmp387;
1040
 
                         tmp405 = (K923879532 * tmp403) - (K382683432 * tmp404);
1041
 
                         tmp409 = (K923879532 * tmp404) + (K382683432 * tmp403);
 
1052
                         tmp405 =
 
1053
                             (K923879532 * tmp403) - (K382683432 * tmp404);
 
1054
                         tmp409 =
 
1055
                             (K923879532 * tmp404) + (K382683432 * tmp403);
1042
1056
                    }
1043
1057
               }
1044
1058
               {
1370
1384
}
1371
1385
 
1372
1386
static const int twiddle_order[] =
1373
 
{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31};
1374
 
fftw_codelet_desc fftw_twiddle_32_desc =
1375
 
{
 
1387
    { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
 
1388
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 };
 
1389
fftw_codelet_desc fftw_twiddle_32_desc = {
1376
1390
     "fftw_twiddle_32",
1377
1391
     (void (*)()) fftw_twiddle_32,
1378
1392
     32,