fixes to simpleMatrix ... really need to get a copy from work

src/OpenFOAM/matrices/simpleMatrix/simpleMatrix.C

@@ -31,11 +31,24 @@ License
 template<class Type>
 Foam::simpleMatrix<Type>::simpleMatrix(const label mSize)
 :
-    scalarSquareMatrix(mSize),
+    scalarSquareMatrix(mSize, mSize, pTraits<scalar>::zero),
     source_(mSize, pTraits<Type>::zero)
 {}
 
 
+template<class Type>
+Foam::simpleMatrix<Type>::simpleMatrix
+(
+    const label mSize,
+    const scalar t,
+    const Field<Type>& source
+)
+:
+    scalarSquareMatrix(mSize, mSize, t),
+    source_(mSize, source)
+{}
+
+
 template<class Type>
 Foam::simpleMatrix<Type>::simpleMatrix
 (

src/OpenFOAM/matrices/simpleMatrix/simpleMatrix.H

@@ -74,9 +74,12 @@ public:
 
     // Constructors
 
-        //- Construct given size
+        //- Construct given size, initialize coefficients to zero.
         simpleMatrix(const label);
 
+        //- Construct given size and initial coefficient value
+        simpleMatrix(const label, const scalar, const Field<Type>&);
+
         //- Construct from components
         simpleMatrix(const scalarSquareMatrix&, const Field<Type>&);