show_parabolic

PURPOSE ^

SHOW_PARABOLIC Plots the solution on a domain consisting of triangles and parallelograms

SYNOPSIS ^

function show_parabolic(Elements3,Elements4,Coordinates,U,N)

DESCRIPTION ^

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.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

Generated on Tue 21-Jul-2009 10:55:32 by m2html © 2003