Построение гистограмм, графиков, тепловых карт

[Iryna Kapran]

С++ делает подсчёт и выдаёт данные, но возникает трудность в создании гистограммы, графика и тепловой карты на основе полученных результатов. Вот код для получения данных:

#include <iostream>
#include <iomanip>
#include <cmath>

int main() {
int S[20*20];
for (int i = 0; i < 20; i++) {
for (int j = 0; j < 20; j++) {
if (std::rand() % 100 >= 50) { S[20 * i + j] = 1; }
else { S[20 * i + j] = -1; }
}
}
for (int i = 0; i < 20; i++) { // Here we show an array.
for (int j = 0; j < 20; j++) {
std::cout << std::setw(3) << S[20 * i + j] << " ";
}
std::cout << std::endl;
}
int T = 273; // Temperature.
for (unsigned int k = 0; k < 10; k++) {
int J = 0, r = 0;
r = std::rand()%400;
if ((r + 1 >= 0) && (r + 1 <= 400)) {
if (S[r] == S[r + 1]) { J--; }
else if (S[r] != S[r + 1]) { J++; }
}
if ((r - 1 >= 0) && (r - 1 <= 400)) {
if (S[r] == S[r - 1]) { J--; }
else if (S[r] != S[r - 1]) { J++; }
}
if ((r + 20 >= 0) && (r + 20 <= 400)) {
if (S[r] == S[r + 20]) { J--; }
else if (S[r] != S[r + 20]) { J++; }
}
if ((r - 20 >= 0) && (r - 20 <= 400)) {
if (S[r] == S[r - 20]) { J--; }
else if (S[r] != S[r - 20]) { J++; }
}
int E1 = J; J = 0;
if ((r + 1 >= 0) && (r + 1 <= 400)) {
if (S[r] == -S[r + 1]) { J--; }
else if (S[r] != -S[r + 1]) { J++; }
}
if ((r - 1 >= 0) && (r - 1 <= 400)) {
if (S[r] == -S[r - 1]) { J--; }
else if (S[r] != -S[r - 1]) { J++; }
}
if ((r + 20 >= 0) && (r + 20 <= 400)) {
if (S[r] == -S[r + 20]) { J--; }
else if (S[r] != -S[r + 20]) { J++; }
}
if ((r - 20 >= 0) && (r - 20 <= 400)) {
if (S[r] == -S[r - 20]) { J--; }
else if (S[r] != -S[r - 20]) { J++; }
}
if (J - E1 < 0) { S[r] = -S[r]; }
else if (J - E1 > 0) {
if (std::rand() % 100 >= (3 ^ ((E1 - J) / T))) {
S[r] = -S[r];
}
}
std::cout << std::endl;
for (int i = 0; i < 20; i++) { // Here we show an array.
for (int j = 0; j < 20; j++) {
std::cout << std::setw(3) << S[20 * i + j] << " ";
}
std::cout << std::endl;
}
}
}

[Пётр Седов]

Сообщение от Iryna Kapran

С++ делает подсчёт и выдаёт данные, но возникает трудность в создании гистограммы, графика и тепловой карты на основе полученных результатов.

Как это всё должно быть нарисовано по массиву S, в котором только 1 и -1? Непонятно же. Тут программисты ошиваются, а не физики.

Сообщение от Iryna Kapran

Вот код для получения данных:

Который надо окружать tag-ами [ CODE ] [ / CODE ] (без пробелов).

Сообщение от Iryna Kapran

		if ((r + 1 >= 0) && (r + 1 <= 400)) {

Второе граничное условие должно быть «r + 1 < 400», потому что в C++ индекс должен быть строго меньше длины массива. Вообще, здесь лучше сделать функцию:

bool index_is_valid(int index) {
  return (0 <= index) && (index < 400);
}

и потом её использовать во всех таких if-ах.

Сообщение от Iryna Kapran

			if (S[r] == S[r + 1]) { J--; } 
			else if (S[r] != S[r + 1]) { J++; }

Так не надо писать, надо так:

if (S[r] == S[r + 1]) { J--; } 
else { J++; }

В else мы как раз и попадём, если элементы не равны.

Сообщение от Iryna Kapran

			if (std::rand() % 100 >= (3 ^ ((E1 - J) / T))) {

Это вы так пытаетесь число 3 возвести в степень «(E1 - J) / T»? В C++, ^ -- это битовый xor, а возведение в степень делается функцией pow из <math.h> (она же std::pow из <cmath>). Кстати, E1, J и T -- это переменные типа int, поэтому деление будет целочисленное, с округлением.