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
|
#define CATCH_CONFIG_RUNNER
#include <catch2/catch.hpp>
#include <ocfunc.h>
#include <code.h>
#include <section.h>
#include <neuron.h>
extern void nrn_threads_create(int, int);
extern void nrn_threads_free();
extern int ivocmain_session(int, const char**, const char**, int);
extern int nrn_main_launch;
extern int nrn_nobanner_;
/// Needed for compilation
void modl_reg() {}
extern int nrn_nthread;
extern NrnThread* nrn_threads;
extern int nrn_use_fast_imem;
extern int use_cachevec;
int main(int argc, char* argv[]) {
// global setup...
nrn_main_launch = 2;
int argc_nompi = 2;
const char* argv_nompi[] = {"NEURON", "-nogui", nullptr};
nrn_nobanner_ = 1;
ivocmain_session(argc_nompi, argv_nompi, NULL, 0);
#undef run
int result = Catch::Session().run(argc, argv);
#define run hoc_run
// global clean-up...
return result;
}
SCENARIO("Test fast_imem calculation", "[Neuron][fast_imem]") {
GIVEN("A section") {
hoc_oc("create s\n");
WHEN("fast_imem and cachevec is allocated") {
nrn_use_fast_imem = true;
use_cachevec = 1;
nrn_fast_imem_alloc();
THEN("nrn_fast_imem should not be nullptr") {
for (int it = 0; it < nrn_nthread; ++it) {
NrnThread* nt = &nrn_threads[it];
REQUIRE(nt->_nrn_fast_imem != nullptr);
}
}
}
WHEN("fast_imem is created") {
hoc_oc(
"objref cvode\n"
"cvode = new CVode()\n"
"cvode.use_fast_imem(1)\n");
WHEN("iinitialize and run nrn_calc_fast_imem") {
hoc_oc("finitialize(-65)\n");
for (NrnThread* nt = nrn_threads; nt < nrn_threads + nrn_nthread; ++nt) {
nrn_calc_fast_imem(nt);
}
THEN("The current in this section is 0") {
for (NrnThread* nt = nrn_threads; nt < nrn_threads + nrn_nthread; ++nt) {
for (int i = 0; i < nt->end; ++i) {
REQUIRE(nt->_nrn_fast_imem->_nrn_sav_rhs[i] == 0.0);
}
}
}
}
}
hoc_oc("delete_section()");
}
}
TEST_CASE("Test return code of execerror", "[NEURON][execerror]") {
REQUIRE(hoc_oc("execerror(\"test error\")") > 0);
}
// AddressSanitizer seems to intercept the mallinfo[2]() system calls and return
// null values from them.
#ifndef NRN_ASAN_ENABLED
TEST_CASE("Test nrn_mallinfo returns non-zero", "[NEURON][nrn_mallinfo]") {
SECTION("HOC") {
REQUIRE(
hoc_oc(
"if (nrn_mallinfo(0) <= 0) { execerror(\"nrn_mallinfo returned <= 0 in HOC\") }") ==
0);
}
SECTION("C++") {
REQUIRE(nrn_mallinfo(0) > 0);
}
}
#endif
|
