On 10/18/2012 02:04 PM, Alec Kalinin wrote:

Hello SfePy!

I am solving the Poisson problem and during my solution I got the warning:

Warning: /usr/local/lib/python2.7/dist-packages/sfepy/fem/fields_base.py:1258: RuntimeWarning: invalid value encountered in divide data_vertex /= nod_vol[:,nm.newaxis]

The code for this warning: def post_process(out, pb, state, extend = False): # evaluate gradient in nodes grad_data = pb.evaluate('ev_grad.i1.Omega(u)', mode = 'qp')

     grad_field = H1NodalVolumeField('grad', np.float64, (3, ),
                                     pb.domain.regions['Omega'])
     grad_var = FieldVariable('grad', 'parameter', grad_field, 3,
                              primary_var_name='(set-to-None)')
     *grad_var.data_from_qp(grad_data, pb.integrals['i1']) *

This warning is appeared only for the meshes with superfluous vertices. The code and the meshes is in the attachment. I think this warning my be related to the issue with superfluous vertices. See early discussion [1] and git commit [2]

Good catch, thanks for reporting this! It is now fixed in the master repository

• both your meshes work without causing any warnings.

Cheers, r.

[1] https://groups.google.com/forum/?fromgroups=#!topic/sfepy-devel/3Qae2M5bBps [2] https://github.com/sfepy/sfepy/commit/4ef0f5dbc7fac15942ba2031debf1296a35eb4...

Sincerely, Alec