generalize error scripts to take an arbitrary matrix permeability

scripts/compute_1d_error.m

-function [totalRelativeFractureError] = compute_1d_error( coarseFile,referenceFile )
+function [totalRelativeFractureError] = compute_1d_error( coarseFile, referenceFile, kmat )
 %compute1dError Evaluates the norm in each of the fractures
 % 
 % The coarse file must contain a nFractures x 1 cell array with fracXAndP for
@@ -11,7 +11,10 @@ load(referenceFile,'p','t','x','k')
 load(coarseFile,'fracXAndPCell')
 maxp = max(p);
 minp = min(p);
-matrixIndices = find(k == 1);
+if nargin == 2
+    kmat = 1;
+end
+matrixIndices = find(k == kmat);
 t(matrixIndices) = [];
 p(matrixIndices) = [];
 

scripts/compute_2d_error.m

-function [relativeMatrixError] = compute_2d_error( coarseFile,referenceFile )
+function [relativeMatrixError] = compute_2d_error( coarseFile, referenceFile, kmat )
 %compute_2d_error Evaluates the norm in the entire matrix
 %   Detailed explanation goes here
 load(coarseFile,'T','P','Points')
@@ -12,7 +12,10 @@ load(referenceFile,'t','p','x','k','domainarea')
 
 points          = x;
 deltaP = max(p)-min(p);
-fracIndices = find(k ~= 1);
+if nargin == 2
+    kmat = 1;
+end
+fracIndices = find(k ~= kmat);
 t(fracIndices) = [];
 p(fracIndices) = [];
 nCells          = length(t);