Commit abbd5df6 authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

one less but still not ok

parent e81d56f0
Pipeline #9678 passed with stages
in 3 minutes and 5 seconds
...@@ -513,12 +513,12 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, const double ...@@ -513,12 +513,12 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, const double
double s[fields->n], sold[fields->n], ds[fields->n][D], dsold[fields->n][D]; double s[fields->n], sold[fields->n], ds[fields->n][D], dsold[fields->n][D];
double sf[fields->local_size], dsf[fields->local_size][D]; double sf[fields->local_size], dsf[fields->local_size][D];
fe_fields_f(fields, xiel, sf); fe_fields_f(fields, xiel, sf);
fe_fields_df(fields, QP[iqp], dxidx, dsf); fe_fields_df(fields, xiel, dxidx, dsf);
fe_fields_eval_grad_sf(fields, mesh, iel, sf, dsf, problem->solution, s, ds); fe_fields_eval_grad_sf(fields, mesh, iel, sf, dsf, problem->solution, s, ds);
fe_fields_eval_grad_sf(fields, mesh, iel, sf, dsf, solution_old, sold, dsold); fe_fields_eval_grad_sf(fields, mesh, iel, sf, dsf, solution_old, sold, dsold);
double dc[D], c, sfporosity[problem->field_porosity->local_size], dsfporosity[problem->field_porosity->local_size][D] ; double dc[D], c, sfporosity[problem->field_porosity->local_size], dsfporosity[problem->field_porosity->local_size][D] ;
fe_fields_f(problem->field_porosity, QP[iqp], sfporosity); fe_fields_f(problem->field_porosity, xiel, sfporosity);
fe_fields_df(problem->field_porosity, QP[iqp], dxidx, dsfporosity); fe_fields_df(problem->field_porosity, xiel, dxidx, dsfporosity);
fe_fields_eval_grad_sf(problem->field_porosity, mesh, iel, sfporosity, dsfporosity, problem->porosity, &c, &dc); fe_fields_eval_grad_sf(problem->field_porosity, mesh, iel, sfporosity, dsfporosity, problem->porosity, &c, &dc);
double f0[fields->n], f00[fields->n*fields->n],f01[fields->n*fields->n][D]; double f0[fields->n], f00[fields->n*fields->n],f01[fields->n*fields->n][D];
for (int i = 0; i < fields->n; ++i) { for (int i = 0; i < fields->n; ++i) {
...@@ -535,12 +535,12 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, const double ...@@ -535,12 +535,12 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, const double
double a = 0; double a = 0;
if (problem->n_fluids==2) { if (problem->n_fluids==2) {
double sfconcentration[problem->field_concentration->local_size]; double sfconcentration[problem->field_concentration->local_size];
fe_fields_f(problem->field_concentration, QP[iqp], sfconcentration); fe_fields_f(problem->field_concentration, xiel, sfconcentration);
fe_fields_eval_sf(problem->field_concentration, mesh, iel, sfconcentration, problem->concentration, &a); fe_fields_eval_sf(problem->field_concentration, mesh, iel, sfconcentration, problem->concentration, &a);
} }
const FEElement *mesh_element = problem->mesh->element; const FEElement *mesh_element = problem->mesh->element;
double meshf[mesh_element->nlocal]; double meshf[mesh_element->nlocal];
mesh_element->f(QP[iqp], meshf); mesh_element->f(xiel, meshf);
double umesh[D] = {0}; double umesh[D] = {0};
for (int i = 0; i < mesh_element->nlocal; ++i) { for (int i = 0; i < mesh_element->nlocal; ++i) {
for (int j=0; j<D; ++j) { for (int j=0; j<D; ++j) {
......
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