Skip to content
Snippets Groups Projects

Feature matrix free solution of linear PDEs

Merged Feature matrix free solution of linear PDEs
feature/matrix-free into master
Merged René Heß requested to merge feature/matrix-free into master

This implements matrix free solution of linear PDEs. Some notes:

  1. Right now the code simply uses BiCGSTAB from istl from istl for solving the linear systems. A more subtle solver setup should be postponed until the istl-ini factory arrived at PDELab.
  2. Poisson DG doesn't converge. Tests with the convectiondiffusiondg lop from PDELab show that the hand written lop has the same behaviour. Because of that I turned off DG tests.
  3. Right now nonlinear PDE's are not covered since PDELab doesn't seem to support nonlinear matrix free solutions. We probably have to adjust Newton first.

This merge request is related to #26 (closed) and I will update the issue accordingly.

Merge request reports

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
  • Dominic Kempf
    Dominic Kempf @dominic ·
    Owner

    filetag should be driver. I can only guess, but I think this will result in the include also being present in the operator (doesnt hurt though).

    By Dominic Kempf on 2016-08-17T15:18:27 (imported from GitLab project)

  • Dominic Kempf
    Dominic Kempf @dominic ·
    Owner

    As said this morning, I would just use Coefficient(element, count=1) here to keep the amount of stuff in dune.perftool.ufl.execution minimal. A user would never use ApplyFunction, so there is no need to expose it to him.

    By Dominic Kempf on 2016-08-17T15:26:04 (imported from GitLab project)

  • Dominic Kempf
    Dominic Kempf @dominic ·
    Owner

    This could/(should?) be redesigned to only take 2 arguments: argnumber and coeffcount, which default to None. Only those modified terminals are returned that match the specified condition.

    By Dominic Kempf on 2016-08-17T15:30:06 (imported from GitLab project)

  • Dominic Kempf
    Dominic Kempf @dominic ·
    Owner

    Looks great! I was looking at it very thorough and couldn't find more to complain :D

    By Dominic Kempf on 2016-08-17T15:31:47 (imported from GitLab project)

  • René Heß Added 1 commit:

    Added 1 commit:

    • 6917da50 - Some small improvements and bugfixes
  • René Heß
    René Heß @rhess ·
    Author Maintainer

    This should fix your complains. Thanks for looking at the code.

  • René Heß
    René Heß @rhess ·
    Author Maintainer

    Totally forgot about this one after lunch ;).

  • Dominic Kempf
    Dominic Kempf @dominic ·
    Owner

    Great! Merging.

    By Dominic Kempf on 2016-08-18T08:35:58 (imported from GitLab project)

  • Dominic Kempf Status changed to merged

    Status changed to merged

    By Dominic Kempf on 2016-08-18T08:36:04 (imported from GitLab project)

  • Dominic Kempf mentioned in commit cbe4bf96

    mentioned in commit cbe4bf96

    By Dominic Kempf on 2016-08-18T08:36:04 (imported from GitLab project)

Please register or sign in to reply