






'''simple docstring'''


import random
import unittest

import numpy as np
import torch

from diffusers import DDIMScheduler, LDMSuperResolutionPipeline, UNetaDModel, VQModel
from diffusers.utils import PIL_INTERPOLATION, floats_tensor, load_image, slow, torch_device
from diffusers.utils.testing_utils import enable_full_determinism, require_torch


enable_full_determinism()




class 							lowerCAmelCase					(							unittest.TestCase   ):







   """simple docstring"""







   @property
   def    __magic_name__    ( self			)    ->   str:
    __a    :	Optional[Any]									=						1
    __a    :	Dict									=						3
    __a    :	Union[str, Any]									=						(32, 32)

    __a    :	Union[str, Any]									=						floats_tensor((batch_size, num_channels) + sizes   ,   rng=random.Random(0			)			).to(_A			)
    return image



   @property
   def    __magic_name__    ( self			)    ->   Optional[Any]:
    torch.manual_seed(0			)
    __a    :	List[str]									=						UNetaDModel(
        block_out_channels=(32, 64)   ,   layers_per_block=2   ,   sample_size=32   ,   in_channels=6   ,   out_channels=3   ,   down_block_types=('DownBlock2D', 'AttnDownBlock2D')   ,   up_block_types=('AttnUpBlock2D', 'UpBlock2D')   ,   )
    return model



   @property
   def    __magic_name__    ( self			)    ->   int:
    torch.manual_seed(0			)
    __a    :	List[str]									=						VQModel(
        block_out_channels=[32, 64]   ,   in_channels=3   ,   out_channels=3   ,   down_block_types=['DownEncoderBlock2D', 'DownEncoderBlock2D']   ,   up_block_types=['UpDecoderBlock2D', 'UpDecoderBlock2D']   ,   latent_channels=3   ,   )
    return model



   def    __magic_name__    ( self			)    ->   List[str]:
    __a    :	Tuple									=						'cpu'
    __a    :	int									=						self.dummy_uncond_unet
    __a    :	int									=						DDIMScheduler()
    __a    :	int									=						self.dummy_vq_model

    __a    :	Tuple									=						LDMSuperResolutionPipeline(unet=_A   ,   vqvae=_A   ,   scheduler=_A			)
    ldm.to(_A			)
    ldm.set_progress_bar_config(disable=_A			)

    __a    :	int									=						self.dummy_image.to(_A			)

    __a    :	List[Any]									=						torch.Generator(device=_A			).manual_seed(0			)
    __a    :	str									=						ldm(image=_A   ,   generator=_A   ,   num_inference_steps=2   ,   output_type='numpy'			).images

    __a    :	List[str]									=						image[0, -3:, -3:, -1]

    assert image.shape == (1, 64, 64, 3)
    __a    :	int									=						np.array([0.8_678, 0.8_245, 0.6_381, 0.6_830, 0.4_385, 0.5_599, 0.4_641, 0.6_201, 0.5_150]			)

    assert np.abs(image_slice.flatten() - expected_slice			).max() < 1E-2



   @unittest.skipIf(torch_device != 'cuda'   ,   'This test requires a GPU'			)
   def    __magic_name__    ( self			)    ->   Any:
    __a    :	List[str]									=						self.dummy_uncond_unet
    __a    :	Tuple									=						DDIMScheduler()
    __a    :	int									=						self.dummy_vq_model

    # put models in fp16
    __a    :	str									=						unet.half()
    __a    :	Optional[Any]									=						vqvae.half()

    __a    :	str									=						LDMSuperResolutionPipeline(unet=_A   ,   vqvae=_A   ,   scheduler=_A			)
    ldm.to(_A			)
    ldm.set_progress_bar_config(disable=_A			)

    __a    :	Optional[Any]									=						self.dummy_image.to(_A			)

    __a    :	str									=						ldm(_A   ,   num_inference_steps=2   ,   output_type='numpy'			).images

    assert image.shape == (1, 64, 64, 3)




@slow
@require_torch
class 							lowerCAmelCase					(							unittest.TestCase   ):







   """simple docstring"""







   def    __magic_name__    ( self			)    ->   List[Any]:
    __a    :	int									=						load_image(
        'https://huggingface.co/datasets/hf-internal-testing/diffusers-images/resolve/main'
        '/vq_diffusion/teddy_bear_pool.png'			)
    __a    :	Tuple									=						init_image.resize((64, 64)   ,   resample=PIL_INTERPOLATION['lanczos']			)

    __a    :	Any									=						LDMSuperResolutionPipeline.from_pretrained('duongna/ldm-super-resolution'   ,   device_map='auto'			)
    ldm.set_progress_bar_config(disable=_A			)

    __a    :	Tuple									=						torch.manual_seed(0			)
    __a    :	Tuple									=						ldm(image=_A   ,   generator=_A   ,   num_inference_steps=20   ,   output_type='numpy'			).images

    __a    :	Optional[int]									=						image[0, -3:, -3:, -1]

    assert image.shape == (1, 256, 256, 3)
    __a    :	Optional[int]									=						np.array([0.7_644, 0.7_679, 0.7_642, 0.7_633, 0.7_666, 0.7_560, 0.7_425, 0.7_257, 0.6_907]			)

    assert np.abs(image_slice.flatten() - expected_slice			).max() < 1E-2












'''simple docstring'''


from __future__ import annotations

import random

from matplotlib import animation
from matplotlib import pyplot as plt

# Frame rate of the animation
SCREAMING_SNAKE_CASE_				      =							2_0

# Time between time steps in seconds
SCREAMING_SNAKE_CASE_				      =							INTERVAL / 1_0_0_0




class 							lowerCAmelCase					:







   """simple docstring"""







   def __init__( self   ,   _A   ,   _A   ,   _A   ,   _A   ,   _A = 1.0   ,   _A = 1.0   ,   _A = "blue"   ,   )    ->   None:
    __a    :	str									=						position_x
    __a    :	Optional[Any]									=						position_y
    __a    :	List[Any]									=						velocity_x
    __a    :	Tuple									=						velocity_y
    __a    :	Dict									=						mass
    __a    :	int									=						size
    __a    :	Tuple									=						color



   @property
   def    __magic_name__    ( self			)    ->   tuple[float, float]:
    return self.position_x, self.position_y



   @property
   def    __magic_name__    ( self			)    ->   tuple[float, float]:
    return self.velocity_x, self.velocity_y



   def    __magic_name__    ( self   ,   _A   ,   _A   ,   _A			)    ->   None:
    self.velocity_x += force_x * delta_time
    self.velocity_y += force_y * delta_time



   def    __magic_name__    ( self   ,   _A			)    ->   None:
    self.position_x += self.velocity_x * delta_time
    self.position_y += self.velocity_y * delta_time




class 							lowerCAmelCase					:







   """simple docstring"""







   def __init__( self   ,   _A   ,   _A = 1.0   ,   _A = 1.0   ,   _A = 0.0   ,   )    ->   None:
    __a    :	List[str]									=						bodies
    __a    :	List[str]									=						gravitation_constant
    __a    :	int									=						time_factor
    __a    :	List[str]									=						softening_factor



   def __len__( self			)    ->   int:
    return len(self.bodies			)



   def    __magic_name__    ( self   ,   _A			)    ->   None:
    for bodya in self.bodies:
     __a    :	Any									=						0.0
     __a    :	str									=						0.0
     for bodya in self.bodies:
      if bodya != bodya:
       __a    :	Optional[Any]									=						bodya.position_x - bodya.position_x
       __a    :	Dict									=						bodya.position_y - bodya.position_y

       # Calculation of the distance using Pythagoras's theorem
       # Extra factor due to the softening technique
       __a    :	Optional[Any]									=						(dif_x**2 + dif_y**2 + self.softening_factor) ** (
           1 / 2
       )

       # Newton's law of universal gravitation.
       force_x += (
           self.gravitation_constant * bodya.mass * dif_x / distance**3
       )
       force_y += (
           self.gravitation_constant * bodya.mass * dif_y / distance**3
       )

            # Update the body's velocity once all the force components have been added
     bodya.update_velocity(_A   ,   _A   ,   delta_time * self.time_factor			)

    # Update the positions only after all the velocities have been updated
    for body in self.bodies:
     body.update_position(delta_time * self.time_factor			)






def   lowerCAmelCase__			(					SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__		):
 body_system.update_system(SCREAMING_SNAKE_CASE__		)

 # Update the positions of the patches
 for patch, body in zip(SCREAMING_SNAKE_CASE__       ,		body_system.bodies		):
  __a    :	Dict									=						(body.position_x, body.position_y)






def   lowerCAmelCase__			(					SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__ = -1       ,		SCREAMING_SNAKE_CASE__ = 1       ,		SCREAMING_SNAKE_CASE__ = -1       ,		SCREAMING_SNAKE_CASE__ = 1       ,		):
 __a    :	int									=						plt.figure()
 fig.canvas.set_window_title(SCREAMING_SNAKE_CASE__		)
 __a    :	int									=						plt.axes(
     xlim=(x_start, x_end)       ,		ylim=(y_start, y_end)		)  # Set section to be plotted
 plt.gca().set_aspect('equal'		)  # Fix aspect ratio

 # Each body is drawn as a patch by the plt-function
 __a    :	str									=						[
     plt.Circle((body.position_x, body.position_y)       ,		body.size       ,		fc=body.color		)
     for body in body_system.bodies
 ]

 for patch in patches:
  ax.add_patch(SCREAMING_SNAKE_CASE__		)

 # Function called at each step of the animation
 def update(SCREAMING_SNAKE_CASE__		) -> list[plt.Circle]:
  update_step(SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__		)
  return patches

 __a    :	Dict									=						animation.FuncAnimation(  # noqa: F841
     SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		interval=SCREAMING_SNAKE_CASE__       ,		blit=SCREAMING_SNAKE_CASE__		)

 plt.show()






def   lowerCAmelCase__			(					):
 __a    :	Dict									=						0.970_0436
 __a    :	List[Any]									=						-0.2430_8753
 __a    :	Union[str, Any]									=						0.4_6620_3685
 __a    :	Any									=						0.4323_6573

 __a    :	Optional[int]									=						[
     Body(SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		size=0.2       ,		color='red'		),
     Body(-position_x       ,		-position_y       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		size=0.2       ,		color='green'		),
     Body(0       ,		0       ,		-2 * velocity_x       ,		-2 * velocity_y       ,		size=0.2       ,		color='blue'		),
 ]
 return BodySystem(SCREAMING_SNAKE_CASE__       ,		time_factor=3		)






def   lowerCAmelCase__			(					):
 __a    :	Optional[int]									=						7.34_76E22
 __a    :	List[Any]									=						5.9_72E24
 __a    :	List[Any]									=						1022
 __a    :	str									=						384399000
 __a    :	Optional[int]									=						6.6_74E-11

 # Calculation of the respective velocities so that total impulse is zero,
 # i.e. the two bodies together don't move
 __a    :	str									=						earth_mass * velocity_dif / (earth_mass + moon_mass)
 __a    :	Optional[int]									=						moon_velocity - velocity_dif

 __a    :	Union[str, Any]									=						Body(-earth_moon_distance       ,		0       ,		0       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		10000000       ,		'grey'		)
 __a    :	Union[str, Any]									=						Body(0       ,		0       ,		0       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		50000000       ,		'blue'		)
 return BodySystem([earth, moon]       ,		SCREAMING_SNAKE_CASE__       ,		time_factor=1000000		)






def   lowerCAmelCase__			(					):
 __a    :	Tuple									=						[]
 for _ in range(10		):
  __a    :	Dict									=						random.uniform(-0.5       ,		0.5		)
  __a    :	Optional[int]									=						random.uniform(-0.5       ,		0.5		)

  # Bodies are created pairwise with opposite velocities so that the
  # total impulse remains zero
  bodies.append(
      Body(
          random.uniform(-0.5       ,		0.5		)       ,		random.uniform(-0.5       ,		0.5		)       ,		SCREAMING_SNAKE_CASE__       ,		SCREAMING_SNAKE_CASE__       ,		size=0.05       ,		)		)
  bodies.append(
      Body(
          random.uniform(-0.5       ,		0.5		)       ,		random.uniform(-0.5       ,		0.5		)       ,		-velocity_x       ,		-velocity_y       ,		size=0.05       ,		)		)
 return BodySystem(SCREAMING_SNAKE_CASE__       ,		0.01       ,		10       ,		0.1		)


if __name__ == "__main__":
  plot("Figure-8 solution to the 3-body-problem", example_a(), -2, 2, -2, 2)
  plot(
      "Moon's orbit around the earth",
      example_a(),
      -4_3_0_0_0_0_0_0_0,
      4_3_0_0_0_0_0_0_0,
      -4_3_0_0_0_0_0_0_0,
      4_3_0_0_0_0_0_0_0,
  )
  plot("Random system with many bodies", example_a(), -1.5, 1.5, -1.5, 1.5)




