C++博客-The Welkin of FlyingBear

个人网站进行中...

FlyingBear — Mon, 11 Jun 2007 11:21:00 GMT

看了一些很酷的个人网站,像曾经的燕尘小筑,感觉以后咱做这一行的,肯定是需要这样一个平台来展示自己的,做这个是早晚的事,试了一些第三方的BLOG,比如这个,总觉得不是自己的东西,用起来就是不爽,所以还是得自己弄一个,不过因为是代表个人的终极网站了,所以也不急,现在完成度大概只有1% :)

另外,俞益州老师来我们学校讲学的时候,我看到他也没做什么准备,电脑上有一个demo文件夹,里面放了各种PPT和演示用的视频,随手拿来就讲,这个很好,也要学,还有什么个人简历,创业计划什么的,能写好,整理好的,都放在一起,什么时候见到个志同道合的,立马拿出来介绍.哈哈

个人网站的设计目标,首先必须要HTML与PHP分离,这可以方便UI设计师的工作,不用去看超级ugly的PHP代码,这可以靠smarty模板引擎做到,小试了一下,用起来还是很直观的:)

另外,最好还需要PHP与逻辑层分离,这个没有仔细研究过,应该是要用一些框架,像Zend Framework,CakePHP什么的,口碑还不错.突然看到两张图,讲怎么开发好的PHP项目的,看过以后豁然开朗,share一下.
(点击看大图)

FlyingBear 2007-06-11 19:21 发表评论

SSE优化一例

FlyingBear — Fri, 04 May 2007 09:17:00 GMT

摘要: 前两天为了加速一段求梯度的代码，用了SSE指令，在实验室PMH大侠的指导下，最终实现了3倍速度提升（极限是4倍，因为4个浮点数一起计算）。在这里写一下心得，欢迎拍砖。SSE加速的几个关键是(1) 用于并行计算的数据结构要16字节对齐(2) 直接写汇编，不要用SSE的Load Store指令(3) 对于SSE本身不提供的三角函数等指令，可以用查表法，但要用SSE来算索引号相比起用GPU加速来说，SS... 阅读全文

FlyingBear 2007-05-04 17:17 发表评论

发现自己土了

FlyingBear — Sun, 29 Apr 2007 14:21:00 GMT

摘要: 写完上一篇日志"静态回调函数到成员函数的转换", 问了一下石教主有什么看法, 结果不问不知道, 一问才发现自己彻底土了, 这个问题用boost::bind可以解决得这么好, 调用的函数参数个数都可以没有限制... 唉, 以后还是要好好学习ps. 不过不知道这个方法运行效率怎么样, 还要继续试验// 下面的内容于2007年5月4日改成最后用的版本了实现的功能:(1) 用于把非多线程的代码改成多线程(... 阅读全文

FlyingBear 2007-04-29 22:21 发表评论

静态回调函数到成员函数的转换

FlyingBear — Sat, 21 Apr 2007 06:23:00 GMT

经常会有一些类库或者API要求传入一个等效于全局函数的函数指针作为回调函数,一个典型的例子是Win32的建立线程

DWORD WINAPI ThreadFunc(LPVOID lpParameter);
HANDLE hThread = CreateThread(NULL, NULL, ThreadFunc, NULL, NULL, NULL);

然而,对于我们自己的工程来说,更希望作为线程的函数是某个类的成员函数,所以需要在这个全局函数里调用类的成员函数,像这样

ClassA a;
DWORD WINAPI ThreadFunc(LPVOID lpParameter)
{
return ((ClassA*) lpParameter)->SomeMethod();
}

HANDLE hThread = CreateThread(NULL, NULL, ThreadFunc, &a, NULL, NULL);

DirectX SDK以前就用过这样的方法,也可以把全局函数换成类的静态成员函数,

但这样用起来很麻烦,每次都要写一个全局函数的Shell,所以这里,我们想办法把这个过程自动化起来.

问题的关键是在静态或者全局的函数里,只能得到一个传进来的参数,而要指定一个成员函数,需要一个this指针,以及一个指向类成员函数的指针,这两个参数没法全部通过lpParameter传进来,除非lpParameter传一个额外写的结构的指针,这个结构里包含this指针和成员函数指针,但这样需要临时分配一个对象,不方便,所以只能牺牲运行时代码的简洁性,用模板参数来传成员函数指针的值,而lpParameter只负责传this,但是另一个问题是用模板参数来传具体的值时,必须类型已知,像下面的代码是不能运作的

template<typename T, PtrToMemThreadFun pFunc>
DWORD WINAPI ThreadFunc(LPVOID lpParameter)
{
return (((T*)lpParameter)->*pFunc)();
}

CreateThread(NULL, NULL, ThreadFunc<ClassA, &ClassA::SomeMethod>, &a, NULL, NULL);

因为PtrToMemThreadFun的类型不定,所以只能先外包一层模板类来确定类型,像这样

template<typename T>
struct ThreadFac
{
    typedef DWORD (T::*PtrToMemThreadFunc)();
    typedef T ClassType;

    template<PtrToMemThreadFunc pFunc>
    static DWORD WINAPI ThreadFunc(LPVOID lpParameter)
    {
        return (((ClassType*)lpParameter)->*pFunc)();
    }
};

调用时,就可以自动化了

// Globally
CreateThread(NULL, NULL, ThreadFac<ClassA>::ThreadFunc<&ClassA::SomeMethod>, &a, NULL, NULL);

// In methods of ClassA
CreateThread(NULL, NULL, ThreadFac<ClassA>::ThreadFunc<&ClassA::SomeMethod>, this, NULL, NULL);

FlyingBear 2007-04-21 14:23 发表评论

日志系统的几则妙用

FlyingBear — Sat, 17 Feb 2007 14:19:00 GMT

1.HTML日志
对日志系统的文本输出部分稍做修改，每一行前加

，后加

，输出结果就成了一个可浏览的HTML，然后，凡HTML能包括的东西，都能让日志系统去LOG，做成一个图文并茂的日志，这对于科研应用来说有很重要的意义，因为经常生成临时的图像，或者需要对某个数据结果可视化，这时候如果只要一个简单的语句如

LogManager::getSingleton().logImage(someImage)

或者用宏进一步简化

LOG(someImage)

就完成了临时图像的自动命名，并保存，然后添加到日志里，那将是很有吸引力的，生成的日志有点像一个Report，可以直接发给Boss看，当工作汇报 :)

为了对LOG进去的媒体内容有一定的处理，如显示的时候缩放到合适的大小等，需要在日志文件的头部加入一些Javascript的代码，这需要日志系统的支持，能在开始记日志的时候判断当前日志文件是否已经存在，只在不存在的时候才加入相应的页首代码。

当然，更进一步可以做成结构化的日志，这个还有待尝试

2.Comment - Report 机制
如果写的程序经常发生变化，那么时间长了以后，日志下来的数据是由什么样的代码产生的，就会不记得了，导致大量的数据没有办法分析，这里我用这样的办法

在产生数据的代码附近用这样一个语句

COMMENT(subject)

subject是一个字符串值，表示这一部分数据的主题，这个宏展开以后是这样的

LogManager::getSingleton().addComment(subject, __FILE__, __LINE__);

这个函数在HTML的日志文件里新增一行，并生成一个超链接目的地，同时LogManager记下这个(string, string, int)组合

在LogManager被关闭的时候，会对所有的 Comment根据 Subject归类，生成一个Techinical Report，像下面这样

Tech Report:
---------------------------------------------------------
* SubjectA:
Comment1 - FileName - LineNo
Code Segment of Comment1
Comment2 - FileName - LineNo
Code Segment of Comment2
...
* SubjectB:
...
---------------------------------------------------------
同一个Subject可以对应多个代码段，单击Report里的 Comment 可以跳转到先前 Comment生成的超链接日的地。
Code Segment由LogManager根据__FILE__和__LINE__从文件里读取

3.变量记录宏

一些有用的宏对记录变量很有好处
#define LOGVAR(x) ToolLib::LOG(#x + std::string(" = ") + TOSTR(x))
#define LOGARR(x, i) ToolLib::LOG(#x + std::string("[") + TOSTR(i) + "] = " + TOSTR(x[i]))

用了这些宏就可以随处记录变量的值

4.直接数据记录

LogManager需要有一个方法给出所用的文件对象,如 ofstream或 FILE*，这样，对于大数据量的日志，比如几千个int值，可以用这个文件对象直接输入，而省去每一行记录日期时间的麻烦

5.替换stdout或cout

经常有些第三方的库用printf或者cout来输出信息，需要想办法把这些内容也catch到自己的日志系统里，对于printf模式的，可以用

FILE* fid = fopen();
*stdout = *fid;

来实现，而对于cout型的，不好办，可以自己写一个新的ostream子类，然后用宏把cout重新define一下，设成自己的ostream，并重写一下sync函数来capture输出。

FlyingBear 2007-02-17 22:19 发表评论

Boost::Serialization试用手记

FlyingBear — Thu, 15 Feb 2007 17:06:00 GMT

摘要: 对于核心数据相对比较集中的应用程序来说，serialization机制可以直接充当文档保存与打开功能的实现工具，这也是很多成熟的应用程序framework都提供serialization支持的原因但是个人认为，serialization最精彩的用处在于保存现场，比如在探索性的科研应用程序开发过程中，很可能一部分算法已经固定下来，其余的有待进一步探索，而确定下来的部分有可能十分time consum... 阅读全文

FlyingBear 2007-02-16 01:06 发表评论

一组简易有效的性能测试函数

FlyingBear — Thu, 15 Feb 2007 13:18:00 GMT

在写与性能相关的应用程序时,经常需要测代码的运行时间,听说Intel Compiler下面用VTune可以做得很好,可惜不会用,受Matlab下的Tic,Toc语句的启发,写了一些类似的函数,希望可以实现这样的功能

HBENCHMARK hid;
Tic(hid);
// some code
Toc(hid);
// show time

如果要测的代码是一个循环体的一部分,又希望测得总时间,像下面这样

HBENCHMARK hid;
Tik(hid);
while()
{
    //  some code
    Trk(hid);
    //  code to benchmark
    Tsk(hid);
    //  some code
}
Tok(hid);
//  show time

这里的Trk表示恢复计时,Tsk表示暂停计时,Tok里面包括了停止计时和显示时间.这种方式同样可以应用于跨过程的计时,只要有统一的Benchmark Handle, 就可以,实现的代码如下

    typedef LARGE_INTEGER HBENCHMARK;

    //////////////////////////////////////////////////////////////////////////
    // Single turn timer routine

    // begin
    void Tic(HBENCHMARK& node);
    // end, display timer.
    void Toc(HBENCHMARK& node, std::string prefix = "");

    //////////////////////////////////////////////////////////////////////////
    // Multiple turns timer routine

    // initialize timer
    void Tik(HBENCHMARK& node);
    // resume timer
    void Trk(HBENCHMARK& node);
    // stop timer
    void Tsk(HBENCHMARK& node);
    // end, display timer
    void Tok(HBENCHMARK& node, std::string prefix = "");

    //-----------------------------------------------------------------------
    void Tic(HBENCHMARK& node)
    {
        QueryPerformanceCounter(&node);
    }
    //-----------------------------------------------------------------------
    void Toc(HBENCHMARK& node, std::string prefix)
    {
        static LARGE_INTEGER endtime;
        QueryPerformanceCounter(&endtime);

        node.QuadPart = endtime.QuadPart - node.QuadPart;

        if(prefix.length() > 0)
            ToolLib::LogManager::getSingleton().logMessage(prefix);
        ToolLib::LogManager::getSingleton().logMessage("Benchmark: " +
            TOSTR(node.HighPart) + "," + TOSTR(node.LowPart));
    }
    //-----------------------------------------------------------------------
    void Tik(HBENCHMARK& node)
    {
        node.QuadPart = 0;
    }
    //-----------------------------------------------------------------------
    void Trk(HBENCHMARK& node)
    {
        static LARGE_INTEGER curTime;
        QueryPerformanceCounter(&curTime);
        node.QuadPart = curTime.QuadPart - node.QuadPart;
    }
    //-----------------------------------------------------------------------
    void Tsk(HBENCHMARK& node)
    {
        static HBENCHMARK endTime;
        QueryPerformanceCounter(&endTime);
        node.QuadPart = endTime.QuadPart - node.QuadPart;
    }
    //-----------------------------------------------------------------------
    void Tok(HBENCHMARK& node, std::string prefix)
    {
        if(prefix.length() > 0)
            ToolLib::LogManager::getSingleton().logMessage(prefix);
        ToolLib::LogManager::getSingleton().logMessage("Benchmark: " +
            TOSTR(node.HighPart) + "," + TOSTR(node.LowPart));
    }

这里使用的是WindowsAPI QueryPerformanceCounter函数,为的是测得精确的时钟周期数

FlyingBear 2007-02-15 21:18 发表评论