注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

c预编译#ifdef命令demo  

2013-01-09 23:34:50|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
#ifdef raw
#include <raw.h>
#else
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <unistd.h>
#endif
#include <math.h>

void begin(void);
float calc_fir(const float *a, float *b, float W, float last);

/* Globals: */
static int numiters = -1;

#ifndef raw
int main(int argc, char **argv)
{
  int option;

  while ((option = getopt(argc, argv, "i:")) != -1)
  {
    switch(option)
    {
    case 'i':
      numiters = atoi(optarg);
    }
  }

  begin();
  return 0;
}
#endif

#define DEPTH 128
#define LENGTH 12

void begin(void)
{
  int val = 0;
  const int limit = 10000;
  int i;
  float W[DEPTH], last[DEPTH];
  float f1[LENGTH], f2[LENGTH];
  float *a, *b, *c;
  
  for (i = 0; i < DEPTH; i++)
  {
    W[i] = (float)(2*i*i) / (float)(i+1);
    last[i] = 0;
  }
  a = f1;
  b = f2;
  
  /* Main loop: */
  while (numiters == -1 || numiters-- > 0)
  {
    for (i = 0; i < LENGTH; i += 2)
    {
      a[i] = 0;
      a[i+1] = val++;
      if (val >= limit) val = 0;
    }
    for (i = 0; i < DEPTH; i++)
    {
      last[i] = calc_fir(a, b, W[i], last[i]);
      c = a;
      a = b;
      b = c;
    }
      
    for (i = 0; i < LENGTH; i += 2)
#ifdef raw
      print_float(a[i]);
#else
      printf("%f\n", a[i]);
#endif
  }
}

float calc_fir(const float *a, float *b, float W, float last)
{
  int i;
  
  for (i = 0; i < LENGTH/2; i++)
  {
    float s = a[i*2];
    b[i*2] = s + last * W;
    b[i*2+1] = last;
    last = a[i*2+1];
  }
  return last;
}

  评论这张
 
阅读(444)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016