for(int kk=0;kk<3;kk++)
CPPUNIT_ASSERT_DOUBLES_EQUAL(expectP3[kk],res[kk],1e-12);
}
if(rank==4)
{
const double expectP4[6]={47.0, 47.0, 47.0, 39.0, 39.0, 31.0};
for(int kk=0;kk<3;kk++)
CPPUNIT_ASSERT_DOUBLES_EQUAL(expectP3[kk],res[kk],1e-12);
}
if(rank==4)
{
const double expectP4[6]={47.0, 47.0, 47.0, 39.0, 39.0, 31.0};