@@ -89,19 +89,19 @@ void ReportEvent::lfp_calc(NrnThread* nt) {
89
89
int count = 0 ;
90
90
double sum = 0.0 ;
91
91
for (const auto & kv: cell_mapping->lfp_factors ) {
92
- int segment_id = kv.first ;
93
- double factor = kv.second ;
94
- if (std::isnan (factor)) {
95
- factor = 0.0 ;
96
- }
97
- double iclamp = 0.0 ;
98
- for (const auto & value: summation_report.currents_ [segment_id]) {
99
- double current_value = *value.first ;
100
- int scale = value.second ;
101
- iclamp += current_value * scale;
102
- }
103
- sum += (fast_imem_rhs[segment_id] + iclamp) * factor;
104
- count++;
92
+ int segment_id = kv.first ;
93
+ double factor = kv.second ;
94
+ if (std::isnan (factor)) {
95
+ factor = 0.0 ;
96
+ }
97
+ double iclamp = 0.0 ;
98
+ for (const auto & value: summation_report.currents_ [segment_id]) {
99
+ double current_value = *value.first ;
100
+ int scale = value.second ;
101
+ iclamp += current_value * scale;
102
+ }
103
+ sum += (fast_imem_rhs[segment_id] + iclamp) * factor;
104
+ count++;
105
105
}
106
106
*(to_report.front ().var_value ) = sum;
107
107
}
@@ -115,8 +115,7 @@ void ReportEvent::deliver(double t, NetCvode* nc, NrnThread* nt) {
115
115
{
116
116
if (report_type == ReportType::SummationReport) {
117
117
summation_alu (nt);
118
- }
119
- else if (report_type == ReportType::LFPReport) {
118
+ } else if (report_type == ReportType::LFPReport) {
120
119
lfp_calc (nt);
121
120
}
122
121
// each thread needs to know its own step
0 commit comments