Mercurial > repos > public > sbplib

diff +multiblock/DiffOpTimeDep.m @ 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 acb58769610e
children
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