Hi py17js,
The error is caused by the mesh - there are vertices that are not in any elements and/or double vertices (two vertices in single spatial position).
The double vertices can be fixed by
python3 script/convert_mesh.py -m big_wire.mesh big_wire-new.mesh
but even after that there are some extra vertices. Nevertheless, the probe seems to work after some modifications:
def probe_hook(data, probe, label, problem): import matplotlib.pyplot as plt import matplotlib.font_manager as fm
def get_it(name, var_name):
var = problem.create_variables([var_name])[var_name]
var.set_data(data[name].data[var.field.vertex_remap_i])
pars, vals = probe(var)
vals = vals.squeeze()
return pars, vals
results = {}
results['Temperature'] = get_it('T', 'T')
print(results)
fig = plt.figure()
plt.clf()
fig.subplots_adjust(hspace=0.4)
plt.subplot(111)
pars, vals = results['Temperature']
plt.plot(pars, vals, label=r'$T$',
lw=1, ls='-', marker='+', ms=3)
plt.ylabel('T')
plt.xlabel('probe %s' % label, fontsize=8)
plt.legend(loc='best', prop=fm.FontProperties(size=10))
return plt.gcf(), resultsYou can check your mesh using
python3 script/show_mesh_info.py -d big_wire.mesh
- among other things it outputs the mesh Euler characteristic, which should be 1 for a simple block mesh. It is 532 for your mesh and 13 for the mesh with double vertices merged. So while the example seems to work, try fixing the mesh first.
r.
On 08. 03. 21 14:33, py17js@leeds.ac.uk wrote:
Hi Robert,
Thanks for the response to my last post. Your help has allowed my to obtain the results i was hoping for.
I have been able to probe some of the file however, i am having difficulting probing temperaure from an electric-thermal problem.
The error i am getting is : ValueError: incompatible data shape! (5693 (variable) == 6224 (data)), which is to do the the DOF of the variable but i can figure out how to solve it.
I have atttached my mesh, file, and probe file.
https://www.mediafire.com/file/n9qj8iox15z00jg/big_wire.mesh/file
https://www.mediafire.com/file/t5qsv2xxyj7oiuk/probe_wire.py/file
https://www.mediafire.com/file/0v35d0v9taloiph/wire.py/file
SfePy mailing list -- sfepy@python.org To unsubscribe send an email to sfepy-leave@python.org https://mail.python.org/mailman3/lists/sfepy.python.org/ Member address: cimrman3@ntc.zcu.cz