Fix handling of quantities in input sequences

[alihamdan]

Fixes #145

This required a few setters to abandon the ureg.wraps decorator but the alternative is very simple so I don't see a problem. I also tested other alternative like Q_.from_sequence but these don't always work as expected.

[alihamdan]

@benoit9126 did you have some time to take a look at this PR? I think we should cut a release after it is merged to fix several warnings from new pandas versions.

[benoit9126]

Yes, and both @Saelyos and I don't manage to reproduce this bug on our computer with Python 3.11 or Python 3.12. Could you please provide us the version of numpy and pint that you are using? The problem may be solved by modifying the dependencies requirements.

[alihamdan]

1. Checkout the develop branch
2. Create a python file bug.py with:

import numpy as np
from roseau.load_flow import *

# Create two buses
source_bus = Bus(id="sb", phases="abcn")
load_bus = Bus(id="lb", phases="abcn")

# Define the reference of potentials to be the neutral of the source bus
ground = Ground(id="gnd")
# Fix the potential of the ground at 0 V
pref = PotentialRef(id="pref", element=ground)
ground.connect(source_bus, phase="n")

# Create a LV source at the first bus
# Volts (phase-to-neutral because the source is connected to the neutral)
un = Q_(0.4 / np.sqrt(3), "kV")
source_voltages = [un, un * np.exp(-2j * np.pi / 3), un * np.exp(2j * np.pi / 3)]
vs = VoltageSource(id="vs", bus=source_bus, voltages=source_voltages)

# Add a load at the second bus
load = PowerLoad(id="load", bus=load_bus, powers=[10e3 + 0j, 10e3, 10e3])  # VA

# Add a LV line between the source bus and the load bus
# R = 0.1 Ohm/km, X = 0
lp = LineParameters("lp", z_line=(0.1 + 0.0j) * np.eye(4, dtype=complex))
line = Line(id="line", bus1=source_bus, bus2=load_bus, parameters=lp, length=2.0)

3. run python bug.py

Versions:

$ python --version
Python 3.11.5
$ poetry show | grep -E '(numpy|pint)'
numpy                         1.26.1       Fundamental package for array co...
pint                          0.22         Physical quantities module

This also reproduced with Python 3.12 and the latest version of pint and numpy 3 weeks ago

[benoit9126]

@alihamdan: @Saelyos is on it. He will rewrite a part of the wrapper to handle the case "array of quantities". Two advantages:

1. the feature will be available everywhere.
2. the improvements in the wrapper (Improve wraps performances hgrecco/pint#1866 and Wraps benchmark hgrecco/pint#1862) will be immediately available in RLF.