e=((p1+p2)+(p1-p2));

     e=((p1+p2)+(p1-p2));

     e.simplify();

     e.simplify();

     buf << "Coeff. of p2 should be 0";

     buf << "Coeff. of p2 should be 0";

     //    std::cout<<e[p1]<<e[p2]<<e[p3]<<std::endl;

     //    std::cout<<e[p1]<<e[p2]<<e[p3]<<std::endl;

     check(e.find(p2)==e.end(), buf.str());

     check(e.find(p2)==e.end(), buf.str());

     e=((p1+p2)+(p1-0.99*p2));

     e=((p1+p2)+(p1-0.99*p2));

     double tolerance=0.001;

     double tolerance=0.001;

     e.simplify(tolerance);

     e.simplify(tolerance);

     buf << "Coeff. of p2 should be 0.01";

     buf << "Coeff. of p2 should be 0.01";

     check(e[p2]>0, buf.str());

     check(e[p2]>0, buf.str());