C++ 调用eigen作矩阵运算

C++ 调用eigen做矩阵运算
  我用C++调用eigen做矩阵运算时候结果出现2.7E+308的错误
void CGet_DCB::Cal_Part_Nbb(MatrixXd &nbb1_,MatrixXd &nbb2_,MatrixXd &nbb3_,MatrixXd &w1_,MatrixXd &w2_)
{
int start_ = num_obs[period-1],end_ = num_obs[period];
MatrixXd b1(end_ - start_,num_site+32);;
MatrixXd b2(end_ - start_,25);
MatrixXd l(end_ - start_,1);
for (int i = start_;i < end_;i++)
{
l(i-start_,0) = l_mat[i];
vector<double> row = b_mat[i];
for (int j = 0;j< num_site + 32;j++)
b1(i-start_,j) = row[j];
for (int k = 0;k < 25; k++)
b2(i-start_,k) = row[num_site+32+k];
}
//nbb2_ = b1.transpose()*b2;
nbb1_ = b1.transpose()*b1;
nbb3_ = b2.transpose()*b2;
w1_ = b1.transpose()*l;
w2_ = b2.transpose()*l;
}
b1、b2、l都没错
求得的nbb1_ 、w1_ 也是正常的nbb2_ 、nbb3_ 、w2_ 部分结果错误，请问这是为什么？
------解决方案--------------------
单步调试和设断点调试是程序员必须掌握的技能之一。

在觉得可能出问题的地方添加输出相关变量的值到日志中。

写日志用这个：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#ifdef WIN32
    #include <windows.h>
    #include <io.h>
#else
    #include <unistd.h>
    #include <sys/time.h>
    #include <pthread.h>
    #define  CRITICAL_SECTION   pthread_mutex_t
    #define  _vsnprintf         vsnprintf
#endif
//Log{
#define MAXLOGSIZE 20000000
#define MAXLINSIZE 16000
#include <time.h>
#include <sys/timeb.h>
#include <stdarg.h>
char logfilename1[]="MyLog1.log";
char logfilename2[]="MyLog2.log";
static char logstr[MAXLINSIZE+1];
char datestr[16];
char timestr[16];
char mss[4];
CRITICAL_SECTION cs_log;
FILE *flog;
#ifdef WIN32
void Lock(CRITICAL_SECTION *l) {
    EnterCriticalSection(l);
}
void Unlock(CRITICAL_SECTION *l) {
    LeaveCriticalSection(l);
}
#else
void Lock(CRITICAL_SECTION *l) {
    pthread_mutex_lock(l);
}
void Unlock(CRITICAL_SECTION *l) {
    pthread_mutex_unlock(l);
}
#endif
void LogV(const char *pszFmt,va_list argp) {
    struct tm *now;
    struct timeb tb;

    if (NULL==pszFmt
------解决方案--------------------
0==pszFmt[0]) return;
    _vsnprintf(logstr,MAXLINSIZE,pszFmt,argp);
    ftime(&tb);
    now=localtime(&tb.time);
    sprintf(datestr,"%04d-%02d-%02d",now->tm_year+1900,now->tm_mon+1,now->tm_mday);
    sprintf(timestr,"%02d:%02d:%02d",now->tm_hour     ,now->tm_min  ,now->tm_sec );
    sprintf(mss,"%03d",tb.millitm);
    printf("%s %s.%s %s",datestr,timestr,mss,logstr);
    flog=fopen(logfilename1,"a");
    if (NULL!=flog) {
        fprintf(flog,"%s %s.%s %s",datestr,timestr,mss,logstr);
        if (ftell(flog)>MAXLOGSIZE) {
            fclose(flog);
            if (rename(logfilename1,logfilename2)) {
                remove(logfilename2);
                rename(logfilename1,logfilename2);
            }
        } else {
            fclose(flog);
        }
    }
}
void Log(const char *pszFmt,...) {
    va_list argp;

    Lock(&cs_log);
    va_start(argp,pszFmt);
    LogV(pszFmt,argp);
    va_end(argp);
    Unlock(&cs_log);
}
//Log}
int main(int argc,char * argv[]) {
    int i;
#ifdef WIN32
    InitializeCriticalSection(&cs_log);
#else
    pthread_mutex_init(&cs_log,NULL);
#endif
    for (i=0;i<10000;i++) {
        Log("This is a Log %04d from FILE:%s LINE:%d\n",i, __FILE__, __LINE__);