tensor_a

PURPOSE ^

TENSOR_A Evaluate conductivity tensor A

SYNOPSIS ^

function a =tensor_a(x, y)

DESCRIPTION ^

TENSOR_A   Evaluate conductivity tensor A
   A = TENSOR_A(X,Y) returns function values at a discrete point given by
   X on the macro and Y on the micro domain (see article for details).
   This input data has to be chosen by the user. 
   X has dimension 1 x 2 and Y has dimension 1 x 2. A is a 2 x 2 matrix.

   See also HMMFEM2D, F, and G.


   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 a =tensor_a(x, y)
0002 %TENSOR_A   Evaluate conductivity tensor A
0003 %   A = TENSOR_A(X,Y) returns function values at a discrete point given by
0004 %   X on the macro and Y on the micro domain (see article for details).
0005 %   This input data has to be chosen by the user.
0006 %   X has dimension 1 x 2 and Y has dimension 1 x 2. A is a 2 x 2 matrix.
0007 %
0008 %   See also HMMFEM2D, F, and G.
0009 %
0010 %
0011 %   The code is available at http://anmc.epfl.ch/ and described in
0012 %   further detail in
0013 %
0014 %   A. Abdulle and A. Nonnenmacher
0015 %   "A short and versatile finite element multiscale code for
0016 %   homogenization problems"
0017 %   Computer Methods in Applied Mechanics and Engineering,
0018 %   http://dx.doi.org/10.1016/j.cma.2009.03.019
0019 %
0020 %   Please cite this article in any publication describing research
0021 %   performed using the software.
0022 %
0023 %
0024 %   Email           : assyr.abdulle@epfl.ch and achim.nonnenmacher@epfl.ch
0025 %   Last updated    : 04/29/2009 with MATLAB 7.4
0026 %
0027 %   FE_HMM2D is Copyright (C) 2009 A. Abdulle and A. Nonnenmacher.
0028 %   The software is provided free for non-commercial use unter the terms of
0029 %   the GNU General Public License. See "copyright.m" for full details.
0030 
0031 %
0032 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0033 
0034 
0035 aa= (1.5+sin(2*pi*y(1)))/(1.5+sin(2*pi*y(2))) + ...
0036     (1.5+sin(2*pi*y(2)))/(1.5+cos(2*pi*y(1))) + ...
0037     sin(4*x(1)*x(2)) + 1;
0038 
0039 a= eye(2).*aa;
0040 
0041 
0042 end

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