# Changeset 521 for branches/Couenne/Couenne/src/util

Ignore:
Timestamp:
Apr 28, 2007 2:10:25 AM (13 years ago)
Message:

introducing better bounds for sin/cos. Tool for drawing cuts

File:
1 edited

Unmodified
Removed
• ## branches/Couenne/Couenne/src/util/drawCuts.cpp

 r515 minY =  COUENNE_INFINITY; if (1 || (img -> code () == COU_EXPRSIN) || if ((img -> code () == COU_EXPRSIN) || (img -> code () == COU_EXPRPOW) || (img -> code () == COU_EXPREXP) || (img -> code () == COU_EXPRCOS)) { printf (" ==> "); w -> print (std::cout); printf ("\n"); //    fprintf (stderr, " ==> "); w -> print (std::cerr); fprintf (stderr, "\n"); expression *lbe, *ube; int xi = indep -> Index (); printf ("looking into w_%d = f (x_%d)\n", w -> Index (), xi); //    fprintf (stderr, "looking into w_%d = f (x_%d)\n", //       w -> Index (), xi); indep -> getBounds (lbe, ube); if (y < minY) minY = y; printf ("#=# %.12e %.12e\n", x, y); fprintf (stderr, "%.12e %.12e\n", x, y); } maxY += (maxY-minY) / 20; minY -= (maxY-minY) / 20; } //lb -= 1; //ub += 1; lb -= (ub-lb) / 20; ub += (ub-lb) / 20; // plot lines defining constraint (only for cuts involving at } printf ("#=# #m=2,S=%d\n", (cs.rowCutPtr (jj) -> sense () == 'L') ? 10:11); fprintf (stderr, "#m=2,S=%d\n", (cs.rowCutPtr (jj) -> sense () == 'L') ? 10:11); printf ("#=# %.12e %.12e\n", lb0, (rhs - el [1] * lb0) / el [0]); printf ("#=# %.12e %.12e\n", ub0, (rhs - el [1] * ub0) / el [0]); fprintf (stderr, "%.12e %.12e\n", lb0, (rhs - el [1] * lb0) / el [0]); fprintf (stderr, "%.12e %.12e\n", ub0, (rhs - el [1] * ub0) / el [0]); } cg -> Problem () -> X () [xi] = curx; exit(0); } }
Note: See TracChangeset for help on using the changeset viewer.