diff --git a/examples/domains/customXY/domainsXY.m b/examples/domains/customXY/domainsXY.m
index 735e23a21..4a54f0493 100644
--- a/examples/domains/customXY/domainsXY.m
+++ b/examples/domains/customXY/domainsXY.m
@@ -27,7 +27,7 @@
 vfWater = 1 - totalVF;
 
 % Now convert the Volume Fractions to SLD's...
-siSLD = vfSilicon * bulk_in;
+siSLD = vfSilicon * bulk_in(contrast);
 oxSLD = vfOxide * 3.41e-6;
 
 % Layer SLD depends on whether we are calculating the domain or not...
diff --git a/targetFunctions/+domainsTF/processCustomFunction.m b/targetFunctions/+domainsTF/processCustomFunction.m
index f4ac574f3..b2d23bb84 100644
--- a/targetFunctions/+domainsTF/processCustomFunction.m
+++ b/targetFunctions/+domainsTF/processCustomFunction.m
@@ -18,11 +18,11 @@
         coder.varsize('output2',[10000 numberOfOutputColumns], [1 1]);
 
         if isnan(str2double(functionHandle))
-            [output1, subRoughs(i)] = callMatlabFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 1);
-            [output2, ~] = callMatlabFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 2);
+            [output1, subRoughs(i)] = callMatlabFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 1);
+            [output2, ~] = callMatlabFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 2);
         else
-            [output1, subRoughs(i)] = callCppFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 1);
-            [output2, ~] = callCppFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 2);
+            [output1, subRoughs(i)] = callCppFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 1);
+            [output2, ~] = callCppFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 2);
         end
 
         % If SLD is real, add dummy imaginary column
diff --git a/targetFunctions/+normalTF/processCustomFunction.m b/targetFunctions/+normalTF/processCustomFunction.m
index 2cc046fe6..0ae83aa6f 100644
--- a/targetFunctions/+normalTF/processCustomFunction.m
+++ b/targetFunctions/+normalTF/processCustomFunction.m
@@ -14,9 +14,9 @@
         coder.varsize('output',[10000 numberOfOutputColumns],[1 1]);
 
         if isnan(str2double(functionHandle))
-            [output, subRoughs(i)] = callMatlabFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 0);
+            [output, subRoughs(i)] = callMatlabFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 0);
         else
-            [output, subRoughs(i)] = callCppFunction(functionHandle, paramValues, bulkIns(i), bulkOuts, i, 0);
+            [output, subRoughs(i)] = callCppFunction(functionHandle, paramValues, bulkIns, bulkOuts, i, 0);
         end
 
         % If SLD is real, add dummy imaginary column