parallel all with tbb, also find bug with mirror boundary type when using integral