SfePy

Download

sfepy@python.org

March 2013

6 participants
12 discussions

New Gallery Page
by Ankit Mahato March 20, 2013

March 20, 2013

Hi all, I have developed a new gallery page for SfePy which is more fluid. All suggestions are welcome before I commit it to my github repository. http://home.iitk.ac.in/~amahato/sfepy-gallery/ Note that on navigating to any particular problem the background image as well as the "CODE" link also changes redirecting to that particular example page. Regards, Ankit

2 1

Using the dw_diffusion_coupling term with inhomogenous materials
by Ben Derrett March 13, 2013

March 13, 2013

Hello, I've been using SfePy for about a week - it's a great piece of work. I'm working on a 2D problem involving a dw_diffusion_coupling term<http://sfepy.org/doc-devel/src/sfepy/terms/termsLaplace.html#sfepy.terms.te…>. For a homogenous material everything works as expected using the material definition 'mat_diffusion_coupling': ({'f': np.array([[0], [0]])},), > For an inhomogenous material, since the shape of this local parameter is (2, 1), I defined a function which takes the array of coordinates - coors - and returns an array of shape (coors.shape[0], 2, 1) like so def diffusion_coupling_coefficient(ts, coors, mode=None, **kwargs): > return {'f': np.zeros((coors.shape[0], 2, 1))} > This gives the error File ".../sfepy/terms/termsLaplace.py", line 199, in dw_fun > status = terms.mulATB_integrate(out, vg.bfg, mat * val, vg) > ValueError: operands could not be broadcast together with shapes > (5151,2,1) (5000,8,1,1) > , so I expect I'm doing something silly! How do I get this to work? (I've attached a minimal version of the code giving the error.) Many thanks, Ben

3 6