← Back to branch summary

~ubuntu-branches/ubuntu/quantal/fftw/quantal

« back to all changes in this revision

Viewing changes to rfftw/fhf_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:
rfftw/fhf_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:52 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 -hc2hc-forward 32 */
 
21
/* Generated on Mon Mar 24 02:08:13 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 -hc2hc-forward 32 */
27
27
 
28
28
/*
29
29
 * This function contains 764 FP additions, 340 FP multiplications,
30
30
 * (or, 618 additions, 194 multiplications, 146 fused multiply/add),
31
31
 * 91 stack variables, and 256 memory accesses
32
32
 */
33
 
static const fftw_real K145142338 = FFTW_KONST(+0.145142338627231183818096187908697637345738139);
34
 
static const fftw_real K478470167 = FFTW_KONST(+0.478470167866104432467898943490134984741424603);
35
 
static const fftw_real K235698368 = FFTW_KONST(+0.235698368412998824278193812952627188828730159);
36
 
static const fftw_real K440960632 = FFTW_KONST(+0.440960632174177514856378431830194174754221310);
37
 
static const fftw_real K317196642 = FFTW_KONST(+0.317196642081822749107585806612746685337843547);
38
 
static const fftw_real K386505226 = FFTW_KONST(+0.386505226681368480405453304879234900485520646);
39
 
static const fftw_real K497592363 = FFTW_KONST(+0.497592363336098443122418476554739960787737434);
40
 
static const fftw_real K049008570 = FFTW_KONST(+0.049008570164780300997097781944320922930568337);
41
 
static const fftw_real K277785116 = FFTW_KONST(+0.277785116509801112371415406974266437187468595);
42
 
static const fftw_real K415734806 = FFTW_KONST(+0.415734806151272618539394188808952878369280406);
43
 
static const fftw_real K097545161 = FFTW_KONST(+0.097545161008064133924142434238511120463845809);
44
 
static const fftw_real K490392640 = FFTW_KONST(+0.490392640201615224563091118067119518486966865);
45
 
static const fftw_real K765366864 = FFTW_KONST(+0.765366864730179543456919968060797733522689125);
46
 
static const fftw_real K1_847759065 = FFTW_KONST(+1.847759065022573512256366378793576573644833252);
47
 
static const fftw_real K1_414213562 = FFTW_KONST(+1.414213562373095048801688724209698078569671875);
48
 
static const fftw_real K2_000000000 = FFTW_KONST(+2.000000000000000000000000000000000000000000000);
49
 
static const fftw_real K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
50
 
static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
51
 
static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
52
 
static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
53
 
static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
54
 
static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
55
 
static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
 
33
static const fftw_real K145142338 =
 
34
FFTW_KONST(+0.145142338627231183818096187908697637345738139);
 
35
static const fftw_real K478470167 =
 
36
FFTW_KONST(+0.478470167866104432467898943490134984741424603);
 
37
static const fftw_real K235698368 =
 
38
FFTW_KONST(+0.235698368412998824278193812952627188828730159);
 
39
static const fftw_real K440960632 =
 
40
FFTW_KONST(+0.440960632174177514856378431830194174754221310);
 
41
static const fftw_real K317196642 =
 
42
FFTW_KONST(+0.317196642081822749107585806612746685337843547);
 
43
static const fftw_real K386505226 =
 
44
FFTW_KONST(+0.386505226681368480405453304879234900485520646);
 
45
static const fftw_real K497592363 =
 
46
FFTW_KONST(+0.497592363336098443122418476554739960787737434);
 
47
static const fftw_real K049008570 =
 
48
FFTW_KONST(+0.049008570164780300997097781944320922930568337);
 
49
static const fftw_real K277785116 =
 
50
FFTW_KONST(+0.277785116509801112371415406974266437187468595);
 
51
static const fftw_real K415734806 =
 
52
FFTW_KONST(+0.415734806151272618539394188808952878369280406);
 
53
static const fftw_real K097545161 =
 
54
FFTW_KONST(+0.097545161008064133924142434238511120463845809);
 
55
static const fftw_real K490392640 =
 
56
FFTW_KONST(+0.490392640201615224563091118067119518486966865);
 
57
static const fftw_real K765366864 =
 
58
FFTW_KONST(+0.765366864730179543456919968060797733522689125);
 
59
static const fftw_real K1_847759065 =
 
60
FFTW_KONST(+1.847759065022573512256366378793576573644833252);
 
61
static const fftw_real K1_414213562 =
 
62
FFTW_KONST(+1.414213562373095048801688724209698078569671875);
 
63
static const fftw_real K2_000000000 =
 
64
FFTW_KONST(+2.000000000000000000000000000000000000000000000);
 
65
static const fftw_real K195090322 =
 
66
FFTW_KONST(+0.195090322016128267848284868477022240927691618);
 
67
static const fftw_real K980785280 =
 
68
FFTW_KONST(+0.980785280403230449126182236134239036973933731);
 
69
static const fftw_real K555570233 =
 
70
FFTW_KONST(+0.555570233019602224742830813948532874374937191);
 
71
static const fftw_real K831469612 =
 
72
FFTW_KONST(+0.831469612302545237078788377617905756738560812);
 
73
static const fftw_real K382683432 =
 
74
FFTW_KONST(+0.382683432365089771728459984030398866761344562);
 
75
static const fftw_real K923879532 =
 
76
FFTW_KONST(+0.923879532511286756128183189396788286822416626);
 
77
static const fftw_real K707106781 =
 
78
FFTW_KONST(+0.707106781186547524400844362104849039284835938);
56
79
 
57
80
/*
58
81
 * Generator Id's : 
59
 
 * $Id: exprdag.ml,v 1.41 1999/05/26 15:44:14 fftw Exp $
60
 
 * $Id: fft.ml,v 1.43 1999/05/17 19:44:18 fftw Exp $
61
 
 * $Id: to_c.ml,v 1.25 1999/10/26 21:41:32 stevenj Exp $
 
82
 * $Id: exprdag.ml,v 1.43 2003/03/16 23:43:46 stevenj Exp $
 
83
 * $Id: fft.ml,v 1.44 2003/03/16 23:43:46 stevenj Exp $
 
84
 * $Id: to_c.ml,v 1.26 2003/03/16 23:43:46 stevenj Exp $
62
85
 */
63
86
 
64
 
void fftw_hc2hc_forward_32(fftw_real *A, const fftw_complex *W, int iostride, int m, int dist)
 
87
void fftw_hc2hc_forward_32(fftw_real *A, const fftw_complex *W,
 
88
                           int iostride, int m, int dist)
65
89
{
66
90
     int i;
67
91
     fftw_real *X;
1452
1476
                         tmp643 = K707106781 * (tmp579 - tmp578);
1453
1477
                         tmp554 = tmp552 - tmp553;
1454
1478
                         tmp559 = tmp555 - tmp558;
1455
 
                         tmp560 = (K923879532 * tmp554) + (K382683432 * tmp559);
1456
 
                         tmp574 = (K382683432 * tmp554) - (K923879532 * tmp559);
 
1479
                         tmp560 =
 
1480
                             (K923879532 * tmp554) + (K382683432 * tmp559);
 
1481
                         tmp574 =
 
1482
                             (K382683432 * tmp554) - (K923879532 * tmp559);
1457
1483
                    }
1458
1484
                    tmp638 = tmp636 + tmp637;
1459
1485
                    tmp644 = tmp637 - tmp636;
1469
1495
                         ASSERT_ALIGNED_DOUBLE;
1470
1496
                         tmp565 = tmp561 - tmp564;
1471
1497
                         tmp570 = tmp568 - tmp569;
1472
 
                         tmp571 = (K382683432 * tmp565) - (K923879532 * tmp570);
1473
 
                         tmp575 = (K382683432 * tmp570) + (K923879532 * tmp565);
 
1498
                         tmp571 =
 
1499
                             (K382683432 * tmp565) - (K923879532 * tmp570);
 
1500
                         tmp575 =
 
1501
                             (K382683432 * tmp570) + (K923879532 * tmp565);
1474
1502
                         tmp585 = tmp561 + tmp564;
1475
1503
                         tmp586 = tmp568 + tmp569;
1476
 
                         tmp587 = (K923879532 * tmp585) - (K382683432 * tmp586);
1477
 
                         tmp591 = (K923879532 * tmp586) + (K382683432 * tmp585);
 
1504
                         tmp587 =
 
1505
                             (K923879532 * tmp585) - (K382683432 * tmp586);
 
1506
                         tmp591 =
 
1507
                             (K923879532 * tmp586) + (K382683432 * tmp585);
1478
1508
                    }
1479
1509
               }
1480
1510
               {
1914
1944
                         ASSERT_ALIGNED_DOUBLE;
1915
1945
                         tmp67 = X[3 * iostride];
1916
1946
                         tmp68 = X[19 * iostride];
1917
 
                         tmp69 = (K1_847759065 * tmp67) - (K765366864 * tmp68);
1918
 
                         tmp75 = (K765366864 * tmp67) + (K1_847759065 * tmp68);
 
1947
                         tmp69 =
 
1948
                             (K1_847759065 * tmp67) - (K765366864 * tmp68);
 
1949
                         tmp75 =
 
1950
                             (K765366864 * tmp67) + (K1_847759065 * tmp68);
1919
1951
                         tmp70 = X[11 * iostride];
1920
1952
                         tmp71 = X[27 * iostride];
1921
 
                         tmp72 = (K765366864 * tmp70) - (K1_847759065 * tmp71);
1922
 
                         tmp76 = (K1_847759065 * tmp70) + (K765366864 * tmp71);
 
1953
                         tmp72 =
 
1954
                             (K765366864 * tmp70) - (K1_847759065 * tmp71);
 
1955
                         tmp76 =
 
1956
                             (K1_847759065 * tmp70) + (K765366864 * tmp71);
1923
1957
                    }
1924
1958
               }
1925
1959
               {
2059
2093
                         ASSERT_ALIGNED_DOUBLE;
2060
2094
                         tmp44 = X[5 * iostride];
2061
2095
                         tmp45 = X[21 * iostride];
2062
 
                         tmp46 = (K1_847759065 * tmp44) - (K765366864 * tmp45);
2063
 
                         tmp52 = (K765366864 * tmp44) + (K1_847759065 * tmp45);
 
2096
                         tmp46 =
 
2097
                             (K1_847759065 * tmp44) - (K765366864 * tmp45);
 
2098
                         tmp52 =
 
2099
                             (K765366864 * tmp44) + (K1_847759065 * tmp45);
2064
2100
                         tmp47 = X[13 * iostride];
2065
2101
                         tmp48 = X[29 * iostride];
2066
 
                         tmp49 = (K765366864 * tmp47) - (K1_847759065 * tmp48);
2067
 
                         tmp53 = (K1_847759065 * tmp47) + (K765366864 * tmp48);
 
2102
                         tmp49 =
 
2103
                             (K765366864 * tmp47) - (K1_847759065 * tmp48);
 
2104
                         tmp53 =
 
2105
                             (K1_847759065 * tmp47) + (K765366864 * tmp48);
2068
2106
                    }
2069
2107
               }
2070
2108
               {
2298
2336
}
2299
2337
 
2300
2338
static const int twiddle_order[] =
2301
 
{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};
2302
 
fftw_codelet_desc fftw_hc2hc_forward_32_desc =
2303
 
{
 
2339
    { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
 
2340
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 };
 
2341
fftw_codelet_desc fftw_hc2hc_forward_32_desc = {
2304
2342
     "fftw_hc2hc_forward_32",
2305
2343
     (void (*)()) fftw_hc2hc_forward_32,
2306
2344
     32,