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);