24小时热门版块排行榜

返回列表

本帖产生 1 个程序强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

magic7004

金虫 (职业作家)

程序强帖: 2
应助: 0 (幼儿园)
金币: 444
红花: 4
帖子: 3278
在线: 156.2小时
虫号: 482895
注册: 2007-12-26
性别: GG
专业: 化学工程

[交流] 【交流】VC++, C#, VB答疑专贴已有69人参与

帖主寄语

其实我是业余爱好者，不搞研究也不写论文，专业和工作都和编程基本没关系。只是喜欢胡搞瞎搞而已，所以对Windows编程比较熟悉，数值计算之类的就不懂了。

比较熟悉Windows编程，常用VB、VC++、C#。关于语言、语法和编译的问题，欢迎大家一起来交流讨论~.~

[ Last edited by 波不动 on 2010-1-8 at 19:44 ]

回复此楼

» 本帖已获得的红花（最新10朵）

飞虎

bbettry

谭留香

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有21人回复
不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

流氓不可怕，可怕的是流氓有文化，有文化又BH的流氓无敌~~！

1楼 2009-11-18 14:12:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

swustxxl

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1424.8
散金: 211
红花: 2
帖子: 398
在线: 126.8小时
虫号: 954876
注册: 2010-02-11
专业: 药物设计与药物信息

【求助】VS2005数值计算程序调试报错

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

数据拟合程序，运用现成的动态链接库LBFGS法（带约束的），操作平台是VS2005，将其引入MFC单文档错误迭出——不是让楼主处理数值计算哈！！！！

我将算法程序编入了CLoginDlg.CPP中，点击button1，程序关闭对话框弹出单文档界面，后台同时完成运算。

void CLoginDlg::OnBnClickedButton1()
{
      // TODO: 在此添加控件通知处理程序代码
      _tmain();
      CDialog::OnOK();

}
以上是Cbutton的代码，编译首先报错的是_tmain()位置，说的是不接受0个参数，这个公共函数的原型为
int _tmain(int argc, _TCHAR* argv[])
{
   .........
   .........
   return 0;
}
我干脆直接对原函数做了改动，如下：
int _tmain()
{
   .........
   .........
   return 0;
}
在解决方案中又添加了lbfgsb.lib文件，于是编译、链接都可通过，甚至可以运行，但点击button1，程序就会报错，弹出如下对话框：

算法的最初程序是正确的，可编译可链接可运行，结果正确；我可以新建一个WIN32控制台应用程序重现算法程序的效果。证明算法程序本身无误，应该是动态链接库在加载到MFC程序时出现了错误。
希望有高手能指点一下~~不胜感激

以下是我的MFC程序关于引入算法动态链接库和执行运算的片段，最后附上了WIN32算法程序和我自己的MFC程序：

// LoginDlg.cpp : 实现文件
//

#include "stdafx.h"
#include "xxl_try2.h"
#include "LoginDlg.h"

//#include "stdafx.h"
#include "iostream"
#include "memory.h"
#include "malloc.h"
#include "windows.h"
#include "math.h"

using namespace std;

// CLoginDlg 对话框

IMPLEMENT_DYNAMIC(CLoginDlg, CDialog)

CLoginDlg::CLoginDlg(CWnd* pParent /*=NULL*/)
: CDialog(CLoginDlg::IDD, pParent)
{

}

CLoginDlg::~CLoginDlg()
{
}

void CLoginDlg

DataExchange(CDataExchange* pDX)
{
CDialog

DataExchange(pDX);
}

BEGIN_MESSAGE_MAP(CLoginDlg, CDialog)
      ON_BN_CLICKED(IDC_BUTTON1, &CLoginDlg::OnBnClickedButton1)
END_MESSAGE_MAP()

// CLoginDlg 消息处理程序
// LBFGS.cpp : 定义控制台应用程序的入口点。
//计算梯度保存在g中，计算目标函数值保存在f中
//n是参数个数，x是参数向量
void fcn2(int n,double *x, double*t, double *y,double &f, double *g ){
      f=0;
      g[0]=0;
      g[1]=0;
      g[2]=0;
      g[3]=0;
      g[4]=0;
      g[5]=0;
      for(int k=0;k       {
            double tmp=y[k]-(x[0]*t[k]*exp(x[1]*t[k])+x[2]*t[k]*exp(x[3]*t[k])+x[4]*t[k]*exp(x[5]*t[k]));
            f+=tmp*tmp;
            g[0]+=-2*tmp*t[k]*exp(x[1]*t[k]);
            g[1]+=-2*tmp*t[k]*x[0]*exp(x[1]*t[k])*t[k];
            g[2]+=-2*tmp*t[k]*exp(x[3]*t[k]);
            g[3]+=-2*tmp*t[k]*x[2]*exp(x[3]*t[k])*t[k];
            g[4]+=-2*tmp*t[k]*exp(x[5]*t[k]);
            g[5]+=-2*tmp*t[k]*x[4]*exp(x[5]*t[k])*t[k];
      }

      return;
}

extern "C" void __declspec(dllimport)  __stdcall lbfgsb(int n,int m,double *x, double *l, double *u, int *nbd, double f, double *g, double factr, double pgtol, double *wa, int *iwa,char *task, int iprint,char *csave,bool *lsave,int *isave, double *dsave);

int main()//int argc, _TCHAR* argv[])
{
      const char dllPath[] = "lbfgsb.dll";

      //加载动态库文件
      HINSTANCE hLibrary = LoadLibraryA(dllPath);
      if(hLibrary == NULL){
            int error =GetLastError ();
            cout<<"Error of LoadLibrary = "<
            cout<<"Can't find the dll file."<             FreeLibrary(hLibrary);
            return -1;
      }

      /* （如果获得这个程序的句柄）获得正在运行程序的路径和文件名*/
      bool exePathOut=true;

      if(exePathOut){
            char exePath[1000] = "\0";
            GetModuleFileNameA(hLibrary,exePath,sizeof(exePath));
            cout<       }

      bool *lsave;
      char *task,*csave,*str;
      int n, m,iprint,*nbd,*iwa, *isave;
      double f, factr, pgtol,*x,*l,*u,*g,*wa,*dsave;
      double t[7]={1,1.5,2,2.5,3,3.5,5};
      double y[7]={10.78,80.35,58.81,43.04,23.05,12.35,6.61};

      iprint = 1;
      factr=1.0E+7;
      pgtol=1.0E-5;
      f=0.0;

      n=6;
      m=10;

      nbd=new int[n];
      iwa=new int[3*n];
      x=new double[n];
      l=new double[n];//设置变量下限
      u=new double[n];//设置变量上限
      g=new double[n];
      wa=new double[2*m*n+4*n+12*m*m+12*m];
      dsave=new double[29];
      isave=new int[44];
      str=new char[60];
      csave=new char[60];
      task=new char[60];
      lsave=new bool[16];

      for(int i=0;i       {
            nbd=0;
            l=1.0;
            u=10;
      }

      for(int i=1;i       {
            nbd=0;
            l=1.0;
            u=10;
      }

      for(int i=0;i       {
            x=0.0E0;
      }

      sprintf(task,"%s","START";

      while(1){
            lbfgsb(n,m,x,l,u,nbd,f,g,factr,pgtol,wa,iwa,task,iprint,csave,lsave,isave,dsave);
            strncpy(str,task,2);
            if (strncmp(str,"FG",2)==0)
            {
                     fcn2(7, x,t,y,f, g );
                     continue;
            }
            strncpy(str,task,5);
            if(strncmp(str,"NEW_X",5)!=0)
            {
                     for(int i=0;i                      {
                              cout<<"x["<                      }
                     break;
            }
      }

      delete[] nbd;
      delete[] iwa;
      delete[] x;
      delete[] l;
      delete[] u;
      delete[] g;
      delete[] wa;

      delete[] lsave;
      delete[] task;
      delete[] csave;
      delete[] str;
      delete[] isave;
      delete[] dsave;

      system("pause";
      return 0;
}

void CLoginDlg::OnBnClickedButton1()
{
      // TODO: 在此添加控件通知处理程序代码
      main();
      CDialog::OnOK();

}

这是程序
http://good.gd/767141.htm
这是算法的正确程序
http://good.gd/767144.htm

本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2492928

[ Last edited by swustxxl on 2010-10-20 at 20:39 ]