Implementing assemblePBCTwoPhaseNew() to handle Mrho, fint, etc for pbc. I should keep in mind that the if-block of LR copy should be implemented in this method, as well as in sumIndex family. However, rethink the LR copy for all the methods that would use it.
Rewriting ILAndIR. i) recasted to lists; ii) separated because in 3D theses lists will have variable dimensions since several points of interface might be periodic; iii) pointers changed. iv) CheckPBCOverlapping doesn't need removals in Master/SlaveIndices since iL/iR are also periodic; v) loops of elimination in assemblePBC don't need a special line for iL/iR.
setNormalAndKappaNew to remove 1D treatment.
Inserting elseif-condition of 'wallLeft/Right' in setGenericBC() to save pbc nodes. Remember that the 'wallNoSlip' condition is ascribed as last condition in the method. Hence, this elseif block prevents the pbc boundary of other Dirichlet conditions.
Removing methods from Periodic.
Updating...