Comments (7)
jhale
commented on August 16, 2024
This project sounds very interesting and it would be good if you could keep us informed on how things go as it proceeds.
Regarding your question, FEniCS 2019.1.0 is the legacy version and it does not support mixed/multi-precision operations.
We have built DOLFINx with mixed and multi-precision arithmetic in mind, but we haven't yet explored the potential in terms of implementing real algorithms; for your project I would strongly recommend using the C++ interface to DOLFINx where you have explicit control over the types of both the expansion coefficients and geometry information through templating, rather than Python.
Additionally, we've never tried quad precision, only double and lower precisions.
from dolfinx.
jhale
commented on August 16, 2024
C++ doesn't have a quad precision implementation in its standard library. So you will immediately run into issues with our use of the type concepts std::floating_point and dolfinx::scalar https://github.com/FEniCS/dolfinx/blob/main/cpp/dolfinx/common/types.h#L19
Are you sure you need quad precision? Typically you'd be looking to use doubles as your highest precision and then improve performance through clever use of lower precisions.
from dolfinx.
RamakrishnanChandran
commented on August 16, 2024
Dear Mr. Hale,
Yeah, I will check how I can use the C++ interface to create a driver for
Data extraction for my experiments.
In case you are aware of any tutorials. Could you please give some
suggestions for any tutorial which uses C++ interface directly.
Thanks and regards,
Ramakrishnan C.
…On Fri, Dec 8, 2023 at 11:55 AM Jack S. Hale ***@***.***> wrote:
This project sounds very interesting and it would be good if you could
keep us informed on how things go as it proceeds.
Regarding your question, FEniCS 2019.1.0 is the legacy version and it does
not support mixed/multi-precision operations.
We have built DOLFINx with mixed and multi-precision arithmetic in mind,
but we haven't yet explored the potential in terms of implementing real
algorithms; for your project I would strongly recommend using the C++
interface to DOLFINx where you have explicit control over the types of both
the expansion coefficients and geometry information through templating,
rather than Python.
Additionally, we've never tried quad precision, only lower precisions.
—
Reply to this email directly, view it on GitHub
<#2926 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AZDGVRD2KYLVXZFV3VJACPLYILW2RAVCNFSM6AAAAABAMN3QBWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBWHE2TQNJZGE>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
from dolfinx.
RamakrishnanChandran
commented on August 16, 2024
Dear Mr. Hale,
Yeah, that is true. The AMG solver (part of GaspiLS, a linear solver
library) which is experimented uses Double and single precision in IR
algorithmic setup.
However, f the Matrix and RHS data from FENICS and L2 error norm
measurement if done in higher precision could help in producing results
which are more reliable. That's the reason, I'm looking for Quad precision.
Anyway, the use of concepts coupled with no quad precision in std library
means even if I use C++ interface of FENICSx for Data extraction, I can't
build the code for quad precision right?
Thanks and Regards,
Ramakrishnan C.
…On Fri, Dec 8, 2023 at 12:01 PM Jack S. Hale ***@***.***> wrote:
C++ doesn't have a quad precision implementation in its standard library.
So you will immediately run into issues with our use of the type concepts
std::floating_point and dolfinx::scalar
https://github.com/FEniCS/dolfinx/blob/main/cpp/dolfinx/common/types.h#L19
Are you sure you need quad precision? Typically you'd be looking to use
doubles as your highest precision and then improve performance through
clever use of lower precisions.
—
Reply to this email directly, view it on GitHub
<#2926 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AZDGVREDV2BSTA6OVRWWYUTYILXQ5AVCNFSM6AAAAABAMN3QBWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBWHE3DMNRSGQ>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
from dolfinx.
jhale
commented on August 16, 2024
You can make a fork and modify the concepts to accept the types you need.
from dolfinx.
RamakrishnanChandran
commented on August 16, 2024
Dear Mr. Hale,
Thanks for your hint. I will try to work on that guidelines and update you.
Thanks and Regards,
Ramakrishnan C.
…On Fri, 8 Dec, 2023, 9:49 pm Jack S. Hale, ***@***.***> wrote:
You can make a fork and modify the concepts to accept the types you need.
—
Reply to this email directly, view it on GitHub
<#2926 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AZDGVRFHFCNJY7LC733O5CDYIN4MVAVCNFSM6AAAAABAMN3QBWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBXHAZDMOJRG4>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
from dolfinx.
garth-wells
commented on August 16, 2024
You could try with C++23 fixed width floating-point types (https://en.cppreference.com/w/cpp/types/floating-point). I don't know where compiler support is at.
from dolfinx.
Related Issues (20)
- `eval` is wrong for symmetric tensor elements
- Interpolation of non-symmetric tensor into a symmetric space
- Missing dof_indices method of DirichletBC in python HOT 1
- Add check of element hash
- Interpolation to/from sub-meshes causes parallel deadlock HOT 3
- [BUG]: Wrong dispatchment of field values in the mesh with VTK and VTX writers using a mixed space
- [BUG]: Conda test fails 0.8.0 HOT 3
- Add tests of LinearProblem and Nonlinearproblem
- Change C++ logging library to `spdlog` or `glog` HOT 6
- Support building on Windows
- [BUG]: dolfin / dolfinx /fenics / fenicsx not in debian testing repositories. HOT 4
- Support `dolfinx::fem::Expression` with co-dim 0 submeshes
- [BUG]: `create_entity_permutations` occasionally segfaults for prism meshes HOT 3
- Remove Fides file output HOT 5
- `dolfinx.fem.Form.integral_types` nanobinding is incorrect.
- interior_facet_assembly causes deadlock if there are no interior facets on the process HOT 2
- Test needed for dof permutations on prism and pyramid cells
- Read mixed topology mesh from file HOT 3
- Add mixed topology functionality to mesh partitioning HOT 1
- Write mixed topology mesh to file
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from dolfinx.