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

一粒浮尘

飘渺虚无

 
 
 

日志

 
 

Coding in night ...  

2011-11-28 03:40:46|  分类: work |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
夜已深,我却还在调试程序。。。
写程序可真是个细活。。。

前几天开始重新写一个计算转移函数的程序,主要部分是重新设计光深,重子声速等的计算,以前只是翻译了一下cmbfast,所以那算不上自己的程序。

以前翻译都花了我许多的时间,从不懂fortran,到把cmbfast搞透,很不易啊。。。再翻译成c版本的(其实是提取出自己需要的计算部分)也不容易,因为必须要对计算的细节了如指掌,我想我现在可以做到了。

以前写成功第一个可以准确计算的程序,花了好几个月,当然时间不仅仅浩在写程序本身,更多的其实是用在了先读懂cmbfast,这对一个从未做过数值的我来说真的很难很难,当初真的是经常处于崩溃的边缘,苦苦挣扎了那么久,我终于可以自己独立写出一些像样的程序了。调试程序的经验也越来越丰富,可是今天要写的经历还是给我一次提醒,要细心细心再细心。。。因为我这两天耗费了我大部分时间的,仅仅是几个很不起眼的错误,有两个除号不小心弄成乘号了,对比新旧版本的时候也没看出来,型号还是发现了,不然明天要做更大规模的调试。。。另外一个小错误是少了一个系数,虽不会造成大的误差,但毕竟是存在一定的危险的。

所以以后一定要更加小心~

对比了这次编程的过程,比以前成熟多了!就几天的功夫,写出了两个版本,第一个版本就存在上面提到的错误,但是没检查出来。。。丢人阿。。。然后我以为是GSL的插值程序的计算效率不高,于是就把用GSL插值函数的部分大幅削减,可是程序运行依旧不怎么高效,知道真正的原因被发现,呵呵

时间差不多了,该去睡觉了~
明天可以继续添加新的东西了~

  评论这张
 
阅读(141)| 评论(2)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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