@@ -11,10 +11,6 @@ export SparseMatrixAssembler
11
11
export assemble
12
12
export assemble!
13
13
14
- # @fverdugo for the moment the abstract interface of Assembler
15
- # (and therefore its concrete implementations)
16
- # assumes single field, and single term
17
-
18
14
"""
19
15
Abstract assembly operator
20
16
Parametrized by the type of returned matrix and vector
@@ -122,33 +118,4 @@ function assemble!(
122
118
mat. nzval .= m. nzval
123
119
end
124
120
125
- # Draft of multi field assembler
126
- # function _assemble_sparse_matrix_values(mf_vals,mf_rows,mf_cols,I,E)
127
- # aux_row = I[]; aux_col = I[]; aux_val = E[]
128
- # for (mf_rows_c, mf_cols_c, mf_vals_c) in zip(mf_rows,mf_cols,mf_vals)
129
- # for (vals_c, (ifield, jfield)) in eachblock(mf_vals_c)
130
- # rows_c = mf_rows_c[ifield]
131
- # cols_c = mf_cols_c[jfield]
132
- # row_offset = row_offsets[ifield]
133
- # col_offset = col_offsets[jfield]
134
- # _asseble_cell_values!(aux_row,aux_col,aux_val,rows_c,cols_c,vals_c,col_offset,row_offset)
135
- # end
136
- # end
137
- # (aux_row, aux_col, aux_val)
138
- # end
139
- #
140
- # function _asseble_cell_values!(aux_row,aux_col,aux_val,rows_c,cols_c,vals_c,col_offset,row_offset)
141
- # for (j,gidcol) in enumerate(cols_c)
142
- # if gidcol > 0
143
- # for (i,gidrow) in enumerate(rows_c)
144
- # if gidrow > 0
145
- # push!(aux_row, gidrow+row_offset)
146
- # push!(aux_col, gidcol+col_offset)
147
- # push!(aux_val, vals_c[i,j])
148
- # end
149
- # end
150
- # end
151
- # end
152
- # end
153
-
154
121
end # module Assemblers
0 commit comments