Week June 20 -- June 25

This week I implemented a Dirac solver in 1D, I have it in my private branch and I'll get a permission at the end of the summer to get it opensourced. Originally it didn't work at all, but then I've found a bug and it seems to be working now.

Before it I've implemented an example for the y''+k^2*y=0 equation in hermes1d by writing it as two coupled first order equations and solving the eigenproblem for "k".

I also started to implement a Function class that represents a function in 1D using Fekete points and I plan to implement orthogonal projections and adaptive algorithms for refining the mesh, so that I can play with adaptivity in 1D.

Besides that I studied couple articles about using b-splines to solve the radial Dirac equation and also about using variational formulation for FEM.

I helped fix hermes_common and did lots of administration at the lab (10 online courses and other things). We also got an awesome NIF tour.

Hermes1D patches:

Ondrej Certik (26):
      Remove the autogenerated Makefile
      hermes_common updated
      Dirac solver added
      Plot only 10 eigenfunctions
      Fix the number of equations bug
      Generate the Gauss Lobatto points
      Add the autogenerated file
      Don't use strings as dict keys
      Fekete points manipulation implemented
      fekete: solve the coefficients, tests pass
      More tests added
      system_sin_eigen example added
      Add a comment about the Lagrange interpolation
      Function.project_onto() implemented
      Function.plot() implemented
      Comparisons implemented
      Fix precision problems
      Improve tests
      Start implementing the adaptivity
      fekete: add a debug code
      Removing the Dirac solver for now
      Use smaller basis, while still getting reasonable results
      Add harmonic oscillator option
      hermes_common updated
      Use jsplot if it is available
      Better demo

hermes_common patches:

Ondrej Certik (6):
      Revert "fixed pxd"
      Update the generated .cpp file
      Increase the precision in tests to 10^-10
      Add tests for solve_linear_system_dense_lu()
      Enable scipy tests
      Test for solve_linear_system_cg() added