code: def test_to_joint_gaussian(self):
      jgd1 = self.phi1.to_joint_gaussian()
      jgd2 = self.phi2.to_joint_gaussian()
  
      self.assertEqual(jgd1.variables, self.gauss_phi1.variables)
      np_test.assert_almost_equal(
          jgd1.covariance, self.gauss_phi1.covariance
      )
      np_test.assert_almost_equal(jgd1.mean, self.gauss_phi1.mean)
  
      self.assertEqual(jgd2.variables, self.gauss_phi2.variables)
      np_test.assert_almost_equal(
          jgd2.covariance, self.gauss_phi2.covariance
      )
      np_test.assert_almost_equal(jgd2.mean, self.gauss_phi2.mean)
  
file_path: pgmpy/pgmpy/tests/test_factors/test_continuous/test_Canonical_Factor.py
imports: import unittest
  import numpy
  import numpy.testing
  from pgmpy.factors.distributions import GaussianDistribution
  from pgmpy.factors.continuous import CanonicalDistribution
label_desc: Use image node_id 3 for calling the TestJGDMethods obj's underlying member method code with example usage: obj.test_to_joint_gaussian() without return types