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

| main | Tutorials | Functions | website |