Line data Source code
1 : // Copyright 2024, UChicago Argonne, LLC
2 : // All Rights Reserved
3 : // Software Name: NEML2 -- the New Engineering material Model Library, version 2
4 : // By: Argonne National Laboratory
5 : // OPEN SOURCE LICENSE (MIT)
6 : //
7 : // Permission is hereby granted, free of charge, to any person obtaining a copy
8 : // of this software and associated documentation files (the "Software"), to deal
9 : // in the Software without restriction, including without limitation the rights
10 : // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
11 : // copies of the Software, and to permit persons to whom the Software is
12 : // furnished to do so, subject to the following conditions:
13 : //
14 : // The above copyright notice and this permission notice shall be included in
15 : // all copies or substantial portions of the Software.
16 : //
17 : // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
18 : // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
19 : // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
20 : // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
21 : // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
22 : // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
23 : // THE SOFTWARE.
24 :
25 : #include "neml2/tensors/Transformable.h"
26 :
27 : #include "neml2/tensors/tensors.h"
28 :
29 : namespace neml2
30 : {
31 :
32 : R2
33 76 : transform_from_quaternion(const Quaternion & q)
34 : {
35 76 : return q.to_R2();
36 : }
37 :
38 : R2
39 2 : identity_transform(const TensorOptions & options)
40 : {
41 2 : return R2::identity(options);
42 : }
43 :
44 : R2
45 2 : proper_rotation_transform(const Rot & rot)
46 : {
47 2 : return rot.euler_rodrigues();
48 : }
49 :
50 : R2
51 3 : improper_rotation_transform(const Rot & rot)
52 : {
53 3 : Vec v = rot / rot.norm();
54 6 : return rot.euler_rodrigues() * (R2::identity(rot.options()) - 2 * v.outer(v));
55 3 : }
56 :
57 : R2
58 0 : reflection_transform(const Vec & v)
59 : {
60 0 : return R2::identity(v.options()) - 2 * v.outer(v);
61 : }
62 :
63 : R2
64 2 : inversion_transform(const TensorOptions & option)
65 : {
66 2 : return R2::fill(-1.0, option);
67 : }
68 :
69 : } // namespace neml2
|