Mercurial > repos > public > sbplib

changeset 883:76efb6a7b466 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add Jonatans bugfix to bcSetup.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 02 Nov 2018 10:44:59 -0700
parents 14fee299ada2
children 7d4f57725192
files +scheme/bcSetup.m
diffstat 1 files changed, 7 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/bcSetup.m	Fri Oct 26 16:35:23 2018 -0700
+++ b/+scheme/bcSetup.m	Fri Nov 02 10:44:59 2018 -0700
@@ -28,14 +28,14 @@
         [localClosure, penalty] = diffOp.boundary_condition(bcs{i}.boundary, bcs{i}.type);
         closure = closure + localClosure;
 
-        [ok, isSym, data] = parseData(bcs{i}, penalty, diffOp.grid);
+        [ok, isSymbolic, data] = parseData(bcs{i}, penalty, diffOp.grid);
 
         if ~ok
             % There was no data
             continue
         end
 
-        if isSym
+        if isSymbolic
             symbolicData{end+1} = data;
         else
             gridData{end+1} = data;
@@ -85,8 +85,10 @@
     end
 end
 
-function [ok, isSym, dataStruct] = parseData(bc, penalty, grid)
+function [ok, isSymbolic, dataStruct] = parseData(bc, penalty, grid)
     if ~isfield(bc,'data') || isempty(bc.data)
+        isSymbolic = [];
+        dataStruct = struct();
         ok = false;
         return
     end
@@ -96,14 +98,14 @@
 
     if nArg > 1
         % Symbolic data
-        isSym = true;
+        isSymbolic = true;
         coord = grid.getBoundary(bc.boundary);
         dataStruct.penalty = penalty;
         dataStruct.func = bc.data;
         dataStruct.coords = num2cell(coord, 1);
     else
         % Grid data
-        isSym = false;
+        isSymbolic = false;
         dataStruct.penalty = penalty;
         dataStruct.func = bcs{i}.data;
     end