CalculiX  2.13
A Free Software Three-Dimensional Structural Finite Element Program
Functions/Subroutines
drfftf.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine radf2 (IDO, L1, CC, CH, WA1)
 
subroutine radf3 (IDO, L1, CC, CH, WA1, WA2)
 
subroutine radf4 (IDO, L1, CC, CH, WA1, WA2, WA3)
 
subroutine radf5 (IDO, L1, CC, CH, WA1, WA2, WA3, WA4)
 
subroutine radfg (IDO, IP, L1, IDL1, CC, C1, C2, CH, CH2, WA)
 
subroutine drfftf (N, R, WSAVE, isave)
 
subroutine drfftf1 (N, C, CH, WA, IFAC)
 
subroutine drffti (N, WSAVE, isave)
 
subroutine drffti1 (N, WA, IFAC)
 

Function/Subroutine Documentation

◆ drfftf()

subroutine drfftf (   N,
dimension(1)  R,
dimension(*)  WSAVE,
dimension(*)  isave 
)
467  implicit real*8(a-h,o-z)
468  dimension r(1) ,wsave(*),isave(*)
469  IF (n .EQ. 1) RETURN
470  CALL drfftf1 (n,r,wsave,wsave(n+1),isave)
471  RETURN
drfftf1
subroutine drfftf1(N, C, CH, WA, IFAC)
Definition: drfftf.f:474

◆ drfftf1()

subroutine drfftf1 (   N,
dimension(1)  C,
dimension(1)  CH,
dimension(1)  WA,
dimension(*)  IFAC 
)
474  implicit real*8(a-h,o-z)
475  dimension ch(1) ,c(1) ,wa(1) ,ifac(*)
476  nf = ifac(2)
477  na = 1
478  l2 = n
479  iw = n
480  DO 111 k1=1,nf
481  kh = nf-k1
482  ip = ifac(kh+3)
483  l1 = l2/ip
484  ido = n/l2
485  idl1 = ido*l1
486  iw = iw-(ip-1)*ido
487  na = 1-na
488  IF (ip .NE. 4) GO TO 102
489  ix2 = iw+ido
490  ix3 = ix2+ido
491  IF (na .NE. 0) GO TO 101
492  CALL radf4 (ido,l1,c,ch,wa(iw),wa(ix2),wa(ix3))
493  GO TO 110
494  101 CALL radf4 (ido,l1,ch,c,wa(iw),wa(ix2),wa(ix3))
495  GO TO 110
496  102 IF (ip .NE. 2) GO TO 104
497  IF (na .NE. 0) GO TO 103
498  CALL radf2 (ido,l1,c,ch,wa(iw))
499  GO TO 110
500  103 CALL radf2 (ido,l1,ch,c,wa(iw))
501  GO TO 110
502  104 IF (ip .NE. 3) GO TO 106
503  ix2 = iw+ido
504  IF (na .NE. 0) GO TO 105
505  CALL radf3 (ido,l1,c,ch,wa(iw),wa(ix2))
506  GO TO 110
507  105 CALL radf3 (ido,l1,ch,c,wa(iw),wa(ix2))
508  GO TO 110
509  106 IF (ip .NE. 5) GO TO 108
510  ix2 = iw+ido
511  ix3 = ix2+ido
512  ix4 = ix3+ido
513  IF (na .NE. 0) GO TO 107
514  CALL radf5 (ido,l1,c,ch,wa(iw),wa(ix2),wa(ix3),wa(ix4))
515  GO TO 110
516  107 CALL radf5 (ido,l1,ch,c,wa(iw),wa(ix2),wa(ix3),wa(ix4))
517  GO TO 110
518  108 IF (ido .EQ. 1) na = 1-na
519  IF (na .NE. 0) GO TO 109
520  CALL radfg (ido,ip,l1,idl1,c,c,c,ch,ch,wa(iw))
521  na = 1
522  GO TO 110
523  109 CALL radfg (ido,ip,l1,idl1,ch,ch,ch,c,c,wa(iw))
524  na = 0
525  110 l2 = l1
526  111 CONTINUE
527  IF (na .EQ. 1) RETURN
528  DO 112 i=1,n
529  c(i) = ch(i)
530  112 CONTINUE
531  RETURN
radf5
subroutine radf5(IDO, L1, CC, CH, WA1, WA2, WA3, WA4)
Definition: drfftf.f:240
radf2
subroutine radf2(IDO, L1, CC, CH, WA1)
Definition: drfftf.f:122
radf3
subroutine radf3(IDO, L1, CC, CH, WA1, WA2)
Definition: drfftf.f:150
radfg
subroutine radfg(IDO, IP, L1, IDL1, CC, C1, C2, CH, CH2, WA)
Definition: drfftf.f:300
radf4
subroutine radf4(IDO, L1, CC, CH, WA1, WA2, WA3)
Definition: drfftf.f:186

◆ drffti()

subroutine drffti (   N,
dimension(*)  WSAVE,
dimension(*)  isave 
)
534  implicit real*8(a-h,o-z)
535  dimension wsave(*),isave(*)
536  IF (n .EQ. 1) RETURN
537  CALL drffti1 (n,wsave(n+1),isave)
538  RETURN
drffti1
subroutine drffti1(N, WA, IFAC)
Definition: drfftf.f:541

◆ drffti1()

subroutine drffti1 (   N,
dimension(1)  WA,
dimension(*)  IFAC 
)
541  implicit real*8(a-h,o-z)
542  dimension wa(1) ,ifac(*) ,ntryh(4)
543  DATA ntryh(1),ntryh(2),ntryh(3),ntryh(4)/4,2,3,5/
544  nl = n
545  nf = 0
546  j = 0
547  101 j = j+1
548  IF (j-4) 102,102,103
549  102 ntry = ntryh(j)
550  GO TO 104
551  103 ntry = ntry+2
552  104 nq = nl/ntry
553  nr = nl-ntry*nq
554  IF (nr) 101,105,101
555  105 nf = nf+1
556  ifac(nf+2) = ntry
557  nl = nq
558  IF (ntry .NE. 2) GO TO 107
559  IF (nf .EQ. 1) GO TO 107
560  DO 106 i=2,nf
561  ib = nf-i+2
562  ifac(ib+2) = ifac(ib+1)
563  106 CONTINUE
564  ifac(3) = 2
565  107 IF (nl .NE. 1) GO TO 104
566  ifac(1) = n
567  ifac(2) = nf
568  tpi = 6.28318530717959
569  argh = tpi/float(n)
570  is = 0
571  nfm1 = nf-1
572  l1 = 1
573  IF (nfm1 .EQ. 0) RETURN
574  DO 110 k1=1,nfm1
575  ip = ifac(k1+2)
576  ld = 0
577  l2 = l1*ip
578  ido = n/l2
579  ipm = ip-1
580  DO 109 j=1,ipm
581  ld = ld+l1
582  i = is
583  argld = float(ld)*argh
584  fi = 0.
585  DO 108 ii=3,ido,2
586  i = i+2
587  fi = fi+1.
588  arg = fi*argld
589  wa(i-1) = cos(arg)
590  wa(i) = sin(arg)
591  108 CONTINUE
592  is = is+ido
593  109 CONTINUE
594  l1 = l2
595  110 CONTINUE
596  RETURN

◆ radf2()

subroutine radf2 (   IDO,
  L1,
dimension(ido,l1,2)  CC,
dimension(ido,2,l1)  CH,
dimension(1)  WA1 
)
122  implicit real*8(a-h,o-z)
123  dimension ch(ido,2,l1) ,cc(ido,l1,2) ,
124  1 wa1(1)
125  DO 101 k=1,l1
126  ch(1,1,k) = cc(1,k,1)+cc(1,k,2)
127  ch(ido,2,k) = cc(1,k,1)-cc(1,k,2)
128  101 CONTINUE
129  IF (ido-2) 107,105,102
130  102 idp2 = ido+2
131  DO 104 k=1,l1
132  DO 103 i=3,ido,2
133  ic = idp2-i
134  tr2 = wa1(i-2)*cc(i-1,k,2)+wa1(i-1)*cc(i,k,2)
135  ti2 = wa1(i-2)*cc(i,k,2)-wa1(i-1)*cc(i-1,k,2)
136  ch(i,1,k) = cc(i,k,1)+ti2
137  ch(ic,2,k) = ti2-cc(i,k,1)
138  ch(i-1,1,k) = cc(i-1,k,1)+tr2
139  ch(ic-1,2,k) = cc(i-1,k,1)-tr2
140  103 CONTINUE
141  104 CONTINUE
142  IF (mod(ido,2) .EQ. 1) RETURN
143  105 DO 106 k=1,l1
144  ch(1,2,k) = -cc(ido,k,2)
145  ch(ido,1,k) = cc(ido,k,1)
146  106 CONTINUE
147  107 RETURN

◆ radf3()

subroutine radf3 (   IDO,
  L1,
dimension(ido,l1,3)  CC,
dimension(ido,3,l1)  CH,
dimension(1)  WA1,
dimension(1)  WA2 
)
150  implicit real*8(a-h,o-z)
151  dimension ch(ido,3,l1) ,cc(ido,l1,3) ,
152  1 wa1(1) ,wa2(1)
153  DATA taur,taui /-.5,.866025403784439/
154  DO 101 k=1,l1
155  cr2 = cc(1,k,2)+cc(1,k,3)
156  ch(1,1,k) = cc(1,k,1)+cr2
157  ch(1,3,k) = taui*(cc(1,k,3)-cc(1,k,2))
158  ch(ido,2,k) = cc(1,k,1)+taur*cr2
159  101 CONTINUE
160  IF (ido .EQ. 1) RETURN
161  idp2 = ido+2
162  DO 103 k=1,l1
163  DO 102 i=3,ido,2
164  ic = idp2-i
165  dr2 = wa1(i-2)*cc(i-1,k,2)+wa1(i-1)*cc(i,k,2)
166  di2 = wa1(i-2)*cc(i,k,2)-wa1(i-1)*cc(i-1,k,2)
167  dr3 = wa2(i-2)*cc(i-1,k,3)+wa2(i-1)*cc(i,k,3)
168  di3 = wa2(i-2)*cc(i,k,3)-wa2(i-1)*cc(i-1,k,3)
169  cr2 = dr2+dr3
170  ci2 = di2+di3
171  ch(i-1,1,k) = cc(i-1,k,1)+cr2
172  ch(i,1,k) = cc(i,k,1)+ci2
173  tr2 = cc(i-1,k,1)+taur*cr2
174  ti2 = cc(i,k,1)+taur*ci2
175  tr3 = taui*(di2-di3)
176  ti3 = taui*(dr3-dr2)
177  ch(i-1,3,k) = tr2+tr3
178  ch(ic-1,2,k) = tr2-tr3
179  ch(i,3,k) = ti2+ti3
180  ch(ic,2,k) = ti3-ti2
181  102 CONTINUE
182  103 CONTINUE
183  RETURN

◆ radf4()

subroutine radf4 (   IDO,
  L1,
dimension(ido,l1,4)  CC,
dimension(ido,4,l1)  CH,
dimension(1)  WA1,
dimension(1)  WA2,
dimension(1)  WA3 
)
186  implicit real*8(a-h,o-z)
187  dimension cc(ido,l1,4) ,ch(ido,4,l1) ,
188  1 wa1(1) ,wa2(1) ,wa3(1)
189  DATA hsqt2 /.7071067811865475/
190  DO 101 k=1,l1
191  tr1 = cc(1,k,2)+cc(1,k,4)
192  tr2 = cc(1,k,1)+cc(1,k,3)
193  ch(1,1,k) = tr1+tr2
194  ch(ido,4,k) = tr2-tr1
195  ch(ido,2,k) = cc(1,k,1)-cc(1,k,3)
196  ch(1,3,k) = cc(1,k,4)-cc(1,k,2)
197  101 CONTINUE
198  IF (ido-2) 107,105,102
199  102 idp2 = ido+2
200  DO 104 k=1,l1
201  DO 103 i=3,ido,2
202  ic = idp2-i
203  cr2 = wa1(i-2)*cc(i-1,k,2)+wa1(i-1)*cc(i,k,2)
204  ci2 = wa1(i-2)*cc(i,k,2)-wa1(i-1)*cc(i-1,k,2)
205  cr3 = wa2(i-2)*cc(i-1,k,3)+wa2(i-1)*cc(i,k,3)
206  ci3 = wa2(i-2)*cc(i,k,3)-wa2(i-1)*cc(i-1,k,3)
207  cr4 = wa3(i-2)*cc(i-1,k,4)+wa3(i-1)*cc(i,k,4)
208  ci4 = wa3(i-2)*cc(i,k,4)-wa3(i-1)*cc(i-1,k,4)
209  tr1 = cr2+cr4
210  tr4 = cr4-cr2
211  ti1 = ci2+ci4
212  ti4 = ci2-ci4
213  ti2 = cc(i,k,1)+ci3
214  ti3 = cc(i,k,1)-ci3
215  tr2 = cc(i-1,k,1)+cr3
216  tr3 = cc(i-1,k,1)-cr3
217  ch(i-1,1,k) = tr1+tr2
218  ch(ic-1,4,k) = tr2-tr1
219  ch(i,1,k) = ti1+ti2
220  ch(ic,4,k) = ti1-ti2
221  ch(i-1,3,k) = ti4+tr3
222  ch(ic-1,2,k) = tr3-ti4
223  ch(i,3,k) = tr4+ti3
224  ch(ic,2,k) = tr4-ti3
225  103 CONTINUE
226  104 CONTINUE
227  IF (mod(ido,2) .EQ. 1) RETURN
228  105 CONTINUE
229  DO 106 k=1,l1
230  ti1 = -hsqt2*(cc(ido,k,2)+cc(ido,k,4))
231  tr1 = hsqt2*(cc(ido,k,2)-cc(ido,k,4))
232  ch(ido,1,k) = tr1+cc(ido,k,1)
233  ch(ido,3,k) = cc(ido,k,1)-tr1
234  ch(1,2,k) = ti1-cc(ido,k,3)
235  ch(1,4,k) = ti1+cc(ido,k,3)
236  106 CONTINUE
237  107 RETURN

◆ radf5()

subroutine radf5 (   IDO,
  L1,
dimension(ido,l1,5)  CC,
dimension(ido,5,l1)  CH,
dimension(1)  WA1,
dimension(1)  WA2,
dimension(1)  WA3,
dimension(1)  WA4 
)
240  implicit real*8(a-h,o-z)
241  dimension cc(ido,l1,5) ,ch(ido,5,l1) ,
242  1 wa1(1) ,wa2(1) ,wa3(1) ,wa4(1)
243  DATA tr11,ti11,tr12,ti12 /.309016994374947,.951056516295154,
244  1-.809016994374947,.587785252292473/
245  DO 101 k=1,l1
246  cr2 = cc(1,k,5)+cc(1,k,2)
247  ci5 = cc(1,k,5)-cc(1,k,2)
248  cr3 = cc(1,k,4)+cc(1,k,3)
249  ci4 = cc(1,k,4)-cc(1,k,3)
250  ch(1,1,k) = cc(1,k,1)+cr2+cr3
251  ch(ido,2,k) = cc(1,k,1)+tr11*cr2+tr12*cr3
252  ch(1,3,k) = ti11*ci5+ti12*ci4
253  ch(ido,4,k) = cc(1,k,1)+tr12*cr2+tr11*cr3
254  ch(1,5,k) = ti12*ci5-ti11*ci4
255  101 CONTINUE
256  IF (ido .EQ. 1) RETURN
257  idp2 = ido+2
258  DO 103 k=1,l1
259  DO 102 i=3,ido,2
260  ic = idp2-i
261  dr2 = wa1(i-2)*cc(i-1,k,2)+wa1(i-1)*cc(i,k,2)
262  di2 = wa1(i-2)*cc(i,k,2)-wa1(i-1)*cc(i-1,k,2)
263  dr3 = wa2(i-2)*cc(i-1,k,3)+wa2(i-1)*cc(i,k,3)
264  di3 = wa2(i-2)*cc(i,k,3)-wa2(i-1)*cc(i-1,k,3)
265  dr4 = wa3(i-2)*cc(i-1,k,4)+wa3(i-1)*cc(i,k,4)
266  di4 = wa3(i-2)*cc(i,k,4)-wa3(i-1)*cc(i-1,k,4)
267  dr5 = wa4(i-2)*cc(i-1,k,5)+wa4(i-1)*cc(i,k,5)
268  di5 = wa4(i-2)*cc(i,k,5)-wa4(i-1)*cc(i-1,k,5)
269  cr2 = dr2+dr5
270  ci5 = dr5-dr2
271  cr5 = di2-di5
272  ci2 = di2+di5
273  cr3 = dr3+dr4
274  ci4 = dr4-dr3
275  cr4 = di3-di4
276  ci3 = di3+di4
277  ch(i-1,1,k) = cc(i-1,k,1)+cr2+cr3
278  ch(i,1,k) = cc(i,k,1)+ci2+ci3
279  tr2 = cc(i-1,k,1)+tr11*cr2+tr12*cr3
280  ti2 = cc(i,k,1)+tr11*ci2+tr12*ci3
281  tr3 = cc(i-1,k,1)+tr12*cr2+tr11*cr3
282  ti3 = cc(i,k,1)+tr12*ci2+tr11*ci3
283  tr5 = ti11*cr5+ti12*cr4
284  ti5 = ti11*ci5+ti12*ci4
285  tr4 = ti12*cr5-ti11*cr4
286  ti4 = ti12*ci5-ti11*ci4
287  ch(i-1,3,k) = tr2+tr5
288  ch(ic-1,2,k) = tr2-tr5
289  ch(i,3,k) = ti2+ti5
290  ch(ic,2,k) = ti5-ti2
291  ch(i-1,5,k) = tr3+tr4
292  ch(ic-1,4,k) = tr3-tr4
293  ch(i,5,k) = ti3+ti4
294  ch(ic,4,k) = ti4-ti3
295  102 CONTINUE
296  103 CONTINUE
297  RETURN

◆ radfg()

subroutine radfg (   IDO,
  IP,
  L1,
  IDL1,
dimension(ido,ip,l1)  CC,
dimension(ido,l1,ip)  C1,
dimension(idl1,ip)  C2,
dimension(ido,l1,ip)  CH,
dimension(idl1,ip)  CH2,
dimension(1)  WA 
)
300  implicit real*8(a-h,o-z)
301  dimension ch(ido,l1,ip) ,cc(ido,ip,l1) ,
302  1 c1(ido,l1,ip) ,c2(idl1,ip),
303  2 ch2(idl1,ip) ,wa(1)
304  DATA tpi/6.28318530717959/
305  arg = tpi/float(ip)
306  dcp = cos(arg)
307  dsp = sin(arg)
308  ipph = (ip+1)/2
309  ipp2 = ip+2
310  idp2 = ido+2
311  nbd = (ido-1)/2
312  IF (ido .EQ. 1) GO TO 119
313  DO 101 ik=1,idl1
314  ch2(ik,1) = c2(ik,1)
315  101 CONTINUE
316  DO 103 j=2,ip
317  DO 102 k=1,l1
318  ch(1,k,j) = c1(1,k,j)
319  102 CONTINUE
320  103 CONTINUE
321  IF (nbd .GT. l1) GO TO 107
322  is = -ido
323  DO 106 j=2,ip
324  is = is+ido
325  idij = is
326  DO 105 i=3,ido,2
327  idij = idij+2
328  DO 104 k=1,l1
329  ch(i-1,k,j) = wa(idij-1)*c1(i-1,k,j)+wa(idij)*c1(i,k,j)
330  ch(i,k,j) = wa(idij-1)*c1(i,k,j)-wa(idij)*c1(i-1,k,j)
331  104 CONTINUE
332  105 CONTINUE
333  106 CONTINUE
334  GO TO 111
335  107 is = -ido
336  DO 110 j=2,ip
337  is = is+ido
338  DO 109 k=1,l1
339  idij = is
340  DO 108 i=3,ido,2
341  idij = idij+2
342  ch(i-1,k,j) = wa(idij-1)*c1(i-1,k,j)+wa(idij)*c1(i,k,j)
343  ch(i,k,j) = wa(idij-1)*c1(i,k,j)-wa(idij)*c1(i-1,k,j)
344  108 CONTINUE
345  109 CONTINUE
346  110 CONTINUE
347  111 IF (nbd .LT. l1) GO TO 115
348  DO 114 j=2,ipph
349  jc = ipp2-j
350  DO 113 k=1,l1
351  DO 112 i=3,ido,2
352  c1(i-1,k,j) = ch(i-1,k,j)+ch(i-1,k,jc)
353  c1(i-1,k,jc) = ch(i,k,j)-ch(i,k,jc)
354  c1(i,k,j) = ch(i,k,j)+ch(i,k,jc)
355  c1(i,k,jc) = ch(i-1,k,jc)-ch(i-1,k,j)
356  112 CONTINUE
357  113 CONTINUE
358  114 CONTINUE
359  GO TO 121
360  115 DO 118 j=2,ipph
361  jc = ipp2-j
362  DO 117 i=3,ido,2
363  DO 116 k=1,l1
364  c1(i-1,k,j) = ch(i-1,k,j)+ch(i-1,k,jc)
365  c1(i-1,k,jc) = ch(i,k,j)-ch(i,k,jc)
366  c1(i,k,j) = ch(i,k,j)+ch(i,k,jc)
367  c1(i,k,jc) = ch(i-1,k,jc)-ch(i-1,k,j)
368  116 CONTINUE
369  117 CONTINUE
370  118 CONTINUE
371  GO TO 121
372  119 DO 120 ik=1,idl1
373  c2(ik,1) = ch2(ik,1)
374  120 CONTINUE
375  121 DO 123 j=2,ipph
376  jc = ipp2-j
377  DO 122 k=1,l1
378  c1(1,k,j) = ch(1,k,j)+ch(1,k,jc)
379  c1(1,k,jc) = ch(1,k,jc)-ch(1,k,j)
380  122 CONTINUE
381  123 CONTINUE
382 C
383  ar1 = 1.
384  ai1 = 0.
385  DO 127 l=2,ipph
386  lc = ipp2-l
387  ar1h = dcp*ar1-dsp*ai1
388  ai1 = dcp*ai1+dsp*ar1
389  ar1 = ar1h
390  DO 124 ik=1,idl1
391  ch2(ik,l) = c2(ik,1)+ar1*c2(ik,2)
392  ch2(ik,lc) = ai1*c2(ik,ip)
393  124 CONTINUE
394  dc2 = ar1
395  ds2 = ai1
396  ar2 = ar1
397  ai2 = ai1
398  DO 126 j=3,ipph
399  jc = ipp2-j
400  ar2h = dc2*ar2-ds2*ai2
401  ai2 = dc2*ai2+ds2*ar2
402  ar2 = ar2h
403  DO 125 ik=1,idl1
404  ch2(ik,l) = ch2(ik,l)+ar2*c2(ik,j)
405  ch2(ik,lc) = ch2(ik,lc)+ai2*c2(ik,jc)
406  125 CONTINUE
407  126 CONTINUE
408  127 CONTINUE
409  DO 129 j=2,ipph
410  DO 128 ik=1,idl1
411  ch2(ik,1) = ch2(ik,1)+c2(ik,j)
412  128 CONTINUE
413  129 CONTINUE
414 C
415  IF (ido .LT. l1) GO TO 132
416  DO 131 k=1,l1
417  DO 130 i=1,ido
418  cc(i,1,k) = ch(i,k,1)
419  130 CONTINUE
420  131 CONTINUE
421  GO TO 135
422  132 DO 134 i=1,ido
423  DO 133 k=1,l1
424  cc(i,1,k) = ch(i,k,1)
425  133 CONTINUE
426  134 CONTINUE
427  135 DO 137 j=2,ipph
428  jc = ipp2-j
429  j2 = j+j
430  DO 136 k=1,l1
431  cc(ido,j2-2,k) = ch(1,k,j)
432  cc(1,j2-1,k) = ch(1,k,jc)
433  136 CONTINUE
434  137 CONTINUE
435  IF (ido .EQ. 1) RETURN
436  IF (nbd .LT. l1) GO TO 141
437  DO 140 j=2,ipph
438  jc = ipp2-j
439  j2 = j+j
440  DO 139 k=1,l1
441  DO 138 i=3,ido,2
442  ic = idp2-i
443  cc(i-1,j2-1,k) = ch(i-1,k,j)+ch(i-1,k,jc)
444  cc(ic-1,j2-2,k) = ch(i-1,k,j)-ch(i-1,k,jc)
445  cc(i,j2-1,k) = ch(i,k,j)+ch(i,k,jc)
446  cc(ic,j2-2,k) = ch(i,k,jc)-ch(i,k,j)
447  138 CONTINUE
448  139 CONTINUE
449  140 CONTINUE
450  RETURN
451  141 DO 144 j=2,ipph
452  jc = ipp2-j
453  j2 = j+j
454  DO 143 i=3,ido,2
455  ic = idp2-i
456  DO 142 k=1,l1
457  cc(i-1,j2-1,k) = ch(i-1,k,j)+ch(i-1,k,jc)
458  cc(ic-1,j2-2,k) = ch(i-1,k,j)-ch(i-1,k,jc)
459  cc(i,j2-1,k) = ch(i,k,j)+ch(i,k,jc)
460  cc(ic,j2-2,k) = ch(i,k,jc)-ch(i,k,j)
461  142 CONTINUE
462  143 CONTINUE
463  144 CONTINUE
464  RETURN
c1
static double * c1
Definition: mafillvcompmain.c:30
Hosted by OpenAircraft.com, (Michigan UAV, LLC)