Ошибка сегментирования (сделан дамп памяти)

@Dionisius · Регистрация: 19.10.2013

Студворк — интернет-сервис помощи студентам

Доброго времени суток. На текущий момент имеется программа, позволяющая выводить различные вариации кристаллической решетки 3*4 элементов (таких комбинаций должно быть 27 штук).

Кликните здесь для просмотра всего текста

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
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
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <random>
#include <tbb/tbb.h>
#include <tuple>
//#include "mpi.h"
 
using namespace std;
 
//построение кристаллической решетки
double BuildLattice(double a0)
{
int x=-1;
int y=-1;
int z=-1;
double Rij=0;
 
std::vector<std::vector<double>>arr(108, vector<double>(3));
//создание различных вариаций КР
double rasst[27][12];
 
for (int i=0; i<4; i++){
for (int j=0; j<3; j++){
if (i==0) arr[i][j]=0;
arr[i+1][i]=0;
}}
arr[1][1]=0.5; arr[1][2]=0.5; arr[2][0]=0.5; arr[2][2]=0.5; arr[3][0]=0.5; arr[3][1]=0.5; cout<<endl;
std::vector<std::vector<double>>vct(4, vector<double>(3));
 
for (int i1=1; i1<=27; i1++){
int i4=0;
for (int i2=0; i2<4; i2++){
for (int i3=0; i3<3; i3++){
std::cout<<arr[i2][i3]<<" ";
if(i2==i3) rasst[i1][i4]=0;
if ((i2-i3 ==1) || (i3-i2 == 1)) rasst[i2][i4]=1;
if ((i2-i3 ==2) || (i3-i2 == 2)) rasst[i2][i4]=2;
if ((i2-i3 ==3) || (i3-i2 == 3)) rasst[i2][i4]=3;
i4++;
}
std::cout<<endl;
}
cout<<endl;
if ((i1==9) || (i1==18)) {
arr[0][0]=arr[0][0]+1; arr[1][0]=arr[1][0]+1; arr[2][0]=arr[2][0]+1; arr[3][0]=arr[3][0]+1;
}
if ((i1%9==3) || (i1%9 == 6)) {
arr[0][1]=arr[0][1]+1; arr[1][1]=arr[1][1]+1; arr[2][1]=arr[2][1]+1; arr[3][1]=arr[3][1]+1;
}
if (i1%9==0) {
arr[0][1]=0; arr[2][1]=0; arr[1][1]=0.5; arr[3][1]=0.5;
}
 
if ((i1%3 == 1) || (i1%3==2)) {
arr[0][2]=arr[0][2]+1; arr[1][2]=arr[1][2]+1; arr[2][2]=arr[2][2]+1; arr[3][2]=arr[3][2]+1;
}
if ((i1%9==0) || (i1%3 == 0)) {
arr[0][2]=0; arr[1][2]=0.5; arr[3][2]=0; arr[2][2]=0.5;
}
if (i1%9==0) { //arr[2][1]=0;
}
if (i1%3==0) {
//arr[0][2]=0;
//arr[3][2]=0;
}}
 
return rasst[27][12];
free (rasst);
}
 
int main(int argc, char *argv[]){
 
//РєРѕР»РёС‡РµСЃС‚РІРѕ Р°С‚РѕРјРѕРІ
int N = 108;
 
//РѕР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё 3*3*3
int Vol = 27; 
 
//РєРѕР»-РІРѕ С€Р°РіРѕРІ
int st = std::stoi(argv[1]);
 
//РїР°СЂР°РјРµС‚СЂС‹ РїРѕС‚РµРЅС†РёР°Р»РѕРІ
double a1; double r0; double a0; double p;
double q; double ksi; double B; double x;
double Rij; double alfa = std::stod(argv[2]);
double C11, C12, C44;
std::cout<< "Р’РІРµРґРёС‚Рµ Р·РЅР°С‡РµРЅРёРµ СЂР°СЃСЃС‚РѕСЏРЅРёСЏ alfa: = ";
std::cin>>alfa;
std::cout<< "Р’РІРµРґРёС‚Рµ РєРѕРѕСЂРґРёРЅР°С‚С‹ РЅР°С‡Р°Р»СЊРЅРѕР№ С‚РѕС‡РєРё:" <<endl;
double z1,z2,z3;
std::cout<<"X = ";
cin>>z1;
std::cout<<"Y = ";
cin>>z2;
std::cout<<"Z = ";
cin>>z3;
 
double D0[3][3];
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i==j) D0[i][j] = 1+alfa;
else D0[i][j] = 0;
}
}
 
double D1[3][3];
D1[0][0]=1+alfa;
D1[1][1]=1-alfa;
D1[2][2]=1;
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i!=j) D1[i][j]=0;
}}
 
double D2[3][3];
D2[0][0]=1+alfa; D2[1][1]=1-alfa; D2[2][2]=0;
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i!=j) D2[i][j]=0;
}}
 
double D3[3][3];
D3[0][0]=1; D3[1][1]=1; D3[2][2]=1/(1-(alfa*alfa));
D3[0][1]=alfa; D3[1][0]=D3[0][1]; D3[0][2]=0; D3[1][2]=0;
D3[2][0]=D3[0][2]; D3[2][1]=D3[1][2];
 
//Ag
double a0Ag = 4.085;
double EcAg=-2.96;
double BAg = 1.08;
double C11Ag = 1.32;
double C12Ag = 0.97;
double C44Ag = 0.51;
 
double ksiAg = 1.178;
double pAg = 10.928;
double qAg = 3.139;
 
 
double V=a0Ag*a0Ag*a0Ag;
std::cout<< "РћР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё V0 (Ag) = "<<V<<endl;
 
//Ni
double a0Ni = 3.523;
double EcNi = -4.435;
double EsolNi = 0.539;
double Ein_dimNi = 0.06;
double Eon_dimNi = -0.56;
double aNi = 0.0376;
double ksiNi = 1.07;
double pNi = 16.999;
double qNi = 1.189;
double C11Ni = 2.57;
double C12Ni = 1.55;
double C44Ni = 1.36;
 
V=a0Ni*a0Ni*a0Ni;
std::cout<< "РћР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё V0 (Ni) = "<<V<<endl;
 
//Ag и Ni
double aAgNi = (a0Ni + a0Ag) /2;
double C11AgNi = (C11Ag+C11Ni)/2;
double C12AgNi = (C12Ag+C12Ni)/2;
double C44AgNi = (C44Ag+C44Ni)/2;
 
 
 
int myrank, ranksize;
double Eo, Eplus;
 
cout<<"Р’РІРµРґРёС‚Рµ РЅР°С‡Р°Р»СЊРЅС‹Р№ РїСЂРµРґРµР»: ";
int k;
cin >> k;
cout<<"Р’РІРµРґРёС‚Рµ РєРѕРЅРµС‡РЅС‹Р№ РїСЂРµРґРµР»: ";
int l;
cin >> l;
 
std::cout<<"Р&#152;СЃРїРѕР»СЊР·СѓРµРјС‹Рµ СЌР»РµРјРµРЅС‚С‹: Ag Рё Ni."<<endl;
std::cout<<"Р¦РµР»СЊ: РѕРїСЂРµРґРµР»РёС‚СЊ РїРѕС‚РµРЅС†РёР°Р»С‹ РјРµР¶Р°С‚РѕРјРЅРѕРіРѕ РІР·Р°РёРјРѕРґРµР№СЃС‚РІРёСЏ РґР»СЏ СЃРёСЃС‚РµРјС‹ Ag/Ni (001)."<<endl;
std::cout<<"a(A) - РїР°СЂР°РјРµС‚СЂ СЂРµС€РµС‚РєРё (СЃ РёСЃРїРѕР»СЊР·РѕРІР°РЅРёРµРј СѓСЂР°РІРЅРµРЅРёСЏ СЃРѕСЃС‚РѕСЏРЅРёСЏ Р‘РµСЂС‡Р°-РњСѓСЂРЅР°РіР°РЅР°" << endl;
 
/*std::cout<<"Р’С‹С‡РёСЃР»РµРЅРёРµ СЌРЅРµСЂРіРёРё:"<<endl;
 
std::cout<<"РџРѕ Р°С‚РѕРјР°Рј k-l:"<<endl;
for (int i=k; i<=l; i++)
{
std::cout<<i<<") EiR = "<<ottenergy(a0Ag,EcAg,r0,pAg)<<endl;
std::cout<<i<<") EiB = "<<prenergy(qAg,ksiAg,r0,a0Ag)<<endl;
std::cout<<i<<") Ecoh = "<<fullenergy(ottenergy(a1,a0,r0,p), prenergy(q,ksi,Rij,r0))<<endl;
}*/
 
std::cout<<"AG"<<endl;
B = 4/(9*Vol) * (2*Eo+Eplus)/d2E(alfa)* 1.602;
C11 = (1/Vol) * (d2E(alfa)*C11 + d2E(alfa) *C12) * 1.602;
C12 = (1/Vol) * (d2E(alfa)*C11 - d2E(alfa)*C12) * 1.602;
C44 = (1/Vol) * (d2E(alfa)*C44) * 1.602;
 
 
cout<<"\n\n\n"<<BuildLattice(a0Ag);
 
 
//MPI_Finalize();
return 0;
 
}

Но при компиляции выводится следующая ошибка:
*** Error in './a.out': free(): invalid next size (fast): 0x000000000133cf20 ***
Ошибка сегментирования (сделан дамп памяти)

От чего зависит появление этой ошибки? И кто ранее с этим сталкивался?

@grgdvo · 05.03.2018, 13:18

строка 69: выход за пределы выделенного массива при обращении к элементу массива
строка 70: free не нужен, вы массива до этого объявляли непосредственно (с размещением в стеке), он удалится сам, по окончанию функции.

@Dionisius · 05.03.2018, 20:54 **[ТС]**

Сообщение от grgdvo

строка 69: выход за пределы выделенного массива при обращении к элементу массива
строка 70: free не нужен, вы массива до этого объявляли непосредственно (с размещением в стеке), он удалится сам, по окончанию функции.

Эти две строки я убрал, но теперь та же самая ошибка, но код у нее другой - 0x0000000001b6df20 , которая также связана с памятью.

@grgdvo · 06.03.2018, 06:42

строки 32 и 37: выход за пределы выделенного массива при обращении к элементу массива.

чтобы дальше не гадать, воспользуйтесь valgrind.
он покажет, где вы могли еще накосячить при работе с памятью.

@Dionisius · 06.03.2018, 08:29 **[ТС]**

Сообщение от grgdvo

чтобы дальше не гадать, воспользуйтесь valgrind.
он покажет, где вы могли еще накосячить при работе с памятью.

Уже установил этот valgrind, буду разбираться.

А пока вот то, на что ругается valgrind:

Bash
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
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
==23430== Memcheck, a memory error detector
==23430== Copyright (C) 2002-2015, and GNU GPL'd, by Julian Seward et al.
==23430== Using Valgrind-3.11.0 and LibVEX; rerun with -h for copyright info
==23430== Command: ./a.out 100 0.01
==23430== 
--23430-- Valgrind options:
--23430--    --leak-check=full
--23430--    -v
--23430-- Contents of /proc/version:
--23430--   Linux version 4.13.0-32-generic (buildd@lgw01-amd64-004) (gcc version 5.4.0 20160609 (Ubuntu 5.4.0-6ubuntu1~16.04.5)) #35~16.04.1-Ubuntu SMP Thu Jan 25 10:13:43 UTC 2018
--23430-- 
--23430-- Arch and hwcaps: AMD64, LittleEndian, amd64-sse3
--23430-- Page sizes: currently 4096, max supported 4096
--23430-- Valgrind library directory: /usr/lib/valgrind
--23430-- Reading syms from /home/den/a.out
--23430-- Reading syms from /lib/x86_64-linux-gnu/ld-2.23.so
--23430--   Considering /lib/x86_64-linux-gnu/ld-2.23.so ..
--23430--   .. CRC mismatch (computed aa979a42 wanted 9019bbb7)
--23430--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/ld-2.23.so ..
--23430--   .. CRC is valid
--23430-- Reading syms from /usr/lib/valgrind/memcheck-amd64-linux
--23430--   Considering /usr/lib/valgrind/memcheck-amd64-linux ..
--23430--   .. CRC mismatch (computed eea41ea9 wanted 2009db78)
--23430--    object doesn't have a symbol table
--23430--    object doesn't have a dynamic symbol table
--23430-- Scheduler: using generic scheduler lock implementation.
--23430-- Reading suppressions file: /usr/lib/valgrind/default.supp
==23430== embedded gdbserver: reading from /tmp/vgdb-pipe-from-vgdb-to-23430-by-den-on-???
==23430== embedded gdbserver: writing to   /tmp/vgdb-pipe-to-vgdb-from-23430-by-den-on-???
==23430== embedded gdbserver: shared mem   /tmp/vgdb-pipe-shared-mem-vgdb-23430-by-den-on-???
==23430== 
==23430== TO CONTROL THIS PROCESS USING vgdb (which you probably
==23430== don't want to do, unless you know exactly what you're doing,
==23430== or are doing some strange experiment):
==23430==   /usr/lib/valgrind/../../bin/vgdb --pid=23430 ...command...
==23430== 
==23430== TO DEBUG THIS PROCESS USING GDB: start GDB like this
==23430==   /path/to/gdb ./a.out
==23430== and then give GDB the following command
==23430==   target remote | /usr/lib/valgrind/../../bin/vgdb --pid=23430
==23430== --pid is optional if only one valgrind process is running
==23430== 
--23430-- REDIR: 0x401cfd0 (ld-linux-x86-64.so.2:strlen) redirected to 0x3809e181 (???)
--23430-- Reading syms from /usr/lib/valgrind/vgpreload_core-amd64-linux.so
--23430--   Considering /usr/lib/valgrind/vgpreload_core-amd64-linux.so ..
--23430--   .. CRC mismatch (computed 2567ccf6 wanted 49420590)
--23430--    object doesn't have a symbol table
--23430-- Reading syms from /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so
--23430--   Considering /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so ..
--23430--   .. CRC mismatch (computed 0e27c9a8 wanted ac585421)
--23430--    object doesn't have a symbol table
==23430== WARNING: new redirection conflicts with existing -- ignoring it
--23430--     old: 0x0401cfd0 (strlen              ) R-> (0000.0) 0x3809e181 ???
--23430--     new: 0x0401cfd0 (strlen              ) R-> (2007.0) 0x04c31020 strlen
--23430-- REDIR: 0x401b920 (ld-linux-x86-64.so.2:index) redirected to 0x4c30bc0 (index)
--23430-- REDIR: 0x401bb40 (ld-linux-x86-64.so.2:strcmp) redirected to 0x4c320d0 (strcmp)
--23430-- REDIR: 0x401dd30 (ld-linux-x86-64.so.2:mempcpy) redirected to 0x4c35270 (mempcpy)
--23430-- Reading syms from /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.21
--23430--   Considering /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.21 ..
--23430--   .. CRC mismatch (computed 9adedc07 wanted 319402e2)
--23430--    object doesn't have a symbol table
--23430-- Reading syms from /lib/x86_64-linux-gnu/libm-2.23.so
--23430--   Considering /lib/x86_64-linux-gnu/libm-2.23.so ..
--23430--   .. CRC mismatch (computed e8c3647b wanted c3efddac)
--23430--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/libm-2.23.so ..
--23430--   .. CRC is valid
--23430-- Reading syms from /lib/x86_64-linux-gnu/libgcc_s.so.1
--23430--   Considering /lib/x86_64-linux-gnu/libgcc_s.so.1 ..
--23430--   .. CRC mismatch (computed b9a68419 wanted 29d51b00)
--23430--    object doesn't have a symbol table
--23430-- Reading syms from /lib/x86_64-linux-gnu/libc-2.23.so
--23430--   Considering /lib/x86_64-linux-gnu/libc-2.23.so ..
--23430--   .. CRC mismatch (computed 7a8ee3e4 wanted a5190ac4)
--23430--   Considering /usr/lib/debug/lib/x86_64-linux-gnu/libc-2.23.so ..
--23430--   .. CRC is valid
--23430-- REDIR: 0x576aa00 (libc.so.6:strcasecmp) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x5766280 (libc.so.6:strcspn) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x576ccf0 (libc.so.6:strncasecmp) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x57686f0 (libc.so.6:strpbrk) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x5768a80 (libc.so.6:strspn) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x576a14b (libc.so.6:memcpy@GLIBC_2.2.5) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x5768400 (libc.so.6:rindex) redirected to 0x4c308a0 (rindex)
--23430-- REDIR: 0x575f130 (libc.so.6:malloc) redirected to 0x4c2db20 (malloc)
--23430-- REDIR: 0x5766720 (libc.so.6:strlen) redirected to 0x4c30f60 (strlen)
--23430-- REDIR: 0x5769bf0 (libc.so.6:__GI_memcmp) redirected to 0x4c33b90 (__GI_memcmp)
--23430-- REDIR: 0x5764cd0 (libc.so.6:strcmp) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x5764d10 (libc.so.6:__GI_strcmp) redirected to 0x4c31fe0 (__GI_strcmp)
--23430-- REDIR: 0x576f3f0 (libc.so.6:memcpy@@GLIBC_2.14) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
--23430-- REDIR: 0x577a820 (libc.so.6:__memcpy_sse2_unaligned) redirected to 0x4c324a0 (memcpy@@GLIBC_2.14)
--23430-- REDIR: 0x576f470 (libc.so.6:__GI_memcpy) redirected to 0x4c32b00 (__GI_memcpy)
Р’РІРµРґРёС‚Рµ Р·РЅР°С‡РµРЅРёРµ СЂР°СЃСЃС‚РѕСЏРЅРёСЏ alfa: = 0.01
--23430-- REDIR: 0x4ec7e60 (libstdc++.so.6:operator new(unsigned long)) redirected to 0x4c2e080 (operator new(unsigned long))
--23430-- REDIR: 0x4ec5f10 (libstdc++.so.6:operator delete(void*)) redirected to 0x4c2f1e0 (operator delete(void*))
--23430-- REDIR: 0x576a3b0 (libc.so.6:__GI_mempcpy) redirected to 0x4c34fa0 (__GI_mempcpy)
--23430-- REDIR: 0x5769bb0 (libc.so.6:bcmp) redirected to 0x4a286f0 (_vgnU_ifunc_wrapper)
 
--23430-- REDIR: 0x5771760 (libc.so.6:strchrnul) redirected to 0x4c34da0 (strchrnul)
--23430-- REDIR: 0x5769860 (libc.so.6:memchr) redirected to 0x4c32170 (memchr)
 
--23430-- REDIR: 0x5769ff0 (libc.so.6:__GI_memmove) redirected to 0x4c347e0 (__GI_memmove)
==23430== Invalid write of size 8
==23430==    at 0x40136F: BuildLattice(double) (in /home/den/a.out)
==23430==    by 0x401F71: ottenergy(double, double, double, double) (in /home/den/a.out)
==23430==    by 0x40298D: main (in /home/den/a.out)
==23430==  Address 0x5ab8358 is 0 bytes after a block of size 24 alloc'd
==23430==    at 0x4C2E0EF: operator new(unsigned long) (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==23430==    by 0x403CEF: __gnu_cxx::new_allocator<double>::allocate(unsigned long, void const*) (in /home/den/a.out)
==23430==    by 0x403B20: std::allocator_traits<std::allocator<double> >::allocate(std::allocator<double>&, unsigned long) (in /home/den/a.out)
==23430==    by 0x403989: std::_Vector_base<double, std::allocator<double> >::_M_allocate(unsigned long) (in /home/den/a.out)
==23430==    by 0x403744: std::_Vector_base<double, std::allocator<double> >::_M_create_storage(unsigned long) (in /home/den/a.out)
==23430==    by 0x403474: std::_Vector_base<double, std::allocator<double> >::_Vector_base(unsigned long, std::allocator<double> const&) (in /home/den/a.out)
==23430==    by 0x403EE3: std::vector<double, std::allocator<double> >::vector(std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403DCA: void std::_Construct<std::vector<double, std::allocator<double> >, std::vector<double, std::allocator<double> > const&>(std::vector<double, std::allocator<double> >*, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403C3A: std::vector<double, std::allocator<double> >* std::__uninitialized_fill_n<false>::__uninit_fill_n<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403AAF: std::vector<double, std::allocator<double> >* std::uninitialized_fill_n<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403911: std::vector<double, std::allocator<double> >* std::__uninitialized_fill_n_a<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> >, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&, std::allocator<std::vector<double, std::allocator<double> > >&) (in /home/den/a.out)
==23430==    by 0x403691: std::vector<std::vector<double, std::allocator<double> >, std::allocator<std::vector<double, std::allocator<double> > > >::_M_fill_initialize(unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
 
==23430== Invalid read of size 8
==23430==    at 0x401F72: ottenergy(double, double, double, double) (in /home/den/a.out)
==23430==    by 0xFFEFFFA00: ???
==23430==    by 0x60625F: ??? (in /home/den/a.out)
==23430==  Address 0xffffffffffffffd8 is not stack'd, malloc'd or (recently) free'd
==23430== 
==23430== 
==23430== Process terminating with default action of signal 11 (SIGSEGV)
==23430==  Access not within mapped region at address 0xFFFFFFFFFFFFFFD8
==23430==    at 0x401F72: ottenergy(double, double, double, double) (in /home/den/a.out)
==23430==    by 0xFFEFFFA00: ???
==23430==    by 0x60625F: ??? (in /home/den/a.out)
==23430==  If you believe this happened as a result of a stack
==23430==  overflow in your program's main thread (unlikely but
==23430==  possible), you can try to increase the size of the
==23430==  main thread stack using the --main-stacksize= flag.
==23430==  The main thread stack size used in this run was 8388608.
--23430-- REDIR: 0x575f4f0 (libc.so.6:free) redirected to 0x4c2ed80 (free)
==23430== 
==23430== HEAP SUMMARY:
==23430==     in use at exit: 72,704 bytes in 1 blocks
==23430==   total heap usage: 123 allocs, 122 frees, 80,404 bytes allocated
==23430== 
==23430== Searching for pointers to 1 not-freed blocks
==23430== Checked 107,800 bytes
==23430== 
==23430== LEAK SUMMARY:
==23430==    definitely lost: 0 bytes in 0 blocks
==23430==    indirectly lost: 0 bytes in 0 blocks
==23430==      possibly lost: 0 bytes in 0 blocks
==23430==    still reachable: 72,704 bytes in 1 blocks
==23430==         suppressed: 0 bytes in 0 blocks
==23430== Reachable blocks (those to which a pointer was found) are not shown.
==23430== To see them, rerun with: --leak-check=full --show-leak-kinds=all
==23430== 
==23430== ERROR SUMMARY: 4 errors from 2 contexts (suppressed: 0 from 0)
==23430== 
==23430== 1 errors in context 1 of 2:
==23430== Invalid read of size 8
==23430==    at 0x401F72: ottenergy(double, double, double, double) (in /home/den/a.out)
==23430==    by 0xFFEFFFA00: ???
==23430==    by 0x60625F: ??? (in /home/den/a.out)
==23430==  Address 0xffffffffffffffd8 is not stack'd, malloc'd or (recently) free'd
==23430== 
==23430== 
==23430== 3 errors in context 2 of 2:
==23430== Invalid write of size 8
==23430==    at 0x40136F: BuildLattice(double) (in /home/den/a.out)
==23430==    by 0x401F71: ottenergy(double, double, double, double) (in /home/den/a.out)
==23430==    by 0x40298D: main (in /home/den/a.out)
==23430==  Address 0x5ab8358 is 0 bytes after a block of size 24 alloc'd
==23430==    at 0x4C2E0EF: operator new(unsigned long) (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==23430==    by 0x403CEF: __gnu_cxx::new_allocator<double>::allocate(unsigned long, void const*) (in /home/den/a.out)
==23430==    by 0x403B20: std::allocator_traits<std::allocator<double> >::allocate(std::allocator<double>&, unsigned long) (in /home/den/a.out)
==23430==    by 0x403989: std::_Vector_base<double, std::allocator<double> >::_M_allocate(unsigned long) (in /home/den/a.out)
==23430==    by 0x403744: std::_Vector_base<double, std::allocator<double> >::_M_create_storage(unsigned long) (in /home/den/a.out)
==23430==    by 0x403474: std::_Vector_base<double, std::allocator<double> >::_Vector_base(unsigned long, std::allocator<double> const&) (in /home/den/a.out)
==23430==    by 0x403EE3: std::vector<double, std::allocator<double> >::vector(std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403DCA: void std::_Construct<std::vector<double, std::allocator<double> >, std::vector<double, std::allocator<double> > const&>(std::vector<double, std::allocator<double> >*, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403C3A: std::vector<double, std::allocator<double> >* std::__uninitialized_fill_n<false>::__uninit_fill_n<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403AAF: std::vector<double, std::allocator<double> >* std::uninitialized_fill_n<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430==    by 0x403911: std::vector<double, std::allocator<double> >* std::__uninitialized_fill_n_a<std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> >, std::vector<double, std::allocator<double> > >(std::vector<double, std::allocator<double> >*, unsigned long, std::vector<double, std::allocator<double> > const&, std::allocator<std::vector<double, std::allocator<double> > >&) (in /home/den/a.out)
==23430==    by 0x403691: std::vector<std::vector<double, std::allocator<double> >, std::allocator<std::vector<double, std::allocator<double> > > >::_M_fill_initialize(unsigned long, std::vector<double, std::allocator<double> > const&) (in /home/den/a.out)
==23430== 
==23430== ERROR SUMMARY: 4 errors from 2 contexts (suppressed: 0 from 0)

@grgdvo · 06.03.2018, 15:01

1) компилируйте с отладочной информацией, тогда valgrind должен показывать функции и номера строк, где происходит ошибка
2) на всякий случай можно принудительно компилятору отключить всякую оптимизацию. пусть код будет тупой, но зато найдете ошибку
3) еще ошибка в строке 27. цикл по i пробегает от 0 до 3 => 4 значения (0,1,2,3), а массив вы выделили всего 3 ячейки. налицо выход за пределы массива с попаданием в чужую память, как итог SIGSEGV.

@Dionisius · 06.03.2018, 20:16 **[ТС]**

Сообщение от grgdvo

компилируйте с отладочной информацией, тогда valgrind должен показывать функции и номера строк, где происходит ошибка

Я так понимаю, что команды valgrind --leak-check=full --show-leak-kinds=all -v ./a.out 100 0.01 будет недостаточно?

Сообщение от grgdvo

на всякий случай можно принудительно компилятору отключить всякую оптимизацию.

Каким образом можно это сделать?

@grgdvo · 07.03.2018, 05:58

Bash
1
gcc -O0 -g my.cpp

@Dionisius · 08.03.2018, 21:52 **[ТС]**

delete, т.к. это дубль ответа ниже

@Dionisius · 08.03.2018, 21:52 **[ТС]**

Сообщение от grgdvo

еще ошибка в строке 27. цикл по i пробегает от 0 до 3 => 4 значения (0,1,2,3), а массив вы выделили всего 3 ячейки. налицо выход за пределы массива с попаданием в чужую память, как итог SIGSEGV.

Эту две строки

C++
1
2
if (i==0) arr[i][j]=0;
arr[i+1][i]=0;

я заменил на эту:

C++
1
arr[i][j]=0;

И вместо 4 ошибок valgrind выдаёт одну:

Bash
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
==26764== Invalid read of size 8
==26764==    at 0x401F21: ottenergy(double, double, double, double) (kursach.cpp:150)
==26764==    by 0xFFEFFFA00: ???
==26764==    by 0x60625F: ??? (in /home/den/a.out)
==26764==  Address 0xffffffffffffffd8 is not stack'd, malloc'd or (recently) free'd
==26764== 
==26764== 
==26764== Process terminating with default action of signal 11 (SIGSEGV)
==26764==  Access not within mapped region at address 0xFFFFFFFFFFFFFFD8
==26764==    at 0x401F21: ottenergy(double, double, double, double) (kursach.cpp:150)
==26764==    by 0xFFEFFFA00: ???
==26764==    by 0x60625F: ??? (in /home/den/a.out)
==26764==  If you believe this happened as a result of a stack
==26764==  overflow in your program's main thread (unlikely but
==26764==  possible), you can try to increase the size of the
==26764==  main thread stack using the --main-stacksize= flag.
==26764==  The main thread stack size used in this run was 8388608.
--26764-- REDIR: 0x575f4f0 (libc.so.6:free) redirected to 0x4c2ed80 (free)
==26764== 
==26764== HEAP SUMMARY:
==26764==     in use at exit: 72,704 bytes in 1 blocks
==26764==   total heap usage: 123 allocs, 122 frees, 80,404 bytes allocated
==26764== 
==26764== Searching for pointers to 1 not-freed blocks
==26764== Checked 107,808 bytes
==26764== 
==26764== 72,704 bytes in 1 blocks are still reachable in loss record 1 of 1
==26764==    at 0x4C2DB8F: malloc (in /usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==26764==    by 0x4EC3EFF: ??? (in /usr/lib/x86_64-linux-gnu/libstdc++.so.6.0.21)
==26764==    by 0x40106B9: call_init.part.0 (dl-init.c:72)
==26764==    by 0x40107CA: call_init (dl-init.c:30)
==26764==    by 0x40107CA: _dl_init (dl-init.c:120)
==26764==    by 0x4000C69: ??? (in /lib/x86_64-linux-gnu/ld-2.23.so)
==26764==    by 0x2: ???
==26764==    by 0xFFF000232: ???
==26764==    by 0xFFF00023A: ???
==26764==    by 0xFFF00023E: ???
==26764== 
==26764== LEAK SUMMARY:
==26764==    definitely lost: 0 bytes in 0 blocks
==26764==    indirectly lost: 0 bytes in 0 blocks
==26764==      possibly lost: 0 bytes in 0 blocks
==26764==    still reachable: 72,704 bytes in 1 blocks
==26764==         suppressed: 0 bytes in 0 blocks
==26764== 
==26764== ERROR SUMMARY: 1 errors from 1 contexts (suppressed: 0 from 0)
==26764== 
==26764== 1 errors in context 1 of 1:
==26764== Invalid read of size 8
==26764==    at 0x401F21: ottenergy(double, double, double, double) (kursach.cpp:150)
==26764==    by 0xFFEFFFA00: ???
==26764==    by 0x60625F: ??? (in /home/den/a.out)
==26764==  Address 0xffffffffffffffd8 is not stack'd, malloc'd or (recently) free'd
==26764== 
==26764== ERROR SUMMARY: 1 errors from 1 contexts (suppressed: 0 from 0)

И теперь исправленный код:

Кликните здесь для просмотра всего текста

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
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
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <random>
#include <tbb/tbb.h>
#include <tuple>
#include <cstring>
#include <iomanip>
//#include "mpi.h"
 
using namespace std;
 
 
//СЃС‚СЂРѕРёС‚РµР»СЊСЃС‚РІРѕ СЂРµС€РµС‚РєРё
double BuildLattice(double a0)
{
int x=-1;
int y=-1;
int z=-1;
double Rij=0;
 
std::vector<std::vector<double>>arr(108, vector<double>(3));
//РјР°СЃСЃРёРІ СЂР°СЃСЃС‚РѕСЏРЅРёР№ РјРµР¶РґСѓ Р°С‚РѕРјР°РјРё
double rasst[27][12];
 
for (int i=0; i<4; i++)
{
for (int j=0; j<3; j++)
{
arr[i][j]=0;
}
}
arr[1][1]=0.5; 
arr[1][2]=0.5;
arr[2][0]=0.5;
arr[2][2]=0.5;
arr[3][0]=0.5;
arr[3][1]=0.5;
 
cout<<endl;
std::vector<std::vector<double>>vct(4, vector<double>(3));
 
int i4;
for (int i1=1; i1<=27; i1++)
{
i4=0;
for (int i2=0; i2<4; i2++)
{
for (int i3=0; i3<3; i3++)
{
std::cout<<arr[i2][i3]<<" ";
if(i2==i3)
rasst[i1][i4]=0;
if ((i2-i3 ==1) || (i3-i2 == 1)) 
rasst[i2][i4]=1;
if ((i2-i3 ==2) || (i3-i2 == 2))
rasst[i2][i4]=2;
if ((i2-i3 ==3) || (i3-i2 == 3))
rasst[i2][i4]=3;
 
i4++;
 
//msv[27*i1+4*i2+3*i3]=rasst[i2][i3];
//cout<<"msv["<<27*i1+4*i2+3*i3<<"] = "<< msv[27*i1+4*i2+3*i3] <<endl;
}
std::cout<<endl;
}
cout<<endl;
if ((i1==9) || (i1==18)) 
{
 
arr[0][0]=arr[0][0]+1;
arr[1][0]=arr[1][0]+1;
arr[2][0]=arr[2][0]+1;
arr[3][0]=arr[3][0]+1;
//vct[0][0]=arr[0][0];
//vct[1][0]=arr[1][0];
//vct[2][0]=arr[2][0];
//vct[3][0]=arr[3][0];
}
if ((i1%9==3) || (i1%9 == 6)) 
{
arr[0][1]=arr[0][1]+1;
arr[1][1]=arr[1][1]+1;
arr[2][1]=arr[2][1]+1;
arr[3][1]=arr[3][1]+1;
//vct[0][1]=arr[0][1];
//vct[1][1]=arr[1][1];
//vct[2][1]=arr[2][1];
//vct[3][1]=arr[3][1];
}
if (i1%9==0) 
{
arr[0][1]=0;
arr[2][1]=0;
arr[1][1]=0.5;
arr[3][1]=0.5;
}
 
if ((i1%3 == 1) || (i1%3==2)) 
{
arr[0][2]=arr[0][2]+1;
arr[1][2]=arr[1][2]+1;
arr[2][2]=arr[2][2]+1;
arr[3][2]=arr[3][2]+1;
//vct[0][2]=arr[0][2];
//vct[1][2]=arr[1][2];
//vct[2][2]=arr[2][2];
//vct[3][2]=arr[3][2];
}
//if (i1%9==1) 
//arr[0][2]=0;
 
if ((i1%9==0) || (i1%3 == 0)) 
{
arr[0][2]=0;
arr[1][2]=0.5;
arr[3][2]=0;
arr[2][2]=0.5;
}
if (i1%9==0) 
{
//arr[2][1]=0;
}
/*if ((i1%3 == 1) || (i1%3==2)) 
{
arr[0][2]=arr[0][2]+1;
arr[1][2]=arr[1][2]+1;
arr[2][2]=arr[2][2]+1;
arr[3][2]=arr[3][2]+1;
 
//if ((arr[0][2] ==3)||(arr[1][2]==3)||(arr[2][2]==3
}*/
if (i1%3==0) {
//arr[0][2]=0;
//arr[3][2]=0;
}
}
 
 
}
 
//СЌРЅРµСЂРіРёСЏ РѕС‚С‚Р°Р»РєРёРІР°РЅРёСЏ
double ottenergy(double a1, double a0, double r0, double p)
{
//std::vector<double>EiR(324);
double EiR;
//Rij=msv[i]
EiR = (a1*(BuildLattice(a0)-r0)+a0)* exp(-p*((BuildLattice(a0)/r0)-1));
 
return EiR;
 
}
 
//СЌРЅРµСЂРіРёСЏ РїСЂРёС‚СЏР¶РµРЅРёСЏ
double prenergy(double q, double ksi, double r0, double a0){
double EiB, EiB1;
 
EiB1 = ksi*ksi*exp(-2*q*(BuildLattice(a0)/r0)-1);
 
EiB = -sqrt(EiB1);
return EiB;
}
 
//РїРѕР»РЅР°СЏ СЌРЅРµСЂРіРёСЏ СЃРёСЃС‚РµРјС‹
double fullenergy(double EiB, double EiR){
double Ecoh;
Ecoh = (EiR+EiB);
 
return Ecoh;
}
 
 
double d2E (double alfa)
{
double Ecoh;
double EiB;
double EiR;
double D2Ealfa = fullenergy(EiB, EiR) - 2 * fullenergy(EiB, EiR)+ fullenergy(alfa,EiR);
return D2Ealfa;
}
 
//РјРµС‚РѕРґ СЃР»СѓС‡Р°Р№РЅРѕРіРѕ РїРѕРёСЃРєР°
void randomsearch(double z1, double z2, double z3)
{
 
}
 
//РІС‹С‡РёСЃР»РµРЅРёРµ СЌРЅРµСЂРіРёРё СЂР°СЃС‚РІРѕСЂРёРјРѕСЃС‚Рё
void rastv()
{
 
}
 
//СЌРЅРµСЂРіРёСЏ СЃРІСЏР·Рё РґРёРјРµСЂР° РІ РїРѕРІРµСЂС…РЅРѕСЃС‚РЅРѕРј СЃР»РѕРµ
void povsl_dimer()
{
 
}
 
 
//СЌРЅРµСЂРіРёСЏ СЃРІСЏР·Рё РґРёРјРµСЂР° РЅР° РїРѕРІРµСЂС…РЅРѕСЃС‚Рё
void povdimer()
{
 
}
 
int main(int argc, char *argv[]){
//int myrank; 
//int rangsize;
//MPI_Init (&argc, &argv);
//MPI_Status status;
//int size = mpi::comm_size();
//int rank = mpi::comm_rank();
//MPI_Comm_rank(MPI_COMM_WORLD, &myrank);
//MPI_Comm_size(MPI_COMM_WORLD, &ranksize);
 
//РєРѕР»РёС‡РµСЃС‚РІРѕ Р°С‚РѕРјРѕРІ
int N = 108;
 
//РѕР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё 3*3*3
int Vol = 27; 
 
//РєРѕР»-РІРѕ С€Р°РіРѕРІ
int st = std::stoi(argv[1]);
 
//РїР°СЂР°РјРµС‚СЂС‹ РїРѕС‚РµРЅС†РёР°Р»РѕРІ
double a1 = 1;
double r0;
double a0;
double p;
double q;
double ksi;
double B;
double x;
double Rij;
double alfa = std::stod(argv[2]);
double C11, C12, C44;
std::cout<< "Р’РІРµРґРёС‚Рµ Р·РЅР°С‡РµРЅРёРµ СЂР°СЃСЃС‚РѕСЏРЅРёСЏ alfa: = ";
std::cin>>alfa;
 
std::cout<< "Р’РІРµРґРёС‚Рµ РєРѕРѕСЂРґРёРЅР°С‚С‹ РЅР°С‡Р°Р»СЊРЅРѕР№ С‚РѕС‡РєРё:" <<endl;
 
double z1,z2,z3;
std::cout<<"X = ";
cin>>z1;
std::cout<<"Y = ";
cin>>z2;
std::cout<<"Z = ";
cin>>z3;
 
double D0[3][3];
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i==j) D0[i][j] = 1+alfa;
else D0[i][j] = 0;
}
}
 
double D1[3][3];
D1[0][0]=1+alfa;
D1[1][1]=1-alfa;
D1[2][2]=1;
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i!=j) D1[i][j]=0;
}
}
 
double D2[3][3];
D2[0][0]=1+alfa;
D2[1][1]=1-alfa;
D2[2][2]=0;
for (int i=0; i<3; i++)
{
for (int j=0; j<3; j++)
{
if (i!=j) D2[i][j]=0;
}
}
 
double D3[3][3];
D3[0][0]=1;
D3[1][1]=1;
D3[2][2]=1/(1-(alfa*alfa));
D3[0][1]=alfa;
D3[1][0]=D3[0][1];
D3[0][2]=0;
D3[1][2]=0;
D3[2][0]=D3[0][2];
D3[2][1]=D3[1][2];
 
 
//Ag
double a0Ag = 4.085;
double EcAg=-2.96;
double BAg = 1.08;
double C11Ag = 1.32;
double C12Ag = 0.97;
double C44Ag = 0.51;
 
double ksiAg = 1.178;
double pAg = 10.928;
double qAg = 3.139;
 
 
double V=a0Ag*a0Ag*a0Ag;
std::cout<< "РћР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё V0 (Ag) = "<<V<<endl;
 
//Ni
double a0Ni = 3.523;
double EcNi = -4.435;
double EsolNi = 0.539;
double Ein_dimNi = 0.06;
double Eon_dimNi = -0.56;
double aNi = 0.0376;
double ksiNi = 1.07;
double pNi = 16.999;
double qNi = 1.189;
double C11Ni = 2.57;
double C12Ni = 1.55;
double C44Ni = 1.36;
 
V=a0Ni*a0Ni*a0Ni;
std::cout<< "РћР±СЉС‘Рј РєСЂРёСЃС‚Р°Р»Р»РёС‡РµСЃРєРѕР№ СЂРµС€РµС‚РєРё V0 (Ni) = "<<V<<endl;
 
//Ag Рё Ni
double aAgNi = (a0Ni + a0Ag) /2;
double C11AgNi = (C11Ag+C11Ni)/2;
double C12AgNi = (C12Ag+C12Ni)/2;
double C44AgNi = (C44Ag+C44Ni)/2;
 
 
 
int myrank, ranksize;
double Eo, Eplus;
 
cout<<"Р’РІРµРґРёС‚Рµ РЅР°С‡Р°Р»СЊРЅС‹Р№ РїСЂРµРґРµР»: ";
int k;
cin >> k;
cout<<"Р’РІРµРґРёС‚Рµ РєРѕРЅРµС‡РЅС‹Р№ РїСЂРµРґРµР»: ";
int l;
cin >> l;
 
std::cout<<"Р?СЃРїРѕР»СЊР·СѓРµРјС‹Рµ СЌР»РµРјРµРЅС‚С‹: Ag Рё Ni."<<endl;
std::cout<<"Р¦РµР»СЊ: РѕРїСЂРµРґРµР»РёС‚СЊ РїРѕС‚РµРЅС†РёР°Р»С‹ РјРµР¶Р°С‚РѕРјРЅРѕРіРѕ РІР·Р°РёРјРѕРґРµР№СЃС‚РІРёСЏ РґР»СЏ СЃРёСЃС‚РµРјС‹ Ag/Ni (001)."<<endl;
std::cout<<"a(A) - РїР°СЂР°РјРµС‚СЂ СЂРµС€РµС‚РєРё (СЃ РёСЃРїРѕР»СЊР·РѕРІР°РЅРёРµРј СѓСЂР°РІРЅРµРЅРёСЏ СЃРѕСЃС‚РѕСЏРЅРёСЏ Р‘РµСЂС‡Р°-РњСѓСЂРЅР°РіР°РЅР°" << endl;
 
std::cout<<"Р’С‹С‡РёСЃР»РµРЅРёРµ СЌРЅРµСЂРіРёРё:"<<endl;
/*for (int i=0; i<324; i++)
{
std::cout<<i<<") EiR = "<<ottenergy(a0Ag,EcAg,r0,pAg)<<endl;
//std::cout<<i<<") EiB = "<<prenergy(qAg,ksiAg,Rij,r0)<<endl;
//std::cout<<i<<") Ecoh = "<<fullenergy(ottenergy(a1,a0,r0,p), prenergy(q,ksi,Rij,r0))<<endl;
 
}
 
for (int i=0; i<324; i++)
{
//std::cout<<i<<") EiR = "<<ottenergy(aAg,EcAg,r0,pAg)<<endl;
std::cout<<i<<") EiB = "<<prenergy(qAg,ksiAg,r0,a0Ag)<<endl;
//std::cout<<i<<") Ecoh = "<<fullenergy(ottenergy(a1,a0,r0,p), prenergy(q,ksi,Rij,r0))<<endl;
 
}*/
 
std::cout<<"РџРѕ Р°С‚РѕРјР°Рј k-l:"<<endl;
for (int i=k; i<=l; i++)
{
std::cout<<i<<") EiR = "<<ottenergy(a0Ag,EcAg,r0,pAg)<<endl;
std::cout<<i<<") EiB = "<<prenergy(qAg,ksiAg,r0,a0Ag)<<endl;
std::cout<<i<<") Ecoh = "<<fullenergy(ottenergy(a1,a0,r0,p), prenergy(q,ksi,Rij,r0))<<endl;
}
 
std::cout<<"AG"<<endl;
B = 4/(9*Vol) * (2*Eo+Eplus)/d2E(alfa)* 1.602;
C11 = (1/Vol) * (d2E(alfa)*C11 + d2E(alfa) *C12) * 1.602;
C12 = (1/Vol) * (d2E(alfa)*C11 - d2E(alfa)*C12) * 1.602;
C44 = (1/Vol) * (d2E(alfa)*C44) * 1.602;
 
//if (C11!= C11Ag) cout<<"\nC11 РЅРµ СЃС…РѕРґРёС‚СЃСЏ!\n"<<endl;
 
//if (C12!= C12Ag) cout<<"\nC12 РЅРµ СЃС…РѕРґРёС‚СЃСЏ!\n"<<endl;
 
//if (C44!= C44Ag) cout<<"\nC44 РЅРµ СЃС…РѕРґРёС‚СЃСЏ!\n"<<endl;
 
for (int i=0; i<324; i++)
{
std::cout<<i<<") EiR = "<<ottenergy(a0Ag,EcAg,r0,pAg)<<endl;
}
 
for (int i=0; i<324; i++)
{
std::cout<<i<<") EiB = "<<prenergy(qNi,ksiNi,r0,a0Ni)<<endl;
 
}
cout<<"\n\n\n"<<BuildLattice(a0Ag);
 
//MPI_Finalize();
return 0;
 
}

Новые блоги и статьи Все статьи Все блоги /
Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это дополнительная запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .

@grgdvo 1264 / 978 / 384 Регистрация: 02.09.2012 Сообщений: 3,024
	05.03.2018, 13:18
	строка 69: выход за пределы выделенного массива при обращении к элементу массива строка 70: free не нужен, вы массива до этого объявляли непосредственно (с размещением в стеке), он удалится сам, по окончанию функции. 1

@grgdvo 1264 / 978 / 384 Регистрация: 02.09.2012 Сообщений: 3,024
	06.03.2018, 06:42
	строки 32 и 37: выход за пределы выделенного массива при обращении к элементу массива. чтобы дальше не гадать, воспользуйтесь valgrind. он покажет, где вы могли еще накосячить при работе с памятью. 0

@grgdvo 1264 / 978 / 384 Регистрация: 02.09.2012 Сообщений: 3,024
	06.03.2018, 15:01
	1) компилируйте с отладочной информацией, тогда valgrind должен показывать функции и номера строк, где происходит ошибка 2) на всякий случай можно принудительно компилятору отключить всякую оптимизацию. пусть код будет тупой, но зато найдете ошибку 3) еще ошибка в строке 27. цикл по i пробегает от 0 до 3 => 4 значения (0,1,2,3), а массив вы выделили всего 3 ячейки. налицо выход за пределы массива с попаданием в чужую память, как итог SIGSEGV. 0

@Dionisius 2 / 2 / 0 Регистрация: 19.10.2013 Сообщений: 159
	08.03.2018, 21:52 [ТС]
	delete, т.к. это дубль ответа ниже 0