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.
Ondrej Certik (26):
Remove the autogenerated Makefile
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
Fix precision problems
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
Use jsplot if it is available
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