u (3bb2d6c0) · Commits · Michele Maris / LSPE_PRM

CHANGELOG

+1 −0

Original line number	Diff line number	Diff line
		2021 dec 01 - Repo created
		2021 dec 02 - Uploaded PyPRM 0.0
		2021 dec 03 - Added example
		2023 jul 12 - PyPRM 1.1.0 : Updated wobble angles operator in the original it was rotZ(omegaVA)rotY(zVAX) in the new is rotZ(omegaVA)rotY(zVAX)rotZ(-omegaVA)

+1 −0

+16 −3

Original line number	Diff line number	Diff line
		from .roto_translation import *

		__VERSION__="1.0.0 - 2021 Dec 03"
		__VERSION__="1.1.0 - 2021 Dec 03 - 2023 jul 12"
		__DESCRIPTION__="""
		M.Maris

		@@ -57,6 +57,7 @@ Issue
		0.0.2 - 2021 Nov 22 : Improved definition of matrices
		0.0.3 - 2021 Dec 2 : Structure improved, added example
		1.0.0 - 2021 Dec 3 : First Issue after review with the Pointing Working Group
		1.1.0 - 2023 Jul 12 : Changed the definition of wobble angles operator
		"""

		import numpy as np
		@@ -228,7 +229,14 @@ class PRM0 :
		self.ctheta0=0.
		self.commit()
		#
		def commit(self) :
		def commit(self,oldWobble=False) :
		""" if oldWobble==True the old version of wobble angle operator is applied

		current version : M3 = matr['omegaVAX'].dot(matr['zVAX'].dot(matr['omegaVAX'].T))
		old version : M3 = matr['omegaVAX'].dot(matr['zVAX'])

		the old version does not handle properly the cases zVAX=0.
		"""
		if self.isAllSet :
		out=self._matr['tilt'][0].dot(self._matr['pan'][0].dot(self._matr['roll'][0]))
		out=self._matr['ctheta0'][0].dot(out)
		@@ -239,7 +247,12 @@ class PRM0 :
		out.shape=(1,3,3)
		self._matr['M2']=out*1
		#
		if oldWobble :
		#old version
		out=self._matr['omegaVAX'][0].dot(self._matr['zVAX'][0])
		else :
		#new version
		out=self._matr['omegaVAX'][0].dot(self._matr['zVAX'][0].dot(self._matr['omegaVAX'][0].T))
		out.shape=(1,3,3)
		self._matr['M3']=out*1
		#