Wrong intersection of polygon and linestring

[vissarion]

The following code

#include <iostream>
#include <vector>
#include <boost/geometry.hpp>

namespace bg = boost::geometry;

int main() {
    // Define the polygon
    bg::model::polygon<bg::model::d2::point_xy<double>> polygon;
    bg::read_wkt("POLYGON((1 1.5,3.5 4.5,4 6,2.5 4.5,1 1.5))", polygon);

    // Define the linestring
    bg::model::linestring<bg::model::d2::point_xy<double>> linestring;
    bg::read_wkt("LINESTRING(1.5 1.5,0.5 1.5)", linestring);

    // Compute the intersection
    std::tuple<
        bg::model::multi_point<bg::model::d2::point_xy<double>>,
        bg::model::multi_linestring<bg::model::linestring<bg::model::d2::point_xy<double>>>
    > output;
    bg::intersection(polygon, linestring, output);

    // Print the results
    auto linestrings = std::get<1>(output);
    auto points = std::get<0>(output);
    std::cout << "Intersection linestrings: " << bg::wkt(linestrings) << std::endl;
    std::cout << "Intersection points: " << bg::wkt(points) << std::endl;

    return 0;
}

Gives the following incorrect output:

Intersection linestrings: MULTILINESTRING((1 1.5,0.5 1.5))
Intersection points: MULTIPOINT()

on current develop with clang-15 on ubuntu 22.

[vissarion]

Issue update: the polygon is invalid (incorrect orientation).
If the corrected polygon is used i.e. POLYGON((1 1.5,3.5 4.5,4 6,2.5 4.5,1 1.5)) then the result is empty (still wrong thought).

Intersection linestrings: MULTILINESTRING()
Intersection points: MULTIPOINT()

[vissarion]

Similar issue, intersection of POLYGON((1 1,3 5,4 6,3 4,1 1)) with LINESTRING(0 0,1 1). Now the geometries intersect in a common vertex.