Luis, The subject filter has an Array Bounds Read purify error. I think the logic for handling boundary faces is flawed. I put a cout in the boundary face loop and printed **inner_it. Some of its values are strange. Hope all is well, Bill