Commit 0eb60316 authored by Jeremy BLEYER's avatar Jeremy BLEYER

Update greenOperatorConductionNumba.py

parent 64871f9b
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import time
import numpy as np
from math import pi
......@@ -36,7 +38,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 (c'est un 'tuple' de taille d o d est la dimension. Ex : en dimension 2, pour une image 128x128, N=(128,128)) et renvoie les lments dont on a besoin pour utiliser la fonction operate_field
#La fonction initializeGreen prend en argument N, la taille de la grille/image (c'est un 'tuple' de taille d o d est la dimension. Ex : en dimension 2, pour une image 128x128, N=(128,128)) 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
......@@ -106,7 +108,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 est la taille de la grille/image. 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
#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 est la taille de la grille/image. 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