Commit c3f3802a authored by Antoine MARTIN's avatar Antoine MARTIN

add comments

parent 31874bd1
......@@ -22,6 +22,7 @@ def rec_fill(t,N):
for ii in prange(N[0]):
rec_fill(t[ii],N)
#La fonction Ref est un algorithme de type Moulinec et Suquet qui permet de calculer la conductivit homognise (kref) partir d'une 'image' (ki, qui est un tableau contenant les conductivits de chaque pixel de l'image). N est la taille de l'image. tau_field est le champ de polarisation qui sert de point de dpart de l'algorithme. Et k0 est la conductivit de rfrence utilise pour le calcul.
def Ref(tau_field,N,k0,ki):
d=len(N)
filter_level=2
......
......@@ -36,6 +36,7 @@ def initialize_filters(N,K,n):
filters[i][ki][p] = (np.sin(pi*z)/(n*np.sin(pi*z/n)))**2
return filters
#La fonction initializeGreen prend en argument N, la taille de la grille/image (grille NxN) et renvoie les lments dont on a besoin pour utiliser la fonction operate_field
def initializeGreen(N,filter_level=2):
d = len(N)
numComp=d
......@@ -105,7 +106,7 @@ def operate_fourier_field(x,y,tupleK,N,frequencies,filter_level,filters,tupleFil
#do not use tuple as arguments in parallel mode -> remove tupleFIlters
#if not enough, remove tupleK and use scalar index to flattened x,y
#La fonction operate_field permet de calculer, pour un champ 'x' (qui reprsente un champ de polarisation), la dformation cre par la polarisation 'x'. C'est--dire, -Gamma(x). N donne la taille de la grille/image (la taille est NxN). L'argument k0 est la conductivit du milieu de rfrence utilis pour le calcul. Les autres arguments sont ceux qui sont issus de la fonction initializeGreen
def operate_field(x,yFourier,fft,ifft,tupleK,N,frequencies,filter_level,filters,tupleFilters,k0):
#start = time.time()
xFourier=fft(x)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment