SHOW_PARABOLIC Plots the solution on a domain consisting of triangles and parallelograms SHOW_PARABOLIC(ELEMENTS3, ELEMENTS4, COORDINATES, U) Plots the solution U on the triangulation descibed by ELEMENTS3, ELEMENTS4 and COORDINATES using a two dimensional spline interpolation. ELEMENTS3 denotes a set of triangles with dimension (no. of triangles) x 3. ELEMENTS4 denotes a set of parallelograms (dimension (no. of parallelograms) x 4. COORDINATES is a (no. of coordinates) x 2 matrix which contains the coordiantes of the nodes of ELEMENTS3 and ELEMENTS4. U is a matrix of size Nx(no. of coordinates) which contains the function value at the nodes i. N is the number of time steps solutions are stored for. This function should not be modified. The code is available at http://anmc.epfl.ch/ and described in further detail in A. Abdulle and A. Nonnenmacher "A short and versatile finite element multiscale code for homogenization problems" Computer Methods in Applied Mechanics and Engineering, http://dx.doi.org/10.1016/j.cma.2009.03.019 Please cite this article in any publication describing research performed using the software. Email : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch Last updated : 04/29/2009 with MATLAB 7.4 FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher. The software is provided free for non-commercial use unter the terms of the GNU General Public License. See "copyright.m" for full details.
0001 function show_parabolic(Elements3,Elements4,Coordinates,U,N) 0002 %SHOW_PARABOLIC Plots the solution on a domain consisting of triangles and parallelograms 0003 % SHOW_PARABOLIC(ELEMENTS3, ELEMENTS4, COORDINATES, U) 0004 % Plots the solution U on the triangulation descibed by ELEMENTS3, 0005 % ELEMENTS4 and COORDINATES using a two dimensional spline interpolation. 0006 % 0007 % ELEMENTS3 denotes a set of triangles with dimension (no. of triangles) x 3. 0008 % 0009 % ELEMENTS4 denotes a set of parallelograms (dimension (no. of parallelograms) x 4. 0010 % 0011 % COORDINATES is a (no. of coordinates) x 2 matrix which contains the 0012 % coordiantes of the nodes of ELEMENTS3 and ELEMENTS4. 0013 % 0014 % U is a matrix of size Nx(no. of coordinates) which contains the function 0015 % value at the nodes i. 0016 % 0017 % N is the number of time steps solutions are stored for. 0018 % 0019 % 0020 % This function should not be modified. 0021 % 0022 % 0023 % The code is available at http://anmc.epfl.ch/ and described in 0024 % further detail in 0025 % 0026 % A. Abdulle and A. Nonnenmacher 0027 % "A short and versatile finite element multiscale code for 0028 % homogenization problems" 0029 % Computer Methods in Applied Mechanics and Engineering, 0030 % http://dx.doi.org/10.1016/j.cma.2009.03.019 0031 % 0032 % Please cite this article in any publication describing research 0033 % performed using the software. 0034 % 0035 % 0036 % Email : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch 0037 % Last updated : 04/29/2009 with MATLAB 7.4 0038 % 0039 % FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher. 0040 % The software is provided free for non-commercial use unter the terms of 0041 % the GNU General Public License. See "copyright.m" for full details. 0042 0043 % 0044 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0045 0046 0047 % graphic representation 0048 u0max=max(U(:,1)); 0049 u0min=min(U(:,1)); 0050 subplot(2,2,1); 0051 show(Elements3,Elements4,Coordinates,full(U(:,1))); 0052 title('u(0)'); zlim([u0min,u0max]);caxis([u0min,u0max]); 0053 0054 subplot(2,2,2); 0055 show(Elements3,Elements4,Coordinates,full(U(:,4))); 0056 title('u(0.3)'); zlim([u0min,u0max]);caxis([u0min,u0max]); 0057 0058 subplot(2,2,3); 0059 show(Elements3,Elements4,Coordinates,full(U(:,8))); 0060 title('u(0.7)'); zlim([u0min,u0max]);caxis([u0min,u0max]); 0061 0062 subplot(2,2,4); 0063 show(Elements3,Elements4,Coordinates,full(U(:,11))); 0064 title('u(1)'); zlim([u0min,u0max]);caxis([u0min,u0max]); 0065 0066 end