Commit Graph

680 Commits

Author SHA1 Message Date
Jukka Aho be2da7dd0a moved continuum elasticity test 2016-02-06 15:28:32 +02:00
Jukka Aho d2a42c3c14 mass consistent normals 2016-02-05 22:26:07 +02:00
Jukka Aho 589370172f more integration rules 2016-02-05 22:22:40 +02:00
Jukka Aho 2f990e5eb1 better mesh 2016-02-05 22:22:14 +02:00
Jukka Aho 13aca13c5a posthooks for field and boundary assembly 2016-02-05 14:03:27 +02:00
Jukka Aho ff6abfc5e8 small bug fixes 2016-02-05 12:27:36 +02:00
Jukka Aho 73472b5e70 drop near nonzero values when creating sparse matrices. function to find nonzero rows from sparse matrix. 2016-02-05 11:32:50 +02:00
Jukka Aho 552b72100e matrices for dual formulation for elements 2016-02-05 11:32:09 +02:00
Jukka Aho 2ac577d3f6 functions to find nodes from dofs and vice versa 2016-02-05 11:31:23 +02:00
Jukka Aho 139eff7f55 update primary field to boundary problems also. 2016-02-04 17:32:51 +02:00
Jukka Aho 90c10d19b8 calculate rotated nodal vector 2016-02-04 07:28:19 +02:00
Jukka Aho 75ca704bec two test models 2016-02-04 06:55:08 +02:00
Jukka Aho ab741411e0 first attempts in solving overdetermined system in a more clever way 2016-02-04 06:54:08 +02:00
Jukka Aho 6652f48ea2 fixes 2016-02-03 20:39:03 +02:00
Jukka Aho 90319b2644 Major changes in data structures:
- Combined FieldAssembly and BoundaryAssembly to Assembly
- Combined FieldProblem and BoundaryProblem to Problem
- FieldProblem and BoundaryProblem are now abstract types
- Renamed stiffness_matrix, mass_matrix and force_vector to K, M, f for
  easier notation
- Problems are no more abstract types but concrete types, see
  elasticity.jl for example
- Combined linear_elasticity.jl and elasticity.jl
- Removed obsolete code directsolver.jl
- Almost all tests probably fail at this point
2016-02-03 06:49:42 +02:00
Jukka Aho 60046942e5 new style concrete problem 2016-02-02 22:22:21 +02:00
Jukka Aho fe98a75cbe nt formulation with gap calculation. 2016-02-02 18:00:17 +02:00
Jukka Aho 5550388dc1 mortar contact tests 2016-02-01 09:16:30 +02:00
Jukka Aho bbf3d2678c test of new solver 2016-02-01 09:15:42 +02:00
Jukka Aho f4a1bbcf3f dual basis formulation 2016-02-01 09:15:13 +02:00
Jukka Aho 43b9977dbe redesign of solver, directsolver will be obsolete 2016-02-01 09:14:43 +02:00
Jukka Aho c37250fe6d moving code 2016-02-01 09:14:08 +02:00
Jukka Aho 991a28ec5d element node mappings.. where to get them. 2016-02-01 09:13:34 +02:00
Jukka Aho 73b4609c0e rewrite solver 2016-02-01 09:13:07 +02:00
Jukka Aho 927cf788dc problem -> assembly 2016-02-01 09:12:41 +02:00
Jukka Aho be74499858 update! for several elements only using Dict 2016-02-01 09:11:33 +02:00
Jukka Aho 3052053cdc symmetric formulation 2016-02-01 09:10:44 +02:00
Jukka Aho 84f4c613a9 flag to control does assembly need updating 2016-02-01 09:10:12 +02:00
Jukka Aho 826db61417 still needs some fixing, gives wrong result 2016-01-02 21:44:45 +02:00
Jukka Aho 730dc10935 i think it's working now 2016-01-02 00:54:43 +02:00
Jukka Aho 5fae7c54fd solve total quantity u for now, set solve_residual=false in solver 2016-01-02 00:17:40 +02:00
Jukka Aho 225a1dc165 bug fix post -> pre 2016-01-02 00:17:01 +02:00
Jukka Aho a8e22ccbe1 add several pre/postprocessors with optional arguments, option to solve total or residual vector, function to check convergence (overridable) 2016-01-02 00:16:36 +02:00
Jukka Aho 74e9282697 its now possible to define several pre and postprocessors for assemblies with custom arguments; see test/test_mortar_2d_contact.jl for concrete example 2016-01-01 20:34:10 +02:00
Jukka Aho 01cbf2516f lets use FieldProblem(ElasticityProblem, ...) instead of FieldProblem{ElasticityProblem}(...) so that we can extend FieldProblem and BoundaryProblem without redefining everything in problems. 2016-01-01 18:52:40 +02:00
Jukka Aho ac54dbf2c8 change default DirichletProblem{StandardBasis} -> DirichletProblem 2016-01-01 18:51:34 +02:00
Jukka Aho 407b0c4849 ... 2016-01-01 17:57:24 +02:00
Jukka Aho ac45c31ce3 give convergence information 2016-01-01 17:56:52 +02:00
Jukka Aho 192953214f more general way to define preprocessors and postprocessors for linear solution 2016-01-01 16:56:17 +02:00
Jukka Aho 9d43635315 possibility to preprocess assembly 2016-01-01 16:55:41 +02:00
Jukka Aho e0091e3246 frictionless 2d contact tests, initial commit 2016-01-01 15:00:50 +02:00
Jukka Aho 52fd3f5258 update all kind of fields using update! function 2016-01-01 15:00:18 +02:00
Jukka Aho 91173dc78e asdf 2015-12-31 12:35:58 +02:00
Jukka Aho a586de6966 template 2015-12-31 12:35:45 +02:00
Jukka Aho b1e4a5b823 skip elements without master elements defined 2015-12-31 12:35:33 +02:00
Jukka Aho 0d11e82e37 element template 2015-12-31 12:35:05 +02:00
Jukka Aho 5c3d372795 progessing 2015-12-31 08:40:29 +02:00
Jukka Aho ffb51d251b calculate_nodal_vector: calculate for all dimensions using single factorization 2015-12-31 08:38:33 +02:00
Jukka Aho 482aca2aeb sparse rhs: calculate several right hand sides at once 2015-12-31 08:37:47 +02:00
Jukka Aho 6fb255db7e fix template 2015-12-31 07:41:46 +02:00