Added method setGenericBCPBCNew(string). Now, the periodic points intersecting with other boundary type (phys. group) are considered non-PBC, i.e., the corner points. This modification reduces the periodic indices vectors and it's giving reasonable output for rising movingFrame. However, the pressure is not converging fast. The addition of the periodic conditions inside setGenericBC standard was not well received.
InOut: adding saveVTKPBC to pass beta gradient value as argument. Model: adding elseif condition for wallNoSlip in setGenericBC() and final else for nothing. This is to remove the exclusion condition of wallNoSlip in case of any B.C. setting.
Simulator: I've uncommented the line sumIndexPBCScalarNew() for the r.h.s. vector of pressure system. This solved the problem of corner points, but not the slow convergence of pressure.
Including coupledPBCNew(). However, the method is still wrong, requiring extraction of matrices for operations. Check.