Ex252 should run again

Author: chmerdon

src/common_restrictions/boundarydata_restriction.jl

@@ -47,7 +47,7 @@ function assemble!(R::BoundaryDataRestriction, sol, SC; kwargs...)
     nvals = length(fixeddofs)
 
     ## assign fixed dofs and vals to restriction
-    n = length(SC.b.entries)
+    n = length(SC.b[R.parameters[:unknown]])
     R.parameters[:matrix] = sparse(fixeddofs, 1:nvals, ones(nvals), n, nvals)
     R.parameters[:rhs] = fixedvals
     R.parameters[:multiplier] = zeros(nvals)

src/common_restrictions/linear_functional_restriction.jl

@@ -65,14 +65,14 @@ end
 
 function assemble!(R::LinearFunctionalRestriction{T, Tv}, sol, SC; kwargs...) where {T, Tv}
 
-    n = length(SC.b.entries)
+    n = length(SC.b[R.parameters[:unknown]])
 
-    b = deepcopy(SC.b)
+    b = copy(SC.b)
     fill!(b.entries, 0.0)
 
     assemble!(nothing, b, sol, R.linear_operator, SC; assemble_rhs = true, kwargs...)
 
-    R.parameters[:matrix] = sparse(reshape(b.entries, n, 1))
+    R.parameters[:matrix] = sparse(reshape(view(b[R.parameters[:unknown]]), n, 1))
     R.parameters[:rhs] = Tv[R.value]
     R.parameters[:multiplier] = zeros(1)
 

src/solvers.jl

@@ -43,6 +43,7 @@ function compute_nonlinear_residual!(residual, A, b, sol, unknowns, PD, SC, free
 
     # add Lagrange residuals
     for rs in PD.restrictions
+        @info size(rs.parameters[:matrix]), length(rs.parameters[:multiplier])
         view(residual[rs.parameters[:unknown]]) .-= rs.parameters[:matrix] * rs.parameters[:multiplier]
     end