Hi,
I am a Ph.D. student at the Indian Institute of Science, Bangalore, India.
I have been working on a project that requires a polyhedron.
The polyhedron has constraints with floating-point variables, as shown below in the following code snippet.
Wher p,u,v and w are floating-point values between 0 and 1, u+v+w =1.
Variable xp0(0);
Variable yp0(1);
Variable zp0(2);
Variable u(0);
Variable v(1);
Variable w(2);
Variable p(0);
int x0 = vertices[insideVertex*3+0];
int y0 = vertices[insideVertex*3+1];
int z0 = vertices[insideVertex*3+2];
int x1 = vertices[outsideVertex*3+0];
int y1 = vertices[outsideVertex*3+1];
int z1 = vertices[outsideVertex*3+2];
pd->add_constraint( u+v+w == 1);
pd->add_constraint( u>=0) ;
pd->add_constraint( v>=0);
pd->add_constraint( w>=0);
pd->add_constraint( p>=0);
// pd->add_constraint( q>=0);
// pd->add_constraint( p+q ==1);
pd->add_constraint( ( (( p*(x0 - xp0)+ (1-p)*(x1 - xp0)) == (u*-35821+v*0+(w)*35821)) );
pd->add_constraint( ((p*(y0 - yp0)+ (1-p)*(y1 - yp0)) == (u*35821+v*0+(w)*35821)) );
pd->add_constraint( ((p*(z0 - zp0)+ (1-p)*(z1 - zp0)) == (u*-100000+v*0+(w)*-100000)));
pd->add_constraint( -( (24*p*(x0-xp0))+(24*(1-p)*(x1-xp0))) >= ((PixelX - 24)*(p*(z0-zp0)+ (1-p)*(z1-zp0)) ) );
pd->add_constraint( -( (24*p*(x0-xp0))+(24*(1-p)*(x1-xp0))) < ((PixelX+1 - 24)*(p*(z0-zp0)+ (1-p)*(z1-zp0)) ) );
pd->add_constraint( ( 24*p*(y0-yp0) + 24* (1-p)*(y1-yp0) ) >= ( (PixelY-24) * (p*(z0-zp0)+ (1-p)*(z1-zp0)) ) );
pd->add_constraint( ( 24*p*(y0-yp0) + 24* (1-p)*(y1-yp0) ) < ( (PixelY-24+1) * (p*(z0-zp0)+ (1-p)*(z1-zp0)) ) );
While compiling the program, I am getting the following error;
error: no match for ‘operator*’ (operand types are ‘Parma_Polyhedra_Library::Variable’ and ‘Parma_Polyhedra_Library::Linear_Expression’)
116 | pd->add_constraint( ( (( p*(x0 - xp0)+ (1-p)*(x1 - xp0)) == (u*-35821+v*0+(w)*35821)) );
| ~^~~~~~~~~~~
| | |
| | Parma_Polyhedra_Library::Linear_Expression
| Parma_Polyhedra_Library::Variable
Is there any way to get rid of this error?
Is it possible to write my constraints using PPL (using either floating-point expression or intervals)?
Thanks
HABEEB
