Mercurial > hg > orthanc-stone
view Resources/Computations/ComputeShearParameters.py @ 1314:9b126de2cde2 broker
Since the observer system now uses shared_ptr and many registrations are
done in the constructors, and since we cannot called shared_from_this() in
the constructors, it is mandatory to split construction from registration.
This has been done by making many ctors protected and replacing them
by factory methods that directly return shared_ptrs + added PostConstructor
method when base classes perform shared_from_this() calls too.
author | Benjamin Golinvaux <bgo@osimis.io> |
---|---|
date | Mon, 16 Mar 2020 11:19:50 +0100 |
parents | 46cb2eedc2e0 |
children |
line wrap: on
line source
#!/usr/bin/python from sympy import * import pprint init_printing(use_unicode=True) s13, s23, s43 = symbols('s13 s23 s43') x, y, z, w = symbols('x y z w') A = Matrix([[ 1, 0, s13, 0 ], [ 0, 1, s23, 0 ], [ 0, 0, 1, 0 ], [ 0, 0, s43, 1 ]]) print('\nLacroute\'s shear matrix (A.14) is:') pprint.pprint(A) # At this point, we can write "print(A*p)". However, we don't care # about the output "z" axis, as it is fixed. So we delete the 3rd row # of A. A.row_del(2) p = Matrix([ x, y, z, 1 ]) v = A*p print('\nAction of Lacroute\'s shear matrix on plane z (taking w=1):\n%s\n' % v) print('Scaling = %s' % (1/v[2])) print('Offset X = %s' % (v[0]/v[2]).subs(x, 0)) print('Offset Y = %s' % (v[1]/v[2]).subs(y, 0))