@@ -89,19 +89,19 @@ void ReportEvent::lfp_calc(NrnThread* nt) {
8989 int count = 0 ;
9090 double sum = 0.0 ;
9191 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++;
105105 }
106106 *(to_report.front ().var_value ) = sum;
107107 }
@@ -115,8 +115,7 @@ void ReportEvent::deliver(double t, NetCvode* nc, NrnThread* nt) {
115115 {
116116 if (report_type == ReportType::SummationReport) {
117117 summation_alu (nt);
118- }
119- else if (report_type == ReportType::LFPReport) {
118+ } else if (report_type == ReportType::LFPReport) {
120119 lfp_calc (nt);
121120 }
122121 // each thread needs to know its own step
0 commit comments