Mercurial > repos > public > sbplib

changeset 865:1cc5a0d26453 bcSetupExperiment

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use existing functions
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 25 Jul 2018 15:43:26 -0700
parents c51b3f945c65
children dda1caa55eaf
files +multiblock/DiffOp.m
diffstat 1 files changed, 2 insertions(+), 9 deletions(-) [+]
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Wed Jul 25 15:24:37 2018 -0700
+++ b/+multiblock/DiffOp.m	Wed Jul 25 15:43:26 2018 -0700
@@ -197,15 +197,8 @@
             [blockClosure, blockPenalty] = obj.diffOps{I}.boundary_condition(name, type);
 
             % Expand to matrix for full domain.
-            div = obj.blockmatrixDiv;
-            temp = blockmatrix.zero(div);
-            temp{I,I} = blockClosure;
-            closure = blockmatrix.toMatrix(temp);
-
-            div{2} = size(blockPenalty, 2); % Penalty is a column vector
-            p = blockmatrix.zero(div);
-            p{I} = blockPenalty;
-            penalty = blockmatrix.toMatrix(p);
+            closure = multiblock.local2globalClosure(blockClosure, obj.blockmatrixDiv, I);
+            penalty = multiblock.local2globalPenalty(blockPenalty, obj.blockmatrixDiv, I);
         end
 
         function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)