随着信息时代的到来，需要处理的信息量越来越庞大，需要解决的问题越来越复杂，使得计算量剧增。通过提高单个处理器的计算速度和采用传统的“顺序（串行）”计算技术已难以胜任。因此，需要有功能更强大的计算机系统和计算机技术来支撑。并行计算机及并行计算技术应运而生。
但由于缺乏实验环境和机器设备，我们普通人很难研究并行算法，即使有了想法也同样面临着无法验证的尴尬。不过，好在像面向对象语言c++，java等都提供了多线程，使我们可以模拟多台处理机。下面，我就一个简单的例子，向大家介绍一下在vc下如何利用多线程模拟多处理机并行求取最大值问题。
题目如下：令n＝2的m次方，A是一个2n维的数组，待求最大值的数存放在A(n),A(n＋1)，……A(2n-1)中，所求得的最大值置于A(1),于是算法描述如下：
输入：n＝2的m次方个数存在数组A(n；2n－1)中；
输出：最大数置于A(1)中。
Begin
For k=m-1 to 0 do
For j=2 to 2n-1 par do
A(j)=max(A(2j)，A(2j+1))
End For
End For
End
|
显然，算法的时间t(n)=O(lgn),总比较次数为O(n)，而最大的处理器数p(n)=n/2。（也即最大的线程数。）
首先，我们建立一个基于对话框的应用程序，然后放置如图所示控件：
然后定义全局变量如下，array用于存放产生的随机数。
int g_nCount=0,j=0;
int *array;
| |
在对话框的初始化函数中，我们随机产生n个数并存于数组array中，程序代码如下：
void CDemoDlg::OnBtnInitial()
{
// 初始化数组
srand((unsigned)time(NULL));
int i,temp;
CString str;
UpdateData(true);
g_nCount=pow(2,m_intCount);
array=new int[2*g_nCount];
//根据用户的输入，产生2的m次方的随机数并存于array[n]…a[2n-1]中
for(i=g_nCount;i<2*g_nCount;i++)
{
temp=rand()/100;
array[i]=temp;
}
//显示产生的n个随机数
for(i=g_nCount;i<2*g_nCount;i++)
{
str.Format("Array[%d]= %d",i,array[i]);
m_strArray+="\r\n";
m_strArray+= str;
m_strArray+="\r\n";
}
m_strArray+="\r\n";
UpdateData(false);
} void CDemoDlg::OnBtnCompute()
{
// 启动线程进行计算
int k;
for(k=m_intCount-1;k>=0;k--)
{
for(j=pow(2,k);j {
AfxBeginThread(ComputeThread,GetSafeHwnd(),
THREAD_PRIORITY_NORMAL);
}
}
//线程体：较array[2j]与array[2j+1]的大小，将大值置于array[j]中
UINT ComputeThread(LPVOID pParam)
{
if(array[2*j]>=array[2*j+1])
{
array[j]=array[2*j];
}
else
{
array[j]=array[2*j+1];
}
return 0;
}
}
//显示计算结果
void CDemoDlg::OnBtnShow()
{
// TODO: Add your control notification handler code here
int k;
CString str;
for(k=1;k<2*g_nCount;k++)
{
str.Format("Array[%d]＝ %d",k,array[k]);
m_strArray+="\r\n";
m_strArray+= str;
m_strArray+="\r\n";
} UpdateData(false);
}
|
运行结果如图：
从运行的结果图上我们可以清楚的看到线程的计算排序过程。
此程序只是一个利用多线程进行并行计算的简单例子，希望对各位进行并行算法的研究有所帮助。
查看本文来源
京公网安备 11010802021500号