Hi jeff, three questions:
* are using the svn version? * if you had the correct values, how they would be plotted? I'm asking you because if I know what you want, it might be a good idea to integrate this into the interface without having to write the code yourself, regards, Luis 2011/3/3 jeff witz <[email protected]> > Hello, > > I almost succeed in using matplotlib for plotting 2 dimensionnal linear > triangle ! > > So I have a two dimensional NonLinearelasticity problem, so I do the > computation and extract the solution > > u = md.variable('u') > > In order to be able to plot solution in matplotlib 1.0.1 I need node table > and element table. So I build this table : > > > NbNod=m.pts().shape[1] > NbElem = m.cvid().shape[0] > NbNodElem=m.pid_from_cvid(0)[0].shape[0] > ElemTab = np.zeros((NbElem,NbNodElem)) > for i in range(NbElem): > for j in range(NbNodElem): > ElemTab[i,j]=m.pid_in_cvids(i)[j] > > NodTab=np.transpose(m.pts()) > > The main issue remain in the fact that I don't understand how the u vector > is sorted. > So when I plot the results the values don't correspond to the computed > solution. Could you explain me how I can sort it. > > Logically this vector is 324 : 2* NbNod > > I try to use mf.basic_dof_nodes(), but I don't really understand the way > output vector > > mf.basic_dof_nodes().shape > Out[100]: (2, 324) > > which is equal to (2,2*NbNod) > > Could you help me to solve this out ? > > I think it would be great to be able to treat the data without leaving > ipython. > > I just want to plot a simple tripcolor (pcolor in matlab) with a quiver > over it. It will be quite simple once I have the good u vector > > I hope I can solve this before the Course I have to give Monday to show the > student that everything can be done simply in python. > > Best regards. > > > > > -- > Jean-François WITZ > > _______________________________________________ > Getfem-users mailing list > [email protected] > https://mail.gna.org/listinfo/getfem-users > >
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
