This commit is contained in:
James E McClure 2018-06-05 20:18:48 -04:00
parent 29eb90084f
commit 7debe96512

View File

@ -1064,10 +1064,11 @@ Array<TYPE, FUN> Array<TYPE, FUN>::reverseDim() const
template<class TYPE, class FUN>
Array<TYPE, FUN> Array<TYPE, FUN>::coarsen( const Array<TYPE, FUN> &filter ) const
{
size_t S2 = size();
auto S2 = size();
for ( size_t i = 0; i < S2.size(); i++ ) {
S2[i] /= filter.size( i );
if ( S2[i] * filter.size( i ) != size( i ) )
size_t S2_i=S2[i];
S2_i /= filter.size( i );
if ( S2_i * filter.size( i ) != size( i ) )
throw std::invalid_argument( "Array must be multiple of filter size" );
}
Array<TYPE, FUN> y( S2 );