pcylinder b1 1 2
dset r sqrt(2)/2
box b2 -1 -1 0 1+r 1+r 1

bop b1 b2
bopcut result

set square 18.5364
