finityelement.cpp

来自「偏微分方程数值解 有限元法 面向对象 变分问题 剖分问题 边值处理 误差分析 椭」· C++ 代码 · 共 75 行

#include "Onedimension.h"
#include "Twodimension.h"
//#include "teaminfo.h"
#include<time.h>
#include <stdio.h>
#include <conio.h>
int main()
{
	char ch;
	clock_t start,end;
	double t=0;
	while (true)
	{
		printf("*********************************************************\n");
		printf("欢迎进入第七实验小组有限元分析演示程序……\n");
		printf("              1:一维问题 \n");
        printf("              2:二维问题 \n");
		printf("              0:退出 \n");
		printf("*********************************************************\n");
		printf("请选择(1, 2 或者 0)\n");
		do
		{
			ch=(char) _getch();
		}while (ch!='1' && ch!='2' && ch!='0');
		system("cls");
		if (ch=='0')
			return 0;
		else if (ch=='1')
		{	
			cout<<"------------------------------一维有限元问题开始";
	        cout<<"--------------------------------"<<endl;
			Onedimension a;
			start=clock();
	        a.Dispart();
	        a.Radix();
	        a.Calculate();
	        a.Exact();
			cout<<endl<<"步长h="<<h<<endl;
			end=clock();
			t=(double)(end-start);
   			cout<<endl<<"计算所用时间t="<<t<<"毫秒"<<endl<<endl;
	        a.Error();
			cout<<endl;
			cout<<"------------------------------一维有限元问题结束";
	        cout<<"--------------------------------"<<endl;
		}
		else
		{
			cout<<"------------------------------二维有限元问题开始";
	        cout<<"--------------------------------"<<endl;
 			Twodimension b;
            start=clock();
	        b.Dispart();
	        b.Radix();
	        b.Calculate();
	        b.Exact();
			cout<<endl<<"步长H="<<H<<endl;
			end=clock();
			t=(double)(end-start);
   			cout<<endl<<"计算所用时间T="<<t<<"毫秒"<<endl<<endl;
	        b.Error();
			cout<<endl;
			cout<<"------------------------------二维有限元问题结束";
	        cout<<"--------------------------------"<<endl;
		}
        //get_teaminfo();
		printf("\nPress any key to continue\n");
		_getch();
		system("cls");
	}
	return 0;
}