21 out <<
"TpcSpaceChargeMatrixContainerv1" << std::endl;
22 out <<
" phibins: " <<
m_phibins << std::endl;
23 out <<
" rbins: " <<
m_rbins << std::endl;
24 out <<
" zbins: " <<
m_zbins << std::endl;
42 if( iphi < 0 || iphi >=
m_phibins )
return -1;
43 if( ir < 0 || ir >=
m_rbins )
return -1;
44 if( iz < 0 || iz >=
m_zbins )
return -1;
69 return m_rhs[cell_index][
i];
79 m_entries = std::vector<int>( totalbins, 0 );
80 m_lhs = std::vector<matrix_t>( totalbins, {{}} );
81 m_rhs = std::vector<column_t>( totalbins, {{}} );
124 std::cout <<
"TpcSpaceChargeMatrixContainerv1::add - inconsistent grid sizes" << std::endl;
129 for(
size_t cell_index = 0; cell_index <
m_lhs.size(); ++cell_index )
133 for(
size_t cell_index = 0; cell_index <
m_lhs.size(); ++cell_index )
139 for(
size_t cell_index = 0; cell_index <
m_lhs.size(); ++cell_index )
149 if( cell_index < 0 || cell_index >= (
int)
m_rhs.size() )
return false;
156 if( cell_index < 0 || cell_index >= (
int)
m_rhs.size() )
return false;
157 if( i < 0 || i >=
m_ncoord )
return false;
164 if( cell_index < 0 || cell_index >= (
int)
m_lhs.size() )
return false;
165 if( i < 0 || i >=
m_ncoord )
return false;
166 if( j < 0 || j >=
m_ncoord )
return false;