def get_matrix_local(self):
		if self._tag == 'OB':
			return self.bdata.matrix_local.copy()
		elif self._tag == 'DP':
			return self._ref.matrix_world.inverted_safe() * self._dupli_matrix
		else:  # 'BO', current pose
			# PoseBone.matrix is in armature space, bring in back in real local one!
			par = self.bdata.parent	
			c_bone = self._ref.pose.bones[self.bdata.name]
			par_mat_inv = self._ref.pose.bones[par.name].matrix.inverted_safe() if par else Matrix()			
		
			# There's a bug with bones matrix, leading to wrong results when multiplying matrices
			# at some specific rotation values and when scale is different from 1.0
			# Extracting translation, rotation and scale matrices then recompose them seems to fix it		
			t_mat = Matrix.Translation(c_bone.matrix.to_translation())
			r_mat = c_bone.matrix.to_quaternion().to_matrix().to_4x4()			
			s_mat = extract_scale_matrix(c_bone.matrix.to_3x3()).to_4x4()
			bone_mat = t_mat * r_mat * s_mat			
			
			return par_mat_inv * bone_mat

def extract_scale_matrix(matrix):
	s_x = Vector((matrix[0][0], matrix[1][0],matrix[2][0])).magnitude
	s_y = Vector((matrix[0][1], matrix[1][1],matrix[2][1])).magnitude
	s_z = Vector((matrix[0][2], matrix[1][2],matrix[2][2])).magnitude		
	s_mat = Matrix(((s_x, 0, 0), (0, s_y, 0), (0, 0, s_z)))
	
	return s_mat