diff options
Diffstat (limited to 'include/phat/helpers/dualize.h')
| -rw-r--r-- | include/phat/helpers/dualize.h | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/include/phat/helpers/dualize.h b/include/phat/helpers/dualize.h index 752103d..5f9a05c 100644 --- a/include/phat/helpers/dualize.h +++ b/include/phat/helpers/dualize.h @@ -31,11 +31,24 @@ namespace phat { index nr_of_columns = boundary_matrix.get_num_cols();
dual_matrix.resize( nr_of_columns );
dual_dims.resize( nr_of_columns );
+
+ std::vector< dimension >& dual_sizes = dual_dims;
+
column temp_col;
for( index cur_col = 0; cur_col < nr_of_columns; cur_col++ ) {
boundary_matrix.get_col( cur_col, temp_col );
for( index idx = 0; idx < (index)temp_col.size(); idx++)
+ dual_sizes[ nr_of_columns - 1 - temp_col[ idx ] ]++;
+ }
+
+ for( index cur_col = 0; cur_col < nr_of_columns; cur_col++ ) {
+ dual_matrix[cur_col].reserve(dual_sizes[cur_col]);
+ }
+
+ for( index cur_col = 0; cur_col < nr_of_columns; cur_col++ ) {
+ boundary_matrix.get_col( cur_col, temp_col );
+ for( index idx = 0; idx < (index)temp_col.size(); idx++)
dual_matrix[ nr_of_columns - 1 - temp_col[ idx ] ].push_back( nr_of_columns - 1 - cur_col );
}
|