Commit 756dbd4d authored by Michel Henry's avatar Michel Henry
Browse files

wip

parent 418d866e
Pipeline #9710 failed with stages
in 60 minutes and 39 seconds
......@@ -89,7 +89,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_triangle_p1;
fe->df = fe_df_triangle_p1;
fe->type = "triangle_p1";
}
else if (strcmp(etype,"triangle_p1dg") == 0) {
fe->nlocal = 3;
......@@ -99,7 +98,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_triangle_p1;
fe->df = fe_df_triangle_p1;
fe->type = "triangle_p1dg";
}
else if (strcmp(etype,"triangle_p2") == 0) {
fe->nlocal = 6;
......@@ -109,7 +107,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_triangle_p2;
fe->df = fe_df_triangle_p2;
fe->type = "triangle_p2";
}
else if (strcmp(etype,"line_p1") == 0) {
fe->nlocal = 2;
......@@ -119,7 +116,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_line_p1;
fe->df = fe_df_line_p1;
fe->type = "line_p1";
}
else if (strcmp(etype,"line_p2") == 0) {
fe->nlocal = 3;
......@@ -129,7 +125,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_line_p2;
fe->df = fe_df_line_p2;
fe->type = "line_p2";
}
else if (strcmp(etype,"tetrahedron_p1") == 0) {
fe->nlocal = 4;
......@@ -139,7 +134,6 @@ FEElement *fe_element_new(const char *etype) {
fe->n[3] = 0;
fe->f = fe_f_tetrahedron_p1;
fe->df = fe_df_tetrahedron_p1;
fe->type = "tetrahedron_p1";
}
else {
......
......@@ -10,7 +10,6 @@ struct FEElementStruct{
int n[DIMENSION+1];
void (*f)(const double *xi, double *f);
void (*df)(const double *xi, const double dxidx[D][D], double f[][D]);
const char *type;
};
FEElement *fe_element_new(const char *etype);
......
......@@ -110,8 +110,7 @@ class DieSwell(unittest.TestCase):
#
swelling = 0
while t < tEnd :
# time_integration.iterate(fluid,None,dt)
fluid.implicit_euler(dt, check_residual_norm=5e-4)
time_integration.iterate(fluid, None, dt, check_residual_norm=5e-4)
if t >= 50.0:
fs.iterate(dt)
swelling = fluid.coordinates()[fs.fs_nodes[-1],1]/H
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment