分类: C/C++
2008-08-01 16:58:51
The Migration of Old C Code to Visual C IDE Abstract: According to the character of C programs, this paper presents some techniques to migrate them to Visual C IDE, as a implemention, it reprograms a set of numerical arithmetic programs for further engineering use.关键词:移植;数值计算;封装;模板类
Key words: Migration; Numerical Arithmetic, Encapsulation, Template Class一、引言
int agaus(a,b,n) int n; double a[],b[]; {……;}参数声明应改为数组形式:
int agaus(double a[],double b[],int n)或者改为指针形式:
int agaus(double* a,double* b,int n);调用方法:
agaus(&a[0][0],&b[0],n); /* a二维双精度型数组、b一维双精度型数组,n整型变量 */C/C 中用下标法和指针法都可以访问一个数组,设有数组a,则a[i]和*(a i)无条件等价。如果指针变量p指向数组中的一个元素,则p 1指向同一数组的下一个元素。若p的初值为&a[0],则p i和a I都是a[i]的地址;*(p i)和*(a i)就是p i或a i所指向的数组元素,即a[i];指向数组的指针变量也可以带下标,如p[i]与*(p i)等价。所以,在实际使用该函数,如果遇到数组作形参,可以将数组第一个元素地址作为实参传值调用函数。
double* v; v=malloc(n*m*sizeof(double));/* Line1 */ v=(double*)malloc(n*m*sizeof(double)); /* Line2 */3、某些算法函数可能要调用一些用户自定义函数,如最佳一致逼近的里米兹方法:
void hremz(a,b,p,n,eps) int n; double a,b,eps,p[]; { extern double hremzf(); … }原方法使程序集与应用程序的耦合程度增加,缺乏灵活性,可以改为:
void hremz(double a,double b,double p[],int n,double eps,double (*hremzf)(double x)) {…}用函数指针作参数,调用时直接将函数名作实参即可: hremz(a,b,p,4,eps,hremzf); /* 假设各参数在主程序文件已定义 */
printf("%c",xy[i][j]);我们可以用形似
sprintf( buffer,"%c",xy[i][j]), strcat( str, buffer );的合并语句(其中str是一个足够大的字符串数组参数)代替
printf("%c",xy[i][j]);例如:
char* buffer; buffer =(char*)malloc(n*sizeof(char)); /*n作为参数传递,例如100 */ sprintf( buffer,"%c",xy[i][j]), strcat( str, buffer ); /*把终端输出字符添加到str 串尾*/ ...... free(buffer)如果用到了它们,调用方法以随机样本分析为例:
char str[1024]; str[0]=''\0'';/*初始化为空串*/ irhis(x,100,x0,h,10,1,&dt[0],&g[0],&q[0],str);现在str数组保存了终端输出文本,可以随意使用它,比如在控制台程序里输出:
puts(str);在使用MFC类库时,str可以直接赋值给一个CString对象的实例。 经过以上的工作,我们得到基于ANSI C标准的程序版本,可以在C和C 开发环境下使用。
template举例说明我们关于异常机制和数组越界的检测方法的思路。TMatrix类的operator[]返回一个嵌套类Row的引用,它用来描述某一给定二维数组的一个特定行。Row类的operator[]则返回该行一个特定位置的T类型值。最终实现还是通过Matrixclass TArray { protected: T* pdata; unsigned int length; public: TArray(); TArray(unsigned int); TArray(TArray const&); virtual ~TArray(); void operator = (TArray&); TArray & operator (TArray&); TArray & operator - (TArray&); T const& operator [] (unsigned int)const; T& operator [](unsigned int); T const* GetData() const; unsigned int GetLenght(); void SetLength(unsigned int,bool=true); }; template class TMatrix { protected: unsigned int numberOfRows; unsigned int numberOfColumns; TArray array; public: class Row { TMatrix & matrix; unsigned int const row; public: Row (TMatrix & _matrix,unsigned int _row):matrix(_matrix),row(_row){} T& operator [](unsigned int column)const {return matrix.Select(row,column);} }; TMatrix(); TMatrix(unsigned int, unsigned int); TMatrix(TMatrix & mat); virtual ~TMatrix(); T& Select(unsigned int, unsigned int); Row & operator[](unsigned int); TMatrix & operator (TMatrix & mat); TMatrix & operator - (TMatrix & mat); TMatrix & operator * (TMatrix & mat); bool operator == (TMatrix & mat); TArray & GetData(); unsigned int GetNumberOfRows(); unsigned int GetNumberOfColumns(); bool LoadFromArray(T [],unsigned int,unsigned int); bool LoadFromString(char*,char,char); bool ResetMatrix(unsigned int, unsigned int); bool ReverseMatrix(); void ZeroMatrix(); void RandomMatrix(int max); };
template应用程序实例:T& TMatrix ::Select(unsigned int i, unsigned int j) { char ch[50]; if(i>=numberOfRows) sprintf(ch," Error -- Invalid row: %d",i), throw (ch); if(j>=numberOfColumns) sprintf(ch," Error -- Invalid colum: %d",j), throw (ch); return array[i*numberOfColumns j]; }
unsigned int i,j; unsigned int m=2,n=3; TMatrix< > mat(m,n);//双精度型矩阵 try { for(i=0; i由于我们对operator[]进行了重载,所以数组模板类(矩阵模板类)完全兼容C/C 一维数组(二维数组)的存取操作,因此旧程序中数组变量直接可以用类实例变量替代。 图一是我们用Visual C 6开发的演示程序界面,左边是所有算法的目录树,右边是文本计算结果输出,下部悬浮窗口是算法程序源代码,可以拷贝粘贴,稍作修改即可重用。终端输出如下(注:类实例mat没有初始化): -6.27744e 066 -6.27744e 066 -6.27744e 066 Error -- Invalid colum: 3只输出一行,根据出错提示,把Line3改为:"for(j=0; j-6.27744e 066 -6.27744e 066 -6.27744e 066 -6.27744e 066 -6.27744e 066 -6.27744e 066