integr_order
| main | Tutorials | Functions | website |
Returns the gaussian integration points and weights for a given type of elements.
NOTE: Need to go through this function again
Version : 1.0
Author : George Kourakos
email: giorgk@gmail.com
web : http://groundwater.ucdavis.edu/msim
Date 18-Mar-2014
Department of Land Air and Water
University of California Davis
Contents
Usage
pint = integr_order(opt)
Input
opt: a structure with the following fields:
dim : dimension of elements 1, 2, 3
el_type: element type e.g. 'triangle','quad','prism','hex'
int_ord: THis defines the number of gaussing quadrature. (The input values are quite weird now and need to be fixed)
Output:
pint : [np x ndim+1] Martix that contains the integration points in local coordinates np is the number of gaussian quadrature points ndim is the dimension of the problem plus the integration weights
Example
Set the options for 2D triangles
opt.dim = 2;
opt.el_type = 'triangle';
For linear triange elements set
opt.int_ord = 3; pint = integr_order(opt)
pint =
0.5000 0 0.3333
0.5000 0.5000 0.3333
0 0.5000 0.3333
For quadratic triangle elements set
opt.int_ord = 6; pint = integr_order(opt)
pint =
0.3333 0.3333 0.2250
0.4701 0.0597 0.1324
0.4701 0.4701 0.1324
0.0597 0.4701 0.1324
0.1013 0.7974 0.1259
0.1013 0.1013 0.1259
0.7974 0.1013 0.1259