Mercurial > repos > public > sbplib

changeset 711:397d1b22cc37 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed closure and penalty ao that sparse matrices are created in DiffOpTimeDep
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 30 Nov 2017 18:01:21 +0100
parents 023a95f63950
children c9147e05d228
files +multiblock/DiffOpTimeDep.m
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/+multiblock/DiffOpTimeDep.m	Tue Nov 21 18:12:00 2017 +0100
+++ b/+multiblock/DiffOpTimeDep.m	Thu Nov 30 18:01:21 2017 +0100
@@ -163,15 +163,15 @@
                     [closure, penalty] = obj.singleBoundaryCondition(boundary, type);
                 case 'multiblock.BoundaryGroup'
                     nBlocks = obj.grid.nBlocks();
-                    %[n,m] = size(obj.D);
+                    [n,m] = size(obj.D{1,1}(0));
                     %closure = sparse(n,m);
                     %penalty = sparse(n,0);
                     % closure =@(t)0;
                     % penalty = @(t)0;
                     for i = 1:nBlocks
                         for j = 1:nBlocks
-                            closure{j,i} = @(t)0;
-                            penalty{j,i} = @(t)0;
+                            closure{j,i} = @(t)0*sparse(n,m);
+                            penalty{j,i} = @(t)0*sparse(n,1);
                         end
                     end