Metric tensor perf geom

cases/condif2/condif.c

@@ -0,0 +1,35 @@
+//--------------------------------------------------------------
+// Solution function for convection-diffusion
+//
+// To compile run, for example:
+// mpicc -shared -o libfun.so -fPIC condif.c
+// gcc   -shared -o libfun.so -fPIC condif.c
+//--------------------------------------------------------------
+
+#include <math.h>
+#include <stdio.h>
+
+
+double sol_phi(double *coord, int dim, double time)
+{
+   return 0.0;
+}
+
+void advection(double *coord, int dim, double time, double *adv, int vdim)
+{
+   double r = sqrt(coord[0]*coord[0] + coord[1]*coord[1]);
+
+   adv[0] = -coord[0]/r;//-sqrt(2.0)/2.0;
+   adv[1] = -coord[1]/r;//-sqrt(2.0)/2.0;
+}
+
+double force(double *coord, int dim, double time)
+{
+   return 1.0;
+}
+
+double mu(double *coord, int dim, double time)
+{
+   return 0.00001;
+}
+

cases/condif2/condif_ms.c

@@ -0,0 +1,56 @@
+//--------------------------------------------------------------
+// Solution function for convection-diffusion -- manufactured solution
+//
+// To compile run, for example:
+// mpicc -shared -o libfun.so -fPIC condif_ms.c
+// gcc   -shared -o libfun.so -fPIC condif_ms.c
+//--------------------------------------------------------------
+
+# TO BE MODIFIED  ???
+bla
+
+#include <math.h>
+#include <stdio.h>
+
+void advection(double *x, int dim, double t, double *adv, int vdim)
+{
+   adv[0] = sqrt(2.0)/2.0;
+   adv[1] = sqrt(2.0)/2.0;
+}
+
+double mu(double *x, int dim, double t)
+{
+   return 0.01;
+}
+
+double sol_phi(double *x, int dim, double t)
+{
+   return sin(M_PI*x[0])*sin(M_PI*x[1]);
+}
+
+void grad_phi(double *x, int dim, double t, double *grad, int vdim)
+{
+   grad[0] = M_PI*cos(M_PI*x[0])*sin(M_PI*x[1]);
+   grad[1] = M_PI*sin(M_PI*x[0])*cos(M_PI*x[1]);
+}
+
+double lap_phi(double *x, int dim, double t)
+{
+   return -2*M_PI*M_PI*sin(M_PI*x[0])*sin(M_PI*x[1]);
+}
+
+double force(double *x, int dim, double t)
+{
+   double a[dim];
+   double grad[dim];
+   advection(x, dim, t, a, dim);
+   grad_phi(x, dim, t, grad, dim);
+
+   double res = - mu(x, dim, t)*lap_phi(x, dim, t);
+   for (int i = 0; i < dim; i++)
+   {
+      res += a[i]*grad[i];
+   }
+
+   return res;
+}

cases/condif2/quart-circle-eps.mesh

@@ -0,0 +1,64 @@
+MFEM NURBS mesh v1.0
+
+#
+# MFEM Geometry Types (see fem/geom.hpp):
+#
+# SEGMENT     = 1
+# SQUARE      = 3
+# CUBE        = 5
+#
+
+dimension
+2
+
+elements
+1
+1 3 0 1 2 3
+
+boundary
+4
+1 1 0 1
+3 1 2 3
+4 1 3 0
+2 1 1 2
+
+edges
+4
+0 0 1
+0 3 2
+1 0 3
+1 1 2
+
+vertices
+4
+
+knotvectors
+2
+2 3 0 0 0 1 1 1
+2 3 0 0 0 1 1 1
+
+weights
+1
+1
+1
+1
+0.707106781
+0.707106781
+1
+1
+0.707106781
+
+FiniteElementSpace
+FiniteElementCollection: NURBS2
+VDim: 2
+Ordering: 1
+
+0 0.001
+0.001 0
+1 0
+0 1
+0.001 0.001
+1 1
+0 0.001
+0.001 0
+0.001 0.001

cases/condif2/quart-circle.mesh

@@ -0,0 +1,64 @@
+MFEM NURBS mesh v1.0
+
+#
+# MFEM Geometry Types (see fem/geom.hpp):
+#
+# SEGMENT     = 1
+# SQUARE      = 3
+# CUBE        = 5
+#
+
+dimension
+2
+
+elements
+1
+1 3 0 1 2 3
+
+boundary
+4
+1 1 0 1
+3 1 2 3
+4 1 3 0
+2 1 1 2
+
+edges
+4
+0 0 1
+0 3 2
+1 0 3
+1 1 2
+
+vertices
+4
+
+knotvectors
+2
+2 3 0 0 0 1 1 1
+2 3 0 0 0 1 1 1
+
+weights
+1
+1
+1
+1
+0.707106781
+0.707106781
+1
+1
+0.707106781
+
+FiniteElementSpace
+FiniteElementCollection: NURBS2
+VDim: 2
+Ordering: 1
+
+0 0.00
+0.00 0
+1 0
+0 1
+0.00 0.00
+1 1
+0 0.5
+0.5 0
+0.5 0.5

cases/condif2/run-nurbs

@@ -0,0 +1,9 @@
+#!/bin/bash
+gcc -shared -o libfun.so -fPIC condif.c
+
+exe=../../build/supg/supg
+mesh=quart-circle.mesh
+
+rm output* log*
+
+mpirun -n 1 $exe -m $mesh -o 2 -r 3  --strong-bdr "1 2 3 4"   -l libfun.so | tee log0

cases/condif3/condif.c

@@ -0,0 +1,35 @@
+//--------------------------------------------------------------
+// Solution function for convection-diffusion
+//
+// To compile run, for example:
+// mpicc -shared -o libfun.so -fPIC condif.c
+// gcc   -shared -o libfun.so -fPIC condif.c
+//--------------------------------------------------------------
+
+#include <math.h>
+#include <stdio.h>
+
+
+double sol_phi(double *coord, int dim, double time)
+{
+   double r = sqrt(coord[0]*coord[0] + coord[1]*coord[1]);
+
+   return 0.5*(1.0 - cos(2*M_PI*fmin(1.0, r)));
+}
+
+void advection(double *coord, int dim, double time, double *adv, int vdim)
+{
+   adv[0] = coord[1];
+   adv[1] = -coord[0];
+}
+
+double force(double *coord, int dim, double time)
+{
+   return 0.0;
+}
+
+double mu(double *coord, int dim, double time)
+{
+   return 0.0;
+}
+

cases/condif3/condif_ms.c

@@ -0,0 +1,55 @@
+//--------------------------------------------------------------
+// Solution function for convection-diffusion -- manufactured solution
+//
+// To compile run, for example:
+// mpicc -shared -o libfun.so -fPIC condif_ms.c
+// gcc   -shared -o libfun.so -fPIC condif_ms.c
+//--------------------------------------------------------------
+
+## ??
+
+#include <math.h>
+#include <stdio.h>
+
+void advection(double *x, int dim, double t, double *adv, int vdim)
+{
+   adv[0] = sqrt(2.0)/2.0;
+   adv[1] = sqrt(2.0)/2.0;
+}
+
+double mu(double *x, int dim, double t)
+{
+   return 0.01;
+}
+
+double sol_phi(double *x, int dim, double t)
+{
+   return sin(M_PI*x[0])*sin(M_PI*x[1]);
+}
+
+void grad_phi(double *x, int dim, double t, double *grad, int vdim)
+{
+   grad[0] = M_PI*cos(M_PI*x[0])*sin(M_PI*x[1]);
+   grad[1] = M_PI*sin(M_PI*x[0])*cos(M_PI*x[1]);
+}
+
+double lap_phi(double *x, int dim, double t)
+{
+   return -2*M_PI*M_PI*sin(M_PI*x[0])*sin(M_PI*x[1]);
+}
+
+double force(double *x, int dim, double t)
+{
+   double a[dim];
+   double grad[dim];
+   advection(x, dim, t, a, dim);
+   grad_phi(x, dim, t, grad, dim);
+
+   double res = - mu(x, dim, t)*lap_phi(x, dim, t);
+   for (int i = 0; i < dim; i++)
+   {
+      res += a[i]*grad[i];
+   }
+
+   return res;
+}

cases/condif3/run-nurbs

@@ -0,0 +1,10 @@
+#!/bin/bash
+gcc -shared -o libfun.so -fPIC condif.c
+
+exe=../../build/supg/supg
+mesh=square.mesh
+mesh=square-slit.mesh
+
+rm output* log*
+
+mpirun -n 1 $exe -m $mesh -o 2 -r 3  --strong-bdr "1 2"   -l libfun.so | tee log0