Re: [Getfem-users] Re : Generic elliptic Brick

[Yves Renard]

Dear Jean-François,

You would better use the new bricks because the old ones won't be improved 
anymore.

You are right, the correct syntax is 

md.add_generic_elliptic_brick(mim, 'varname', 'C')

where you first add the tensor C with

add_initialized_data('A',C)

A priori (I am not completely familiar with arrays in python) C could be a 
simple vector with the vlaues in the fortran order or an array with several 
dimensions.

yves.

On mercredi 9 février 2011, jeff witz wrote:
> I try to use this method anyway, but I don't succeed setting the 'A'
> parameter
> 
> I defined the good tensor for a simple composite case :
> 
> C=gf.MdBrick.get(b0, 'param','A')
> C[:]=0
> 
> # Definition of the isotropic transverse tensor in is orthotropic plan in
> plan stress
> El=120e6
> Et=8e6
> G=4
> Nult=0.3
> Nutl=0.01875
> C[0,0,0,0]=El/(1-Nult*Nutl)
> C[0,0,1,1]=Nult*Et/(1-Nult*Nutl)
> C[1,1,0,0]=C[0,0,1,1]
> C[0,1,0,1]=G
> 
> 
> But once there I don't find any method allowing me to set this tensor to
> 'A' og the brick.
> 
> I also try to use the new way for bricks as to know :
> md.add_generic_elliptic_brick(....)
> 
> but I don't find the good syntaxe.
> Could you give me the good syntaxe ?
> 
> 
> I think if I can use this last method I could set the model with :
> add_initialized_data('A',C)


-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------