Thanks a lot, it works !

Le lundi 4 juillet 2016 22:10:05 UTC+2, Robert Cimrman a écrit :

Hi Nicolas,

you need only to return your values in 'qp' mode (and make them floats). It
runs ok with those:

def get_heat_Coeff(ts, coors, mode=None, **kwargs):
     if mode != 'qp': return

     file=open('heat_Coeff.dat', 'r')
     tab=[]
     tab=file.readlines()
     file.close()
     length=len(tab)
     for i in range(length):
         tab[i] = float(tab[i].strip())
     val=nm.asarray(tab)
     val.shape = (coors.shape[0], 1,  1)
     return {'h' : val}

def get_exchange(ts, coors, mode=None, **kwargs):
     if mode != 'qp': return

     file=open('exchange.dat', 'r')
     tab=[]
     tab=file.readlines()
     file.close()
     length=len(tab)
     for i in range(length):
         tab[i] = float(tab[i].strip())
     val=nm.asarray(tab)
     val.shape = (coors.shape[0], 1, 1)
     return {'g' : val}

BTW. make sure that the order of the values corresponds to the ordering of the
quadrature point coordinates in coors argument.

r.

On 07/04/2016 03:27 PM, Nicolas DELAN wrote:
> Hi,
>
> I have some problem with my case:
>
> I have followed your instructions and i was inspired by the get_pars
> function of the  following example case :
> http://sfepy.org/doc-devel/examples/diffusion/poisson_functions.html
>
> But with me it don't works, the coors type returned is 'none' so the
> coors.shape returned nothing. Stangely i have understand that the .py
> configuration file is readed several time by sfepy, the first time coors is
> recognised and i can print the coordinates. The second time, coors become a
> none type variable and i don't understand why.
>
> Also could you take a look at my case ?
>
> I have attached my files to this post.
>
>