Miladiouss

def index_of_which_bin(bin_sizes, absolute_index, verbose=False):
    """
    Given the absolute index, returns which bin it falls in and which element of that bin it corresponds to.
    """
    # Which class/bin does i fall into?
    accum = np.add.accumulate(bin_sizes)
    if verbose:
        print("accum =", accum)
    bin_index  = len(np.argwhere(accum <= absolute_index))
    if verbose:

Miladiouss

# Author: Miladiouss (Milad Pourrahmani)
# Author's Website: Miladiouss.com
# Date: Sep 26, 2018
# Title: Simple Particle Simulator
# Description:
"""
I am hoping to assemble a series of brief educational computer programs that will capture the simplicity and elegance of laws of physics, both in nature and in apprehension.

This particular educational project (Simple Particle Simulator) is the first of this series. In less than 25 lines, plotting and comments aside, we can obtain the trajectory of a particle in ANY force field. The generated plot displays the simulated and exaggerated perihelion precision of Mercury. 

Miladiouss

# Library import for plotting
import matplotlib.pyplot as plt
%matplotlib inline

# Define Logistic Map
def f(x, r):
    return r * x * (1 - x)

# Define history and parameters
history = [0.5]

Miladiouss

# GPU Detection Tests
from torch.cuda import device_count
import gpustat

print("""
If pytorch raises "RuntimeError: cuda runtime error (30)" after suspension,
run the following commands in linux terminal:

sudo rmmod nvidia_uvm
sudo rmmod nvidia

Miladiouss

from datetime import datetime
timeFormat = '%a %b %d %Y at %H:%M:%S'
import socket
deviceName = socket.gethostname()
import getpass
userName = getpass.getuser()

t1 = datetime.now()
print(t1.strftime(timeFormat) + ' on ' + deviceName + ' by ' + userName)

Miladiouss

-- An HSC-PDR2 SQL query to extract RA, Dec, tract, patch of clean galaxies with a magnitude cut.

-- Useful Links:
-- /https://hsc-release.mtk.nao.ac.jp/hscMap-pdr2/app/#/?_=%7B%22view%22%3A%7B%22a%22%3A5.770914805732617,%22d%22%3A-0.019535944017082533,%22fovy%22%3A0.0002503311005562077,%22roll%22%3A0%7D,%22sspParams%22%3A%7B%22type%22%3A%22SDSS_TRUE_COLOR%22,%22filter%22%3A%5B%22HSC-I%22,%22HSC-R%22,%22HSC-G%22%5D,%22simpleRgb%22%3A%7B%22beta%22%3A22026.465794806718,%22a%22%3A1,%22bias%22%3A0.05,%22b0%22%3A0%7D,%22sdssTrueColor%22%3A%7B%22beta%22%3A26003.369421497522,%22a%22%3A1,%22bias%22%3A0.14439999999999995,%22b0%22%3A0%7D%7D%7D
-- https://hsc-release.mtk.nao.ac.jp/schema/#pdr2.pdr2_wide.forced


SELECT
        object_id 
        , tract

Miladiouss

def appendOrdinality(df, columns, ascending, strictlyIncreasing=False):
    """
    Appends absolute and relative ordinality to a dataframe.
    Oridinality is the normalized position of a row in a sorted dataframe.
    
    df:
        Pandas DataFrame
    columns: 
        List of column names for to be used for sorting (e.g. ['prob_1'])
    ascending: 

Miladiouss

# https://gist.github.com/Miladiouss/332b21e3c9d82ca2085b1396641cc80d
#-----------------------------------------------------------
# IPython
from IPython.core.interactiveshell import InteractiveShell
from IPython.display import clear_output
InteractiveShell.ast_node_interactivity = "all"
from IPython import display

#-----------------------------------------------------------
# Scientific

Miladiouss

from time import sleep, time
import numpy as np
from astropy import wcs

w = wcs.WCS(naxis=2)
w.wcs.ctype = ["RA---SIN","DEC--SIN"]
w.wcs.cdelt = (6.2/3600, 6.2/3600)
w.wcs.crval = (30, 60)
w.wcs.crota = (0, 15)
# ----------------------------------------    

Miladiouss

convert -resize 800x800 -delay 25 -loop 0 image-file-*.png animation.gif