| 10 | 1/1 | 返回列表 |
| 查看: 1950 | 回復: 9 | |||
swustxxl金蟲 (正式寫手)
|
[交流]
【求助】VS2005數(shù)值計算程序調試【急求!。!加金幣】
|
|
數(shù)據(jù)擬合程序,運用現(xiàn)成的動態(tài)鏈接庫LBFGS法(帶約束的),操作平臺是VS2005,將其引入MFC單文檔錯誤迭出 我將算法程序編入了CLoginDlg.CPP中,點擊button1,程序關閉對話框彈出單文檔界面,后臺同時完成運算。 void CLoginDlg::OnBnClickedButton1() { // TODO: 在此添加控件通知處理程序代碼 _tmain(); CDialog::OnOK(); } 以上是Cbutton的代碼,編譯首先報錯的是_tmain()位置,說的是不接受0個參數(shù),這個公共函數(shù)的原型為 int _tmain(int argc, _TCHAR* argv[]) { ......... ......... return 0; } 我干脆直接對原函數(shù)做了改動,如下: int _tmain() { ......... ......... return 0; } 在解決方案中又添加了lbfgsb.lib文件,于是編譯、鏈接都可通過,甚至可以運行,但點擊button1,程序就會報錯,彈出如下對話框: ![]() 算法的最初程序是正確的,可編譯可鏈接可運行,結果正確;我可以新建一個WIN32控制臺應用程序重現(xiàn)算法程序的效果。證明算法程序本身無誤,應該是動態(tài)鏈接庫在加載到MFC程序時出現(xiàn)了錯誤。 希望有高手能指點一下~~不勝感激 以下是我的MFC程序關于引入算法動態(tài)鏈接庫和執(zhí)行運算的片段,最后附上了WIN32算法程序和我自己的MFC程序: // LoginDlg.cpp : 實現(xiàn)文件 // #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中,計算目標函數(shù)值保存在f中 //n是參數(shù)個數(shù),x是參數(shù)向量 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 ] |
金蟲 (正式寫手)


木蟲 (正式寫手)
金蟲 (正式寫手)
|
我在解決方案中添加了一個相關的.lib文件,還是不行。 我把所有的東西還原到了最原始的形態(tài)。那個正確的算法程序采用的是WIN32控制臺應用程序,我新建了一個空的WIN32項目,只把算法的CPP和H文件添進去,把動態(tài)鏈接相關的文件拷貝到相關文件夾中,同樣也添加了.lib文件,一切OK了。 應該是MFC程序和WIN32程序之間的差異導致的報錯,請問現(xiàn)在該如何解決? 我的MFC程序是:MFC單文檔,靜態(tài)庫,不使用UNCOIDE。希望能對你判斷問題所在有所幫助。 [ Last edited by swustxxl on 2010-10-17 at 18:45 ] |
金蟲 (正式寫手)
木蟲 (正式寫手)
|
1、你的問題,可以考慮去CSDN問問, 2、是否排除數(shù)據(jù)類型的問題? 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 ] |
金蟲 (正式寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿a區(qū)211,085601-307分求調劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|---|---|---|---|---|
|
[考研] 285求調劑 +3 | FZAC123 2026-03-30 | 3/150 |
|
|
[考研] 環(huán)境工程 085701,267求調劑 +13 | minht 2026-03-29 | 13/650 |
|
|
[考研] 機械學碩總分317求調劑。。! +6 | Acaciad 2026-03-25 | 6/300 |
|
|
[考研] 340求調劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研] 材料科學與工程調劑 +13 | 深V宿舍吧 2026-03-30 | 14/700 |
|
|
[考研] 08工科,295,接受跨專業(yè)調劑 +4 | lmnlzy 2026-03-31 | 4/200 |
|
|
[考研] 生物學 296 求調劑 +7 | 朵朵- 2026-03-26 | 9/450 |
|
|
[考研] 調劑申請 +6 | 張張張張zy 2026-03-31 | 6/300 |
|
|
[考研] 調劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 一志愿大連理工大學,機械工程學碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[考研] 297 地理學070500 復試求調劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
[考研] 一志愿北京化工大學材料與化工(085600)296求調劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
[考研] 材料科學與工程 317求調劑 +7 | JKSOIID 2026-03-26 | 7/350 |
|
|
[考研] 328求調劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 085600,材料與化工321分求調劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 2026年華南師范大學歡迎化學,化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 315分求調劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 0703化學338求調劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|