| 5 | 1/1 | 返回列表 |
| 查看: 1951 | 回復: 9 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
swustxxl金蟲 (正式寫手)
|
[交流]
【求助】VS2005數值計算程序調試【急求。。。〖咏饚拧
|
||
|
數據擬合程序,運用現成的動態(tài)鏈接庫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控制臺應用程序重現算法程序的效果。證明算法程序本身無誤,應該是動態(tài)鏈接庫在加載到MFC程序時出現了錯誤。 希望有高手能指點一下~~不勝感激 以下是我的MFC程序關于引入算法動態(tài)鏈接庫和執(zhí)行運算的片段,最后附上了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: oDataExchange(CDataExchange* pDX){ CDialog: oDataExchange(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"; //加載動態(tài)庫文件 HINSTANCE hLibrary = LoadLibraryA(dllPath); if(hLibrary == NULL){ int error =GetLastError (); cout<<"Error of LoadLibrary = "< cout<<"Can't find the dll file."< 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 [ Last edited by swustxxl on 2010-10-20 at 10:47 ] |
|
1、你的問題,可以考慮去CSDN問問, 2、是否排除數據類型的問題? 3、也有可能是fortran編譯器版本的問題,如:http://software.intel.com/en-us/forums/showthread.php?t=73197 一般情況,google比baidu好用, 以 Visual Fortran run-time error 和 fortl: severe (38): error during write 搜索,看英文的網頁, 引起這個錯誤的因素較多, 可以一個個排除。 [ Last edited by gd88 on 2010-10-20 at 22:55 ] |


木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 359求調劑 +6 | 鄧鄧鄧書書 2026-03-25 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料科學與工程 317求調劑 +8 | JKSOIID 2026-03-26 | 8/400 |
|
|
[考研] 375求調劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 346求調劑 一志愿070303有機化學 +11 | 蘿卜燉青菜 2026-03-28 | 12/600 |
|
|
[考研] 070300一志愿211,312分求調劑院校 +12 | 小黃鴨寶 2026-03-30 | 12/600 |
|
|
[考研] 254材料與化工求調劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 080500-315分復試調劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
[考研] 271求調劑 +5 | 勒布朗@ 2026-03-31 | 7/350 |
|
|
[考研] 一志愿085600中科院寧波所276分求調劑 +12 | 材料學257求調劑 2026-03-28 | 13/650 |
|
|
[考研]
材料調劑
5+5
|
想要一壺桃花水 2026-03-25 | 11/550 |
|
|
[考研] 085600材料與化工調劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
[考研] 293分求調劑,外語為俄語 +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 359求調劑 +5 | 王了個楠 2026-03-25 | 5/250 |
|
|
[考研] 0703化學321分求調劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研] 26考研-291分-廈門大學(085601)-柔性電子學院材料工程專業(yè)求調劑 +5 | min3 2026-03-24 | 6/300 |
|
|
[考研] 材料專碩 085600求調劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 數一英一271專碩(085401)求調劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 316求調劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|