Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 1 addition & 3 deletions Modelica/Electrical/Polyphase/Basic/MutualInductor.mo
Original file line number Diff line number Diff line change
Expand Up @@ -4,9 +4,7 @@ model MutualInductor "Linear mutual inductor"
parameter Real epsilon=1e-9 "Relative accuracy tolerance of matrix symmetry";
parameter SI.Inductance L[m, m] "Mutual inductance matrix";
initial equation
if abs(Modelica.Math.Matrices.det(L)) < epsilon then
Modelica.Utilities.Streams.print("Warning: mutual inductance matrix singular!");
end if;
assert(abs(Modelica.Math.Matrices.det(L)) >= epsilon, "Mutual inductance matrix singular!", AssertionLevel.warning);
equation
assert(sum(abs(L - transpose(L))) < epsilon*sum(abs(L)),"Mutual inductance matrix is not symmetric");
for j in 1:m loop
Expand Down
6 changes: 2 additions & 4 deletions Modelica/Electrical/QuasiStatic/Polyphase/Basic/MultiDelta.mo
Original file line number Diff line number Diff line change
Expand Up @@ -16,12 +16,10 @@ model MultiDelta
QuasiStatic.Polyphase.Basic.PlugToPins_n plugToPins_n(final m=m)
annotation (Placement(transformation(extent={{80,-10},{60,10}})));
protected
parameter Integer kP=if (mBasic<=2 or kPolygon<1 or kPolygon>integer(mBasic - 1)/2) then 1 else kPolygon;
final parameter Integer kP=if (mBasic<=2 or kPolygon<1 or kPolygon>integer(mBasic - 1)/2) then 1 else kPolygon;
equation
when initial() then
if (mBasic<=2 or kPolygon<1 or kPolygon>integer(mBasic - 1)/2) then
print("MultiDelta: replaced erroneous kPolygon = "+String(kPolygon)+" by kPolygon = 1");
end if;
assert(kP==kPolygon, "MultiDelta: replaced erroneous kPolygon = "+String(kPolygon)+" by kPolygon = 1", AssertionLevel.warning);
end when;
for k in 1:mSystems loop
for j in 1:(mBasic - kP) loop
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -5,9 +5,7 @@ model MutualInductor "Linear mutual inductor"
parameter Real epsilon=1e-9 "Relative accuracy tolerance of matrix symmetry";
parameter SI.Inductance L[m, m] "Mutual inductance matrix";
initial equation
if abs(Modelica.Math.Matrices.det(L)) < epsilon then
Modelica.Utilities.Streams.print("Warning: mutual inductance matrix singular!");
end if;
assert(Modelica.Math.Matrices.det(L) >= epsilon, "Mutual inductance matrix singular!", AssertionLevel.warning);
equation
assert(sum(abs(L - transpose(L))) < epsilon*sum(abs(L)),"Mutual inductance matrix is not symmetric");
for l in 1:m loop
Expand Down
Loading