Max number of border

General Discussion
1

Hi everyone,

I have a code with a surface composed of n quadrilaterals glued together, and because n > 256, the number of borders is greater than 1024, so when I try to run my code, i have the following error :

"Sorry number of parameters > 1024
current line = 5578
Compile error : Sorry number of parameters > 1024
line number :5578, ,
error Compile error : Sorry number of parameters > 1024
line number :5578, ,
code = 1 mpirank: 0 "

Does anyone have an idea how to maybe change the parameters limit, or to lower the number of borders? Heres my code,

Have a great day!

Émile

code_surface_tri_7g_#1_30n.edp (261.5 KB)

2

I think this not the got way to build the mesh

you can use multi border to build the mesh.

like in mesh.edo in tutorail directory exemples.

border G1i ( t=pi + theta1 ,pi + theta2 ;i) { x=( cex+R * cos( t ))+i ;
y= cey + R * sin( t ) ; } ;
border G2i ( t=pi - theta2 ,pi - theta1;i ) { x=( cex+R * cos( t ))+i ;
y= -cey + R * sin( t ) ; } ;
border G3i ( t= 1. , 0. ;i) { x= b*cos(alpha) * t +i ;
y= b*sin(alpha) * t ; } ;
border G4i ( t= 0. , 1. ;i) { x= b*cos(alpha) * t +i ;
y= - b*sin(alpha) * t ; } ;
int[int] NN(504); NN = n;
plot(G1i(NN)+G2i(NN)+G3i(NN)+G4i(NN));
3

Bonjour,

Merci beaucoup pour l’aide, je ne connaissais pas la méthode multi border. Ça marche bien pour faire les mesh, par contre j’ai un “fespace” après avec des conditions périodiques que je n’arrive pas à implémenter, vu que chaque quadrilatère a son propre collage. Voilà mon code original :

“fespace Vh(Th , P2 , periodic =[ [ 3, y ] , [ 8 , -y ] ,
[ 1, y ] , [ 10 , -y ] ,
[ 7, y ] , [ 4 , -y ] ,
[ 5, y ] , [ 34 , -y ] ,
[ 11, y ] , [ 24 , -y ] , etc.”

Est-ce que c’est possible et si oui sauriez vous comment faire référence à des bordures de quadrilatères en particulier dans mes conditions périodiques (quelque chose comme … [G1i(1), y] [G2i(3), -y] , [G3i(2), y] [G4i(1), -y] …)?

Merci encore,

Émile

4

J’ai changé les label en 1,2,3,4 pour les labels avec CL de périodicité
calcula avec votre tableau découpé en 2

maintenu la matrice n’est pas inversible, peut être un pb

ce truc marche si il n’y a pas overlapping sur les cl de périodicité.

Je n’est pas compris votre structure de CL donc je ne peux pas vous aide pls.

code_surface_tri_7g_#1_30n.edp (33.4 KB)

5

Qu’est-ce que vous voulez dire par CL?

6

Bonjour , il me semble que CL est l’abréviation de Conditions Limites

7

Salut,

Voilà la structure de mes conditions limites. En gros, chaque collage de bordures est écrit par paire. Par exemple, le segment [ 3, y ] , [ 8 , -y ] , [ 1, y ] , [ 10 , -y ] signifie que la bordure 3 est collée avec la bordure 8 et la bordure 1 avec la 10, les deux par une réflexion d’axe x (chaque point est collé avec celui ayant la valeur en y opposée). Ça revient à coller la 3e bordure (de type G3i) du premier quadrilatère avec la 4e du 2e quadrilatère, et la 1re du premier quadrilatère avec la 2e du 3e quadrilatère (10=2*4 + 2).

J’espère que c’est plus clair.

Émile