Mercurial > repos > public > sbplib

changeset 862:517d6019d63c bcSetupExperiment

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
I remove cell support in multiblock.DiffOp.boundary_condition
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 25 Jul 2018 15:17:16 -0700
parents 8e4b9a5d7c50
children 77aa40da2008
files +multiblock/DiffOp.m
diffstat 1 files changed, 7 insertions(+), 25 deletions(-) [+]
line wrap: on
line diff
--- a/+multiblock/DiffOp.m	Wed Jul 25 09:34:51 2018 -0700
+++ b/+multiblock/DiffOp.m	Wed Jul 25 15:17:16 2018 -0700
@@ -198,32 +198,14 @@
 
             % Expand to matrix for full domain.
             div = obj.blockmatrixDiv;
-            if ~iscell(blockClosure)
-                temp = blockmatrix.zero(div);
-                temp{I,I} = blockClosure;
-                closure = blockmatrix.toMatrix(temp);
-            else
-                for i = 1:length(blockClosure)
-                    temp = blockmatrix.zero(div);
-                    temp{I,I} = blockClosure{i};
-                    closure{i} = blockmatrix.toMatrix(temp);
-                end
-            end
+            temp = blockmatrix.zero(div);
+            temp{I,I} = blockClosure;
+            closure = blockmatrix.toMatrix(temp);
 
-            if ~iscell(blockPenalty)
-                div{2} = size(blockPenalty, 2); % Penalty is a column vector
-                p = blockmatrix.zero(div);
-                p{I} = blockPenalty;
-                penalty = blockmatrix.toMatrix(p);
-            else
-                % TODO: used by beam equation, should be eliminated. SHould only set one BC per call
-                for i = 1:length(blockPenalty)
-                    div{2} = size(blockPenalty{i}, 2); % Penalty is a column vector
-                    p = blockmatrix.zero(div);
-                    p{I} = blockPenalty{i};
-                    penalty{i} = blockmatrix.toMatrix(p);
-                end
-            end
+            div{2} = size(blockPenalty, 2); % Penalty is a column vector
+            p = blockmatrix.zero(div);
+            p{I} = blockPenalty;
+            penalty = blockmatrix.toMatrix(p);
         end
 
         function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)