{"id":23977451,"url":"https://github.com/thecomputekid/ehl-point-stationary","last_synced_at":"2025-04-14T01:07:57.007Z","repository":{"id":151366007,"uuid":"402899586","full_name":"theComputeKid/ehl-point-stationary","owner":"theComputeKid","description":"Elastohydrodynamic Lubrication Point Contact Solver for MATLAB","archived":false,"fork":false,"pushed_at":"2024-04-26T15:39:04.000Z","size":73,"stargazers_count":25,"open_issues_count":0,"forks_count":6,"subscribers_count":1,"default_branch":"main","last_synced_at":"2025-04-14T01:07:39.144Z","etag":null,"topics":[],"latest_commit_sha":null,"homepage":"","language":"MATLAB","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"other","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/theComputeKid.png","metadata":{"files":{"readme":"README","changelog":null,"contributing":null,"funding":null,"license":"LICENSE","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null,"roadmap":null,"authors":null,"dei":null}},"created_at":"2021-09-03T21:17:33.000Z","updated_at":"2025-03-10T14:06:05.000Z","dependencies_parsed_at":null,"dependency_job_id":"682cf257-4725-4175-a2bf-a2070ed351c7","html_url":"https://github.com/theComputeKid/ehl-point-stationary","commit_stats":null,"previous_names":[],"tags_count":2,"template":false,"template_full_name":null,"repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/theComputeKid%2Fehl-point-stationary","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/theComputeKid%2Fehl-point-stationary/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/theComputeKid%2Fehl-point-stationary/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/theComputeKid%2Fehl-point-stationary/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/theComputeKid","download_url":"https://codeload.github.com/theComputeKid/ehl-point-stationary/tar.gz/refs/heads/main","host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":248804803,"owners_count":21164134,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2022-07-04T15:15:14.044Z","host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":[],"created_at":"2025-01-07T07:37:49.158Z","updated_at":"2025-04-14T01:07:56.976Z","avatar_url":"https://github.com/theComputeKid.png","language":"MATLAB","funding_links":[],"categories":[],"sub_categories":[],"readme":"-----------------------------------------------------------------\n\u003e\u003e Elastohydrodynamic Lubrication Point Contact Solver for MATLAB\n-----------------------------------------------------------------\n\nBased on the theory presented in Venner and Lubrecht's book \"Multilevel Methods in Lubrication\", with the following modifications:\n- Fourier Transforms are used to calculate deformation instead of the Multi-Level Multi-Integration mentioned in the book\n- Parallel line solvers are provided.\n- An option to execute the solver using the GPU is added\n\nThere are multiple execution strategies that can be used:\n- cpu_seq : The CPU will solve and progress in the domain line by line, using a banded matrix solver with 6 bands. This is faithful to the algorithm mentioned in the book.\n- cpu_par : The CPU solves in parallel, solving multiple lines at the same time, provided that the data required does not overlap (similar to a red-black algorithm, except for lines). There is some overhead to launching a parallel pool. Requires Parallel Computing Toolbox.\n- gpu_seq : The GPU will solve and progress in the domain line by line, using a banded matrix solver with 5 bands. Only for debugging purposes because this is slow.\n- gpu_par : The GPU solves in parallel, solving multiple lines at the same time, provided that the data required does not overlap (similar to a red-black algorithm, except for lines). Performance relative to the CPU will scale with the number of lines that can be solved in parallel. \n\nExamples are provided in the +examples folder. To run them, you must be in the project root folder (i.e. where this README is). From there, you can run examples from the examples namespace. E.g. from the command window in the root folder:\n\nCPU Examples:\n\u003e\u003e examples.ex1 : Basic Example with Result Plots\n\u003e\u003e examples.ex2 : Basic Example with Textual Verbosity\n\u003e\u003e examples.ex3 : Non-Isoviscous Fluids with Graphical Verbosity\n\u003e\u003e examples.ex4 : Non-Isoviscous Fluids with Single Precision\n\u003e\u003e examples.ex5 : Uses the parallel line solver, but single grid only\n\u003e\u003e examples.ex6 : Uses the parallel line solver, but with multiple grids\n\nGPU Examples (require MATLAB with Parallel Computing Toolbox and a supported NVIDIA GPU):\n\u003e\u003e examples.ex7 : Sequential GPU solver\n\u003e\u003e examples.ex8 : Parallel GPU line solver, but with single grid\n\u003e\u003e examples.ex9 : Parallel GPU line solver, but with multiple grids\n\nBenchmarks are provided to show the difference in performance between the execution modes. They are located in the +bench folder:\n\u003e\u003e bench.singleGrid : Single Grid benchmarks for 5 iterations.\n\u003e\u003e bench.multiGrid : Multi-Grid benchmarks for 1 iteration.\n\nNote on GPU support: GPU support requires mex building, as a batch pentadiagonal solver is not part of MATLAB (yet). The mex file uses C++17, which needs a (recent) version of MATLAB that uses a version of CUDA that supports it. It is built automatically as part of the solution process when required. No user intervention is required.\n\nTested on MATLAB R2023b. In MATLAB 2024a, the banded solver now only accepts a square full matrix, which requires changes to the project. Previously, supplying only the bands of the matrix was sufficient.\n\nVariable naming conventions are kept similar to those used in \"Multilevel Methods in Lubrication\".\n\nReferences:\nVenner, C.H., Lubrecht, A.A. (Eds.), 2000. Chapter 6 - ElastoHydrodynamic Lubrication, in: Multilevel Methods in Lubrication, Tribology Series. Elsevier, pp. 179–223. https://doi.org/10.1016/S0167-8922(00)80010-3\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fthecomputekid%2Fehl-point-stationary","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fthecomputekid%2Fehl-point-stationary","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fthecomputekid%2Fehl-point-stationary/lists"}