|
13.05.2014, 21:35
#38640286
Ссылка:
Ссылка на сообщение:
Ссылка с названием темы:
|
|
|
|
Где-то в степи, (прилагаю код еще раз)
1. Да код умножает верно,но число процессов желательно уменьшить!
2. Спасибо вам за ссылку, щас буду изучать англоязычные сайты)
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
83.
84.
85.
86.
87.
88.
89.
90.
91.
92.
93.
94.
95.
96.
97.
98.
99.
100.
101.
102.
103.
104.
105.
106.
107.
108.
109.
110.
111.
112.
113.
114.
115.
116.
117.
118.
119.
120.
121.
122.
123.
124.
125.
126.
127.
128.
129.
130.
131.
132.
133.
134.
135.
136.
137.
138.
139.
140.
141.
142.
143.
144.
145.
146.
147.
148.
149.
150.
151.
152.
153.
154.
155.
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using MPI;
using System.IO;
namespace ConsoleApplication1
{
class Program
{
static void Main(string[] args)
{
using (new MPI.Environment(ref args))
{
Intracommunicator comm = Communicator.world;
int rank = comm.Rank;
int size = comm.Size;
//double[,] A = { { 1, 1, 1 }, { 1, 1, 1 }, { 1, 1, 1 } };
//double[,] B = { { 1, 1, 1 }, { 1, 1, 1 }, { 1, 1, 1 } };
double[,] A = SchitIzFaila("a.txt");
double[,] B = SchitIzFaila("b.txt");
int razmer = A.GetLength(0);
double[,] C = new double[razmer, razmer];
double[] tempMatrix = new double [size];
if (rank == 0)
{
Console.WriteLine("Matrica A");
vyvodvconsol(A);
Console.WriteLine("Matrica B");
vyvodvconsol(B);
}
if (rank == 0)
{
int IndexRank = 1;
for (int i = 0; i < razmer; i++)
{
for (int j = 0; j < razmer; j++)
{
{
int[] koordinata = new int[2] { i, j };
comm.Send(koordinata, IndexRank, 0);
IndexRank++;
// comm.Wait();
}
}
}
comm.Gather(-1, 0, ref tempMatrix);
comm.Barrier();
C = PoluchMatricuIzPosledovat(tempMatrix, razmer);
Console.WriteLine("Matr C");
vyvodvconsol(C);
}
else
{
int[] koordinata = new int[2];
comm.Receive(0, 0, ref koordinata);
double resultat = Element(PoluchitStroku(koordinata[0], A), PoluchitStolbec(koordinata[1], B));
comm.Gather(resultat, 0, ref tempMatrix);
comm.Barrier();
}
}
}
public static double[,] SchitIzFaila(string adress)
{
StreamReader sr = new StreamReader(adress);
int columns = 0;
int lines = 0;
while (!sr.EndOfStream)
{
string temp = sr.ReadLine();
columns = temp.Split(' ').Count();
lines++;
}
double[,] result = new double[lines, columns];
sr = new StreamReader(adress);
for (int i = 0; i < result.GetLength(0); i++)
{
string[] temp = sr.ReadLine().Split(' ');
for (int j = 0; j < result.GetLength(1); j++)
{
result[i, j] = Convert.ToDouble(temp[j]);
}
}
return result;
}
public static void vyvodvconsol(double[,] matrix)
{
Console.WriteLine();
for (int i = 0; i < matrix.GetLength(0); i++)
{
for (int j = 0; j < matrix.GetLength(1); j++)
{
Console.Write(matrix[i, j]);
Console.Write(" ");
}
Console.Write('\n');
}
Console.WriteLine();
}
public static double[] PoluchitStroku(int idLine, double[,] M)
{
double[] result = new double[M.GetLength(1)];
for (int i = 0; i < result.Length; i++)
{
result[i] = M[idLine, i];
}
return result;
}
public static double[] PoluchitStolbec(int idColumn, double[,] M)
{
double[] result = new double[M.GetLength(0)];
for (int i = 0; i < result.Length; i++)
{
result[i] = M[i,idColumn];
}
return result;
}
public static double Element(double[] line, double[] column)
{
double result = 0;
for (int i = 0; i < line.Length; i++)
{
result += line[i] * column[i];
}
return result;
}
public static double[,] PoluchMatricuIzPosledovat(double[] matrix, int size)
{
double[,] resmatrix = new double[size, size];
int matrixIndex = 1;
for (int i = 0; i < size; i++)
{
for (int j = 0; j < size; j++)
{
resmatrix[i, j] = matrix[matrixIndex];
matrixIndex++;
}
}
return resmatrix;
}
}
}
|
|
|
