US 2024/0193333 A1 · Multilevel distributed parallel computing method for integrated circuit board simulation

Abstract

A multilevel distributed parallel computing method for integrated circuit board simulation, belonging to the technical field of integrated circuit simulation, includes following steps: S100, calculating number of nodes; S200, allocating tasks; S300, calculating and solving; and S40, judging convergence. By dynamically paralleling multiple computers, the present invention can maximize utilization of computing resources, improve simulation efficiency, and finally achieve a purpose of shortening design cycle. And in the simulation process, through the dynamic result analysis method, frequency points that need to be solved can be actually found, and accurate full-band results with the least number of solutions is completed.

Claims6

Dependent claims

1-5. (canceled)

A multilevel distributed parallel computing method for integrated circuit board simulation, wherein the method comprises following steps:

(S100), calculating a number of nodes;
(S200), allocating tasks;
(S300), calculating and solving;
and
(S400), judging convergence;
in the step (S300), calculating and solving comprises specifically solving distribution of electric fields and magnetic fields in the solution area at each of the computing nodes by a finite element method and extracting performance parameters of a circuit according to distribution of electromagnetic fields in a solution area for subsequent simulations;
wherein in the step (S300) the finite element method is used, specifically a finite element wave equation as shown here:
∫
Ω
W
i
·
[
∇
×
(
1
μ
r
⁢
∇
×
E
)
-
k
0
2
⁢
ε
r
⁢
E
]
⁢
d
⁢
Ω
=
-
jk
0
⁢
Z
0
⁢
∫
Ω
W
i
·
J
imp
⁢
d
⁢
Ω
;
wherein the distribution in the solution area in step S300 comprises grid division of the solution area;
the grid division of the solution area comprises tetrahedron division;
and the tetrahedron division is as follows:
by introducing a barycentric coordinate system and taking a tetrahedron unit year as an example, four vertices of each tetrahedron can be expressed as:
(L1e, L2e, L3e, L4e);
the tetrahedron comprises six edges with serial numbers 1, 2, 3, 4, 5 and 6, and the edges are placed as a vector basis function:
Nie=Wi1i2|ie=(Li1e∇Li1e−Li1e∇Li1e)|ie;
among which, lie represents a length of an ith edge;
a unit vector e1 is defined as a unit vector pointing from a node 1 to a node 2, and a gradient relative to gravity center coordinates can be described by following formulas e1·∇L1e=−1/l12 and e1·∇L2e=1/l1e;
among which, ;ie represents a length of an ith edge;
a unit vector ei is defined as a unit vector pointing from a node 1 to a node 2, and a gradient relative to gravity center coordinates is described by following formulas e1·∇L1e=−1/l12 and e1·∇L2e=1/l1e;
wherein an integral region in the equation is divided into several sub-regions Ωn, and a matrix equation is obtainable by using the finite element analysis method;
a left term of the matrix equation comprises two terms:
E
ij
e
=
∫
∫
∫
V
o
(
∇
×
N
i
e
)
·
(
∇
×
N
j
e
)
⁢
dV
⁢
F
ij
e
=
∫
∫
∫
V
e
N
i
e
·
N
j
e
⁢
dV
a right term of the matrix equation is as follows:
b=−jk0Z0fedgekNk*d1=−jk0Z0ξ
according to the above equations, a sparse matrix equation Ax=b is obtained;
the matrix equation is solved by multi frontal algorithm;
by running solution resources used in this step in a local computer or multi-computer parallel, which is dynamically adjusted according to fitting results;
wherein the multilevel distributed parallel computing method also comprises recovering the solution results and fitting frequency response with a formula
f
⁡
(
s
)
≈
∑
n
=
1
N
c
n
s
-
a
n
+
d
+
sh
,
where cn is a residue term and an is a pole;
a fitting process comprises multiplying a left end and a right end of the defined frequency response by an unknown function
σ
⁡
(
s
)
=
∑
n
=
1
N
s
-
a
n
_
+
1
simultaneously, and an overdetermined matrix equation about unknown quantities cn, d, h, and is obtained by substituting a few frequency points into the unknown function:
(
∑
n
=
1
N
c
n
s
-
a
n
_
)
=
(
∑
n
=
1
N
s
-
a
n
_
+
1
)
⁢
f
⁡
(
s
)
an unknown quantity obtained by solving the unknown function is used to describe a formula of the frequency response:
f
⁡
(
s
)
=
(
σ
⁢
f
)
fit
⁢
(
s
)
σ
fit
(
s
)
=
h
⁢
∑
n
=
1
N
+
1
(
s
-
z
n
)
∑
n
=
1
N
(
s
-
)
using a solution of σ (s)=0 as a new pole, a new overdetermined matrix equation can be generated to solve unknown quantities cn, d, and h.

245 more elements

Dependencies are inferred from the claim text.

Multilevel distributed parallel computing method for integrated circuit board simulation

Life of this patent

Multilevel distributed parallel computing method for integrated circuit board simulation

Life of this patent

Abstract

Claims6