Problem with parametrization of two line fields as a 4-field

[ioannaMitropoulou]

I am trying to integrate a 4-field as two overlayed line fields, but I get some parametrization results that do not make sense, like this:

(I've turned off the seamless integration just to save some time)

As a reference, I would expect this field's integration to look like this:

In the code (main.cpp) I have the following steps:

1. loading the two fields U and Urot,
2. calculate the matching of U,
3. combine them in a CCW 4-field with matching = U_matching * 2
4. and integrate.

I also added a curl-reduction step between steps 2 and 3, but I've commented it out as it makes no difference in this case.

If I calculate the matching of the 4-field (after step 3), then the result gives a correct integration, but it's not what I want because then the matching is not constrained to line-field. (You can test that by toggling the bool only_line_field_matching in main.cpp, line 310).

I'm probably missing something related to the matching, but I can't figure out what it is. If u have some time, could you take a quick look? Thank you very much!!

You can find the code here:
https://drive.google.com/file/d/1IOYkFzYZA8s4vBvat1JcvwjM4DZWXeg-/view?usp=sharing

[avaxman]

Hi! I am almost done with the new version of Directional where this is revamped into new code that might likely solve this issue, so I choose not to relate in the current version. It's going out in about two weeks, so I hope that's fine.

[ioannaMitropoulou]

Sounds good, thanks!