/*--------------------------------------------------------------------------*/
/* ALBERTA: an Adaptive multi Level finite element toolbox using */
/* Bisectioning refinement and Error control by Residual */
/* Techniques for scientific Applications */
/* */
/* file: lagrange_4_3d.c */
/* */
/* description: piecewise quartic Lagrange elements in 3d */
/* */
/*--------------------------------------------------------------------------*/
/* */
/* authors: Alfred Schmidt */
/* Zentrum fuer Technomathematik */
/* Fachbereich 3 Mathematik/Informatik */
/* Universitaet Bremen */
/* Bibliothekstr. 2 */
/* D-28359 Bremen, Germany */
/* */
/* Kunibert G. Siebert */
/* Institut fuer Mathematik */
/* Universitaet Augsburg */
/* Universitaetsstr. 14 */
/* D-86159 Augsburg, Germany */
/* */
/* http://www.mathematik.uni-freiburg.de/IAM/ALBERTA */
/* */
/* (c) by A. Schmidt and K.G. Siebert (1996-2003) */
/* */
/*--------------------------------------------------------------------------*/
#define N_BAS4_3D 35
static const REAL bary4_3d[N_BAS4_3D][N_LAMBDA] = {
{ 1.0, 0.0, 0.0, 0.0},
{ 0.0, 1.0, 0.0, 0.0},
{ 0.0, 0.0, 1.0, 0.0},
{ 0.0, 0.0, 0.0, 1.0},
{3.0/4.0, 1.0/4.0, 0.0, 0.0},
{ 0.5, 0.5, 0.0, 0.0},
{1.0/4.0, 3.0/4.0, 0.0, 0.0},
{3.0/4.0, 0.0, 1.0/4.0, 0.0},
{ 0.5, 0.0, 0.5, 0.0},
{1.0/4.0, 0.0, 3.0/4.0, 0.0},
{3.0/4.0, 0.0, 0.0, 1.0/4.0},
{ 0.5, 0.0, 0.0, 0.5},
{1.0/4.0, 0.0, 0.0, 3.0/4.0},
{ 0.0, 3.0/4.0, 1.0/4.0, 0.0},
{ 0.0, 0.5, 0.5, 0.0},
{ 0.0, 1.0/4.0, 3.0/4.0, 0.0},
{ 0.0, 3.0/4.0, 0.0, 1.0/4.0},
{ 0.0, 0.5, 0.0, 0.5},
{ 0.0, 1.0/4.0, 0.0, 3.0/4.0},
{ 0.0, 0.0, 3.0/4.0, 1.0/4.0},
{ 0.0, 0.0, 0.5, 0.5},
{ 0.0, 0.0, 1.0/4.0, 3.0/4.0},
{ 0.0, 1.0/2.0, 1.0/4.0, 1.0/4.0},
{ 0.0, 1.0/4.0, 1.0/2.0, 1.0/4.0},
{ 0.0, 1.0/4.0, 1.0/4.0, 1.0/2.0},
{1.0/2.0, 0.0, 1.0/4.0, 1.0/4.0},
{1.0/4.0, 0.0, 1.0/2.0, 1.0/4.0},
{1.0/4.0, 0.0, 1.0/4.0, 1.0/2.0},
{1.0/2.0, 1.0/4.0, 0.0, 1.0/4.0},
{1.0/4.0, 1.0/2.0, 0.0, 1.0/4.0},
{1.0/4.0, 1.0/4.0, 0.0, 1.0/2.0},
{1.0/2.0, 1.0/4.0, 1.0/4.0, 0.0},
{1.0/4.0, 1.0/2.0, 1.0/4.0, 0.0},
{1.0/4.0, 1.0/4.0, 1.0/2.0, 0.0},
{1.0/4.0, 1.0/4.0, 1.0/4.0, 1.0/4.0}};
/*--------------------------------------------------------------------------*/
/* basisfunction at vertex 0 */
/*--------------------------------------------------------------------------*/
static REAL phi4v0_3d(const REAL l[N_LAMBDA])
{
return((((32.0*l[0] - 48.0)*l[0] + 22.0)*l[0] - 3.0)*l[0]/3.0);
}
static const REAL *grd_phi4v0_3d(const REAL l[N_LAMBDA])
{
static REAL grd[N_LAMBDA];
grd[0] = ((128.0*l[0] - 144.0)*l[0] + 44.0)*l[0]/3.0 - 1.0;
return((const REAL *) grd);
}
static const REAL (*D2_phi4v0_3d(const REAL l[N_LAMBDA]))[N_LAMBDA]
{
static REAL D2[N_LAMBDA][N_LAMBDA] = {{0.0}};
D2[0][0] = (128.0*l[0] - 96.0)*l[0] + 44.0/3.0;
return((const REAL (*)[N_LAMBDA]) D2);
}
/*--------------------------------------------------------------------------*/
/* basisfunction at vertex 1 */
/*--------------------------------------------------------------------------*/
static REAL phi4v1_3d(const REAL l[N_LAMBDA])
{
return((((32.0*l[1] - 48.0)*l[1] + 22.0)*l[1] - 3.0)*l[1]/3.0);
}
static const REAL *grd_phi4v1_3d(const REAL l[N_LAMBDA])
{
static REAL grd[N_LAMBDA];
grd[1] = ((128.0*l[1] - 144.0)*l[1] + 44.0)*l[1]/3.0 - 1.0;
return((const REAL *) grd);
}
static const REAL (*D2_phi4v1_3d(const REAL l[N_LAMBDA]))[N_LAMBDA]
{
static REAL D2[N_LAMBDA][N_LAMBDA];
D2[1][1] = (128.0*l[1] - 96.0)*l[1] + 44.0/3.0;
return((const REAL (*)[N_LAMBDA]) D2);
}
/*--------------------------------------------------------------------------*/
/* basisfunction at vertex 2 */
/*--------------------------------------------------------------------------*/
static REAL phi4v2_3d(const REAL l[N_LAMBDA])
{
return((((32.0*l[2] - 48.0)*l[2] + 22.0)*l[2] - 3.0)*l[2]/3.0);
}
static const REAL *grd_phi4v2_3d(const REAL l[N_LAMBDA])
{
static REAL grd[N_LAMBDA];
grd[2] = ((128.0*l[2] - 144.0)*l[2] + 44.0)*l[2]/3.0 - 1.0;
return((const REAL *) grd);
}
static const REAL (*D2_phi4v2_3d(const REAL l[N_LAMBDA]))[N_LAMBDA]
{
static REAL D2[N_LAMBDA][N_LAMBDA];
D2[2][2] = (128.0*l[2] - 96.0)*l[2] + 44.0/3.0;
return((const REAL (*)[N_LAMBDA]) D2);
}
/*--------------------------------------------------------------------------*/
/* basisfunction at vertex 3 */
/*--------------------------------------------------------------------------*/
static REAL phi4v3_3d(const REAL l[N_LAMBDA])
{
return((((32.0*l[3] - 48.0)*l[3] + 22.0)*l[3] - 3.0)*l[3]/3.0);
}
static const REAL *grd_phi4v3_3d(const REAL l[N_LAMBDA])
{
static REAL grd[N_LAMBDA];
grd[3] = ((128.0*l[3] - 144.0)*l[3] + 44.0)*l[3]/3.0 - 1.0;
return((const REAL *) grd);
}
static const REAL (*D2_phi4v3_3d(const REAL l[N_LAMBDA]))[N_LAMBDA]
{
static REAL D2[N_LAMBDA][N_LAMBDA];
D2[3][3] = (128.0*l[3] - 96.0)*l[3] + 44.0/3.0;
return((const REAL (*)[N_LAMBDA]) D2);
}
/*--------------------------------------------------------------------------*/
/* basisfunction 0 at edge 0 */
/*--------------------------------------------------------------------------*/
static REAL phi4e00_3d(const REAL l[N_LAMBDA])
{
return(((128.0*l[0] - 96.0)*l[0] + 16.0)*l[0]*l[1]/3.0);
}
static const REAL *grd_phi4e00_3d(const REAL l[N_LAMBDA])
{
static REAL grd[N_LAMBDA];
grd[0] = ((128*l[0] - 64.0)*l[0] + 16.0/3.0)*l[1];
grd[1] = ((128*l[0] - 96.0)*l[0] + 16.0)*l[0]/3.0;
return((const REAL *) grd);
}
static const REAL (*D2_phi4e00_3d(const REAL l[N_LAMBDA]))[N_LAMBDA]
{
static REAL D2[N_LAMBDA][N_LAMBDA];
D2[0][0] = (256.0*l[0] - 64.0)*l[1];
D2[0][1] = D2[1][0] = (128.0*l[0] - 64.0)*l[0] + 16.0/3.0;
return((const REAL (*)[N_LAMBDA]) D2);
}
/*--------------------------------------------------------------------------*/
/* basisfunction 1 at edge 0 */
/*--------------------------------------------------------------------------*/
static REAL phi4e01_3d(const REAL l[N_LAMBDA])
{
return((4.0*l[0] - 1.0)*l[0]*(4.0*l[1] - 1.0)*l[1]*4.0);
}
static const REAL *grd_phi4e01_3d(const REAL l[N_LAMBDA])
{
static REAL grd[
adaptive finite element method （269个子文件） 
configure.ac 9KB Makefile.am 4KB Makefile.am 3KB Makefile.am 2KB Makefile.am 1KB Makefile.am 1KB Makefile.am 1KB Makefile.am 1KB Makefile.am 938B Makefile.am 693B Makefile.am 528B Makefile.am 528B Makefile.am 528B Makefile.am 501B Makefile.am 501B Makefile.am 501B Makefile.am 301B Makefile.am 206B Makefile.am 203B Makefile.am 194B Makefile.am 133B Makefile.am 14B Makefile.am 14B Makefile.am 13B macro-big.amc 632B macro.amc 631B macro_3d.amc 310B halfsphere.amc 295B macro.amc 292B macro-44circ.amc 245B macro-neum.amc 242B macro-circ.amc 238B macro-big.amc 230B macro-34circ.amc 225B macro-1.amc 204B macro-2.amc 194B macro-big_2d.amc 159B macro.amc 155B macro-big.amc 149B AUTHORS 874B AUTHORS 673B AUTHORS 609B AUTHORS 187B lagrange_4_3d.c 141KB dof_admin.c 93KB lagrange_3_3d.c 69KB dxtools.c 68KB numint.c 64KB memory.c 64KB assemble.c 59KB assemble_dowb.c 57KB lagrange_4_2d.c 55KB macro.c 50KB graphXO_2d.c 46KB submesh.c 43KB write_mesh_gmv.c 42KB memory_3d.c 42KB lagrange_3_2d.c 40KB eval.c 37KB HB_precon.c 37KB estimator_dowb.c 32KB estimator.c 32KB adapt.c 31KB lagrange_2_3d.c 30KB gltools.c 30KB submesh_3d.c 30KB getopt.c 30KB parametric_2d.c 29KB read_mesh.c 28KB error.c 28KB refine_3d.c 28KB lagrange_4_1d.c 27KB heat.c 27KB alberta_meshtv.c 27KB parametric_3d.c 26KB lagrange_2_2d.c 26KB lagrange_3_1d.c 25KB lagrange_2_1d.c 23KB traverse_nr.c 23KB coarsen_3d.c 23KB write_mesh.c 22KB MG_s2.c 22KB parameters.c 21KB memory_2d.c 21KB submesh_2d.c 21KB disc_lagrange_2_3d.c 20KB MG_s1.c 20KB slip_bound.c 19KB mesh_interface_common.c 19KB parametric_1d.c 18KB coarsen.c 18KB lagrange_1_1d.c 18KB refine.c 18KB macro_2d.c 17KB l2scp.c 17KB bas_fct.c 17KB ellipt2.c 17KB graphXO.c 17KB gmres.c 17KB level_3d.c 16KB共 269 条
- 1
- 2
- 3
资源评论
