Bugfix fluid and solid solvers with different read and write meshes (#311)

Author: David Blom

src/fsi/fluidSolvers/foamFluidSolver.C

@@ -73,10 +73,10 @@ foamFluidSolver::foamFluidSolver(
         movingPatchesDisplOld[movingPatchIDs[patchI]] = vectorField( size, Foam::vector::zero );
     }
 
-    matrix readPositions;
-    getReadPositions( readPositions );
+    matrix writePositions;
+    getWritePositions( writePositions );
 
-    N = readPositions.rows();
+    N = writePositions.rows();
     dim = mesh.nGeometricD();
     data.resize( N, dim );
     data.setZero();

src/fsi/solidSolvers/foamSolidSolver.C

@@ -71,6 +71,11 @@ foamSolidSolver::foamSolidSolver (
         movingPatchIDs[patchI] = patchIndex;
     }
 
+    matrix writePositions;
+    getWritePositions( writePositions );
+
+    N = writePositions.rows();
+
     dim = mesh.nGeometricD();
     data.resize( N, dim );
     data.setZero();