Форум программистов, компьютерный форум CyberForum.ru

Немного из архитектуры ЭВМ - C++

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 12, средняя оценка - 4.67
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
29.04.2011, 15:52     Немного из архитектуры ЭВМ #1
Пусть заданы две квадратных матрицы A и B размером NxN.
Они созданы с помощью двух подходов:

1 подход:
C++
1
2
3
4
int **A;
A = new int*[N];
for(int i=0;i<N;i++)
A[i] = new int[N];
Доступ к элементу: A[i][j]

2 подход:
C++
1
2
int *A;
A = new int[N*N];
Доступ к элементу: A[i*N+j]

Необходимо сложить эти матрицы и сравнить время вычисления. Нужны выводы по этим подходам: количество обращений к памяти,
вычислений, теоретические оценки времени вычислений (можно считать, что
матрицы не загружаются в кэш, время доступа к одному элементу в памяти 10
тактов, арифметическая операция 2 такта), реальное время вычисления.

Кто свободен, может поможите?.
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
29.04.2011, 18:54  [ТС]     Немного из архитектуры ЭВМ #21
вот это libname.lib ?
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Deviaphan
Делаю внезапно и красиво
Эксперт C++
 Аватар для Deviaphan
1283 / 1217 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
29.04.2011, 18:56     Немного из архитектуры ЭВМ #22
Ну, как я и написал "С некоторой долей уверенности"... Поэтому, доверяй, но проверяй!
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
29.04.2011, 19:48  [ТС]     Немного из архитектуры ЭВМ #23
спасибо Вам всем. Но что-то до меня мало что дошло. Код я поняла, но не пойму почему не работает.
Deviaphan
Делаю внезапно и красиво
Эксперт C++
 Аватар для Deviaphan
1283 / 1217 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
30.04.2011, 10:26     Немного из архитектуры ЭВМ #24
Цитата Сообщение от KuKu Посмотреть сообщение
12979 и 7427
Да как так вообще??? У тебя Спектрум что ли? У меня даже под отладчиком получилось 150 миллисекунд. И первый вариант чуть быстрее (как я и ожидал).
В любом случае, на таком коде проводить сравнение не корректно, т.к. компилятор может чего-нить на оптимизировать и результат будет искажённым. (это исходя из того, что релизная версия отработала за 100 миллисекунд. Малая разница в производительности настораживает)
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
30.04.2011, 17:51  [ТС]     Немного из архитектуры ЭВМ #25
Слушайте, вот че я думаю. Количество обращений к памяти (запрос на чтение из памяти, запись в память или поиск в памяти заданной информации) = 3. Это обращение к массиву С, В и А.
Кол-во вычислений = 1. Только сложение.
Т.к. за один обход цикла обращений к памяти 3, то число тактов обращений 30. Число тактов вычислений - 2. Всего 32 такта за один обход цикла.
N*N*32 - число тактов за весь цикл.

Это правильно?

Добавлено через 1 час 16 минут
хотя вот в принципе вроде прога.

C
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
#include "stdafx.h"
#include "stdio.h"
#include "conio.h"
#include "windows.h"
#define N 1024
 
 
void main(){
    LARGE_INTEGER b_start,b_stop,b_time,freq;
    QueryPerformanceFrequency(&freq);
    printf("t/sec=%lf*10^9\n\n",(double)(freq.QuadPart)/1000000000);
 
    //Первый подход
    int **A1;
    A1 = new int*[N];
    for(int i=0;i<N;i++){
        A1[i] = new int[N];
        for(int j=0;j<N;j++)
            A1[i][j]=j*N+i;
    }
    int **B1;
    B1 = new int*[N];
    for(int i=0;i<N;i++){
        B1[i] = new int[N];
        for(int j=0;j<N;j++)
            B1[i][j]=j*N+i;
    }
    int **C1;
    C1 = new int*[N];
    for(int i=0;i<N;i++)
        C1[i] = new int[N];
 
    QueryPerformanceCounter(&b_start);
    for(int i=0;i<N;i++){
        for(int j=0;j<N;j++){
            C1[i][j]=A1[i][j]+B1[i][j];
        }
    }
    QueryPerformanceCounter(&b_stop);
    b_time.QuadPart = b_stop.QuadPart - b_start.QuadPart;
 
    printf("t=%lf\n",(double)(b_time.QuadPart));
    printf("time= \t%lfms\n", (double)(b_time.QuadPart)*1000.0/(double)(freq.QuadPart));
 
 
 
    //Второй подход
    int *A;
    A = new int[N*N];
    for(int i=0;i<N;i++)
        for(int j=0;j<N;j++)
            A[i*N+j]=i*N+j;
    int *B;
    B = new int[N*N];
    for(int i=0;i<N;i++)
        for(int j=0;j<N;j++)
            B[i*N+j]=j*N+i;
    int *C;
    C = new int[N*N];
 
    
    QueryPerformanceCounter(&b_start);
    for(int i=0;i<N*N;i++){
        C[i]=A[i]+B[i];
    }
    QueryPerformanceCounter(&b_stop);
    b_time.QuadPart = b_stop.QuadPart - b_start.QuadPart;
 
    printf("t=%lf\n",(double)(b_time.QuadPart));
    printf("time= \t%lfms\n\n", (double)(b_time.QuadPart)*1000.0/(double)(freq.QuadPart));
 
    _getch();
}
 Комментарий модератора 
Используйте теги форматирования кода!
prZzZ
Программист
 Аватар для prZzZ
108 / 88 / 5
Регистрация: 22.03.2011
Сообщений: 232
30.04.2011, 18:19     Немного из архитектуры ЭВМ #26
Цитата Сообщение от oksanaBM Посмотреть сообщение
QueryPerformanceCounter
Кто так время измеряет? Или у вас Pentium II? Про счетчик меток реального времени никогда не слышали?
C++
1
2
3
4
unsigned __int64 GetCycleCount()
{
    _asm rdtsc
}
Правда его надо градуировать под каждый камень. Очень уж тут всё аппаратно зависимо
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
30.04.2011, 18:25  [ТС]     Немного из архитектуры ЭВМ #27
по учебной программе нужно использовать QueryPerformanceCounter
prZzZ
Программист
 Аватар для prZzZ
108 / 88 / 5
Регистрация: 22.03.2011
Сообщений: 232
30.04.2011, 18:39     Немного из архитектуры ЭВМ #28
Точности большой тогда у вас не будет. Её и с вышеописанным методом не будет т.к. Windows НЕ является системой реального времени, но такты процессора (над этим, кстати, здесь смеялись) посчитать получилось бы. Извините...
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
30.04.2011, 18:41  [ТС]     Немного из архитектуры ЭВМ #29
получается программа неверная?
prZzZ
Программист
 Аватар для prZzZ
108 / 88 / 5
Регистрация: 22.03.2011
Сообщений: 232
30.04.2011, 18:49     Немного из архитектуры ЭВМ #30
Да, нет... Теоретически верно...
А практически:
1)погрешность операционной системы (тут, увы, мы ничего сделать не можем)
2)погрешность компилятора (на уровень ассемблера опускаться тоже не входит в ваши планы)
3)погрешность инструментальных средств (описанных выше)

Всё поплюсуем... Истины добиться очень не просто...
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
30.04.2011, 19:01     Немного из архитектуры ЭВМ
Еще ссылки по теме:

C++ теория вероятности на эвм
C++ Составление архитектуры программы
ОС и архитектура ЭВМ C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
silent_1991
30.04.2011, 19:01     Немного из архитектуры ЭВМ
  #31
 Комментарий модератора 
Поскольку от ТСа поступила информация о том, что на все интересующие его вопросы ответы получены - тему закрываю.
Yandex
Объявления
30.04.2011, 19:01     Немного из архитектуры ЭВМ
Закрытая тема Создать тему
Опции темы

Текущее время: 11:43. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru