Skip to content

Conversation

@ESeNonFossiIo
Copy link
Contributor

@ESeNonFossiIo ESeNonFossiIo commented Jun 29, 2016

Improvements:

  • 3D
  • remove nitsche
  • two tests (just to check compatibility with NS in the case d==0, for further tests some work is needed - see below)
    • trivial -> const solution
    • not trivial but d==0
    • not trivial with d!=0

Missing:

  • mappingQEulerian for bc, output, error. (exception: fev_cache)
    • restart -> problems with bc

if (solution.size() == 0)
return StaticMappingQ1<dim,spacedim>::mapping;

if (!mapping)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

why degree == 1 instead of this->get_fe().degree?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

no reason… I change it! :)

On 29 Jun 2016, at 19:23, Luca Heltai [email protected] wrote:

this->get_fe().degree


// jacobian of ALE transformation:
auto J_ale = determinant(F);

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

ocio che auto + sacado può avere comportamenti imprevedibili

)*JxW[quad];
- div_u * p_test;
// - scalar_product(J_ale * F_inv * u, grad_p_test); // TODO: muovere la divergenza a sinistra! se no c'é anche l'integrale sul bordo.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@luca-heltai Is it possible to evaluate scalar_product(div(J_ale * F_inv * u), p_test)?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants