Skip to content
Snippets Groups Projects
Metavoxel GPU Validation 2.ipynb 211 KiB
Newer Older
Amira Abdel-Rahman's avatar
Amira Abdel-Rahman committed
3001 3002 3003 3004 3005 3006 3007 3008 3009 3010 3011 3012 3013 3014 3015 3016 3017 3018 3019 3020 3021 3022 3023 3024 3025 3026 3027 3028 3029 3030 3031 3032 3033 3034 3035 3036 3037 3038 3039 3040 3041 3042 3043 3044 3045 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056 3057 3058 3059 3060 3061 3062 3063 3064 3065 3066 3067 3068 3069 3070 3071 3072 3073 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083 3084 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098 3099 3100 3101 3102 3103 3104 3105 3106 3107 3108 3109 3110 3111 3112 3113 3114 3115 3116 3117 3118 3119 3120 3121 3122 3123 3124 3125 3126 3127 3128 3129 3130 3131 3132 3133 3134 3135 3136 3137 3138 3139 3140 3141 3142 3143 3144 3145 3146 3147 3148 3149 3150 3151 3152 3153 3154 3155 3156 3157 3158 3159 3160 3161 3162 3163 3164 3165 3166 3167 3168 3169 3170 3171 3172 3173 3174 3175 3176 3177 3178 3179 3180 3181 3182 3183 3184 3185 3186 3187 3188 3189 3190 3191 3192 3193 3194 3195 3196 3197 3198 3199 3200 3201 3202 3203 3204 3205 3206 3207 3208 3209 3210 3211 3212 3213 3214 3215 3216 3217 3218 3219 3220 3221 3222 3223 3224 3225 3226 3227 3228 3229 3230 3231 3232 3233 3234 3235 3236 3237 3238 3239 3240 3241 3242 3243 3244 3245 3246 3247 3248 3249 3250 3251 3252 3253 3254 3255 3256 3257 3258 3259 3260 3261 3262 3263 3264 3265 3266 3267 3268 3269 3270 3271 3272 3273 3274 3275 3276 3277 3278 3279 3280 3281 3282 3283 3284 3285 3286 3287 3288 3289 3290 3291 3292 3293 3294 3295 3296 3297 3298 3299 3300 3301 3302 3303 3304 3305 3306 3307 3308 3309 3310 3311 3312 3313 3314 3315 3316 3317 3318 3319 3320 3321 3322 3323 3324 3325 3326 3327 3328 3329 3330 3331 3332 3333 3334 3335 3336 3337 3338 3339 3340
       "  1661.17,980.458 1661.82,980.458 1662.46,980.458 1663.11,980.458 1663.76,980.458 1664.4,980.458 1665.05,980.458 1665.69,980.458 1666.34,980.458 1666.99,980.458 \n",
       "  1667.63,980.458 1668.28,980.458 1668.93,980.458 1669.57,980.458 1670.22,980.458 1670.87,980.458 1671.51,980.458 1672.16,980.458 1672.8,980.458 1673.45,980.458 \n",
       "  1674.1,980.458 1674.74,980.458 1675.39,980.458 1676.04,980.458 1676.68,980.458 1677.33,980.458 1677.98,980.458 1678.62,980.458 1679.27,980.458 1679.91,980.458 \n",
       "  1680.56,980.458 1681.21,980.458 1681.85,980.458 1682.5,980.458 1683.15,980.458 1683.79,980.458 1684.44,980.458 1685.08,980.458 1685.73,980.458 1686.38,980.458 \n",
       "  1687.02,980.458 1687.67,980.458 1688.32,980.458 1688.96,980.458 1689.61,980.458 1690.26,980.458 1690.9,980.458 1691.55,980.458 1692.19,980.458 1692.84,980.458 \n",
       "  1693.49,980.458 1694.13,980.458 1694.78,980.458 1695.43,980.458 1696.07,980.458 1696.72,980.458 1697.37,980.458 1698.01,980.458 1698.66,980.458 1699.3,980.458 \n",
       "  1699.95,980.458 1700.6,980.458 1701.24,980.458 1701.89,980.458 1702.54,980.458 1703.18,980.458 1703.83,980.458 1704.48,980.458 1705.12,980.458 1705.77,980.458 \n",
       "  1706.41,980.458 1707.06,980.458 1707.71,980.458 1708.35,980.458 1709,980.458 1709.65,980.458 1710.29,980.458 1710.94,980.458 1711.58,980.458 1712.23,980.458 \n",
       "  1712.88,980.458 1713.52,980.458 1714.17,980.458 1714.82,980.458 1715.46,980.458 1716.11,980.458 1716.76,980.458 1717.4,980.458 1718.05,980.458 1718.69,980.458 \n",
       "  1719.34,980.458 1719.99,980.458 1720.63,980.458 1721.28,980.458 1721.93,980.458 1722.57,980.458 1723.22,980.458 1723.87,980.458 1724.51,980.458 1725.16,980.458 \n",
       "  1725.8,980.458 1726.45,980.458 1727.1,980.458 1727.74,980.458 1728.39,980.458 1729.04,980.458 1729.68,980.458 1730.33,980.458 1730.98,980.458 1731.62,980.458 \n",
       "  1732.27,980.458 1732.91,980.458 1733.56,980.458 1734.21,980.458 1734.85,980.458 1735.5,980.458 1736.15,980.458 1736.79,980.458 1737.44,980.458 1738.08,980.458 \n",
       "  1738.73,980.458 1739.38,980.458 1740.02,980.458 1740.67,980.458 1741.32,980.458 1741.96,980.458 1742.61,980.458 1743.26,980.458 1743.9,980.458 1744.55,980.458 \n",
       "  1745.19,980.458 1745.84,980.458 1746.49,980.458 1747.13,980.458 1747.78,980.458 1748.43,980.458 1749.07,980.458 1749.72,980.458 1750.37,980.458 1751.01,980.458 \n",
       "  1751.66,980.458 1752.3,980.458 1752.95,980.458 1753.6,980.458 1754.24,980.458 1754.89,980.458 1755.54,980.458 1756.18,980.458 1756.83,980.458 1757.48,980.458 \n",
       "  1758.12,980.458 1758.77,980.458 1759.41,980.458 1760.06,980.458 1760.71,980.458 1761.35,980.458 1762,980.458 1762.65,980.458 1763.29,980.458 1763.94,980.458 \n",
       "  1764.58,980.458 1765.23,980.458 1765.88,980.458 1766.52,980.458 1767.17,980.458 1767.82,980.458 1768.46,980.458 1769.11,980.458 1769.76,980.458 1770.4,980.458 \n",
       "  1771.05,980.458 1771.69,980.458 1772.34,980.458 1772.99,980.458 1773.63,980.458 1774.28,980.458 1774.93,980.458 1775.57,980.458 1776.22,980.458 1776.87,980.458 \n",
       "  1777.51,980.458 1778.16,980.458 1778.8,980.458 1779.45,980.458 1780.1,980.458 1780.74,980.458 1781.39,980.458 1782.04,980.458 1782.68,980.458 1783.33,980.458 \n",
       "  1783.98,980.458 1784.62,980.458 1785.27,980.458 1785.91,980.458 1786.56,980.458 1787.21,980.458 1787.85,980.458 1788.5,980.458 1789.15,980.458 1789.79,980.458 \n",
       "  1790.44,980.458 1791.08,980.458 1791.73,980.458 1792.38,980.458 1793.02,980.458 1793.67,980.458 1794.32,980.458 1794.96,980.458 1795.61,980.458 1796.26,980.458 \n",
       "  1796.9,980.458 1797.55,980.458 1798.19,980.458 1798.84,980.458 1799.49,980.458 1800.13,980.458 1800.78,980.458 1801.43,980.458 1802.07,980.458 1802.72,980.458 \n",
       "  1803.37,980.458 1804.01,980.458 1804.66,980.458 1805.3,980.458 1805.95,980.458 1806.6,980.458 1807.24,980.458 1807.89,980.458 1808.54,980.458 1809.18,980.458 \n",
       "  1809.83,980.458 1810.48,980.458 1811.12,980.458 1811.77,980.458 1812.41,980.458 1813.06,980.458 1813.71,980.458 1814.35,980.458 1815,980.458 1815.65,980.458 \n",
       "  1816.29,980.458 1816.94,980.458 1817.58,980.458 1818.23,980.458 1818.88,980.458 1819.52,980.458 1820.17,980.458 1820.82,980.458 1821.46,980.458 1822.11,980.458 \n",
       "  1822.76,980.458 1823.4,980.458 1824.05,980.458 1824.69,980.458 1825.34,980.458 1825.99,980.458 1826.63,980.458 1827.28,980.458 1827.93,980.458 1828.57,980.458 \n",
       "  1829.22,980.458 1829.87,980.458 1830.51,980.458 1831.16,980.458 1831.8,980.458 1832.45,980.458 1833.1,980.458 1833.74,980.458 1834.39,980.458 1835.04,980.458 \n",
       "  1835.68,980.458 1836.33,980.458 1836.98,980.458 1837.62,980.458 1838.27,980.458 1838.91,980.458 1839.56,980.458 1840.21,980.458 1840.85,980.458 1841.5,980.458 \n",
       "  1842.15,980.458 1842.79,980.458 1843.44,980.458 1844.08,980.458 1844.73,980.458 1845.38,980.458 1846.02,980.458 1846.67,980.458 1847.32,980.458 1847.96,980.458 \n",
       "  1848.61,980.458 1849.26,980.458 1849.9,980.458 1850.55,980.458 1851.19,980.458 1851.84,980.458 1852.49,980.458 1853.13,980.458 1853.78,980.458 1854.43,980.458 \n",
       "  1855.07,980.458 1855.72,980.458 1856.37,980.458 1857.01,980.458 1857.66,980.458 1858.3,980.458 1858.95,980.458 1859.6,980.458 1860.24,980.458 1860.89,980.458 \n",
       "  1861.54,980.458 1862.18,980.458 1862.83,980.458 1863.48,980.458 1864.12,980.458 1864.77,980.458 1865.41,980.458 1866.06,980.458 1866.71,980.458 1867.35,980.458 \n",
       "  1868,980.458 1868.65,980.458 1869.29,980.458 1869.94,980.458 1870.58,980.458 1871.23,980.458 1871.88,980.458 1872.52,980.458 1873.17,980.458 1873.82,980.458 \n",
       "  1874.46,980.458 1875.11,980.458 1875.76,980.458 1876.4,980.458 1877.05,980.458 1877.69,980.458 1878.34,980.458 1878.99,980.458 1879.63,980.458 1880.28,980.458 \n",
       "  1880.93,980.458 1881.57,980.458 1882.22,980.458 1882.87,980.458 1883.51,980.458 1884.16,980.458 1884.8,980.458 1885.45,980.458 1886.1,980.458 1886.74,980.458 \n",
       "  1887.39,980.458 1888.04,980.458 1888.68,980.458 1889.33,980.458 1889.98,980.458 1890.62,980.458 1891.27,980.458 1891.91,980.458 1892.56,980.458 1893.21,980.458 \n",
       "  1893.85,980.458 1894.5,980.458 1895.15,980.458 1895.79,980.458 1896.44,980.458 1897.08,980.458 1897.73,980.458 1898.38,980.458 1899.02,980.458 1899.67,980.458 \n",
       "  1900.32,980.458 1900.96,980.458 1901.61,980.458 1902.26,980.458 1902.9,980.458 1903.55,980.458 1904.19,980.458 1904.84,980.458 1905.49,980.458 1906.13,980.458 \n",
       "  1906.78,980.458 1907.43,980.458 1908.07,980.458 1908.72,980.458 1909.37,980.458 1910.01,980.458 1910.66,980.458 1911.3,980.458 1911.95,980.458 1912.6,980.458 \n",
       "  1913.24,980.458 1913.89,980.458 1914.54,980.458 1915.18,980.458 1915.83,980.458 1916.48,980.458 1917.12,980.458 1917.77,980.458 1918.41,980.458 1919.06,980.458 \n",
       "  1919.71,980.458 1920.35,980.458 1921,980.458 1921.65,980.458 1922.29,980.458 1922.94,980.458 1923.58,980.458 1924.23,980.458 1924.88,980.458 1925.52,980.458 \n",
       "  1926.17,980.458 1926.82,980.458 1927.46,980.458 1928.11,980.458 1928.76,980.458 1929.4,980.458 1930.05,980.458 1930.69,980.458 1931.34,980.458 1931.99,980.458 \n",
       "  1932.63,980.458 1933.28,980.458 1933.93,980.458 1934.57,980.458 1935.22,980.458 1935.87,980.458 1936.51,980.458 1937.16,980.458 1937.8,980.458 1938.45,980.458 \n",
       "  1939.1,980.458 1939.74,980.458 1940.39,980.458 1941.04,980.458 1941.68,980.458 1942.33,980.458 1942.98,980.458 1943.62,980.458 1944.27,980.458 1944.91,980.458 \n",
       "  1945.56,980.458 1946.21,980.458 1946.85,980.458 1947.5,980.458 1948.15,980.458 1948.79,980.458 1949.44,980.458 1950.08,980.458 1950.73,980.458 1951.38,980.458 \n",
       "  1952.02,980.458 1952.67,980.458 1953.32,980.458 1953.96,980.458 1954.61,980.458 1955.26,980.458 1955.9,980.458 1956.55,980.458 1957.19,980.458 1957.84,980.458 \n",
       "  1958.49,980.458 1959.13,980.458 1959.78,980.458 1960.43,980.458 1961.07,980.458 1961.72,980.458 1962.37,980.458 1963.01,980.458 1963.66,980.458 1964.3,980.458 \n",
       "  1964.95,980.458 1965.6,980.458 1966.24,980.458 1966.89,980.458 1967.54,980.458 1968.18,980.458 1968.83,980.458 1969.48,980.458 1970.12,980.458 1970.77,980.458 \n",
       "  1971.41,980.458 1972.06,980.458 1972.71,980.458 1973.35,980.458 1974,980.458 1974.65,980.458 1975.29,980.458 1975.94,980.458 1976.58,980.458 1977.23,980.458 \n",
       "  1977.88,980.458 1978.52,980.458 1979.17,980.458 1979.82,980.458 1980.46,980.458 1981.11,980.458 1981.76,980.458 1982.4,980.458 1983.05,980.458 1983.69,980.458 \n",
       "  1984.34,980.458 1984.99,980.458 1985.63,980.458 1986.28,980.458 1986.93,980.458 1987.57,980.458 1988.22,980.458 1988.87,980.458 1989.51,980.458 1990.16,980.458 \n",
       "  1990.8,980.458 1991.45,980.458 1992.1,980.458 1992.74,980.458 1993.39,980.458 1994.04,980.458 1994.68,980.458 1995.33,980.458 1995.98,980.458 1996.62,980.458 \n",
       "  1997.27,980.458 1997.91,980.458 1998.56,980.458 1999.21,980.458 1999.85,980.458 2000.5,980.458 2001.15,980.458 2001.79,980.458 2002.44,980.458 2003.08,980.458 \n",
       "  2003.73,980.458 2004.38,980.458 2005.02,980.458 2005.67,980.458 2006.32,980.458 2006.96,980.458 2007.61,980.458 2008.26,980.458 2008.9,980.458 2009.55,980.458 \n",
       "  2010.19,980.458 2010.84,980.458 2011.49,980.458 2012.13,980.458 2012.78,980.458 2013.43,980.458 2014.07,980.458 2014.72,980.458 2015.37,980.458 2016.01,980.458 \n",
       "  2016.66,980.458 2017.3,980.458 2017.95,980.458 2018.6,980.458 2019.24,980.458 2019.89,980.458 2020.54,980.458 2021.18,980.458 2021.83,980.458 2022.48,980.458 \n",
       "  2023.12,980.458 2023.77,980.458 2024.41,980.458 2025.06,980.458 2025.71,980.458 2026.35,980.458 2027,980.458 2027.65,980.458 2028.29,980.458 2028.94,980.458 \n",
       "  2029.58,980.458 2030.23,980.458 2030.88,980.458 2031.52,980.458 2032.17,980.458 2032.82,980.458 2033.46,980.458 2034.11,980.458 2034.76,980.458 2035.4,980.458 \n",
       "  2036.05,980.458 2036.69,980.458 2037.34,980.458 2037.99,980.458 2038.63,980.458 2039.28,980.458 2039.93,980.458 2040.57,980.458 2041.22,980.458 2041.87,980.458 \n",
       "  2042.51,980.458 2043.16,980.458 2043.8,980.458 2044.45,980.458 2045.1,980.458 2045.74,980.458 2046.39,980.458 2047.04,980.458 2047.68,980.458 2048.33,980.458 \n",
       "  2048.98,980.458 2049.62,980.458 2050.27,980.458 2050.91,980.458 2051.56,980.458 2052.21,980.458 2052.85,980.458 2053.5,980.458 2054.15,980.458 2054.79,980.458 \n",
       "  2055.44,980.458 2056.09,980.458 2056.73,980.458 2057.38,980.458 2058.02,980.458 2058.67,980.458 2059.32,980.458 2059.96,980.458 2060.61,980.458 2061.26,980.458 \n",
       "  2061.9,980.458 2062.55,980.458 2063.19,980.458 2063.84,980.458 2064.49,980.458 2065.13,980.458 2065.78,980.458 2066.43,980.458 2067.07,980.458 2067.72,980.458 \n",
       "  2068.37,980.458 2069.01,980.458 2069.66,980.458 2070.3,980.458 2070.95,980.458 2071.6,980.458 2072.24,980.458 2072.89,980.458 2073.54,980.458 2074.18,980.458 \n",
       "  2074.83,980.458 2075.48,980.458 2076.12,980.458 2076.77,980.458 2077.41,980.458 2078.06,980.458 2078.71,980.458 2079.35,980.458 2080,980.458 2080.65,980.458 \n",
       "  2081.29,980.458 2081.94,980.458 2082.59,980.458 2083.23,980.458 2083.88,980.458 2084.52,980.458 2085.17,980.458 2085.82,980.458 2086.46,980.458 2087.11,980.458 \n",
       "  2087.76,980.458 2088.4,980.458 2089.05,980.458 2089.69,980.458 2090.34,980.458 2090.99,980.458 2091.63,980.458 2092.28,980.458 2092.93,980.458 2093.57,980.458 \n",
       "  2094.22,980.458 2094.87,980.458 2095.51,980.458 2096.16,980.458 2096.8,980.458 2097.45,980.458 2098.1,980.458 2098.74,980.458 2099.39,980.458 2100.04,980.458 \n",
       "  2100.68,980.458 2101.33,980.458 2101.98,980.458 2102.62,980.458 2103.27,980.458 2103.91,980.458 2104.56,980.458 2105.21,980.458 2105.85,980.458 2106.5,980.458 \n",
       "  2107.15,980.458 2107.79,980.458 2108.44,980.458 2109.09,980.458 2109.73,980.458 2110.38,980.458 2111.02,980.458 2111.67,980.458 2112.32,980.458 2112.96,980.458 \n",
       "  2113.61,980.458 2114.26,980.458 2114.9,980.458 2115.55,980.458 2116.19,980.458 2116.84,980.458 2117.49,980.458 2118.13,980.458 2118.78,980.458 2119.43,980.458 \n",
       "  2120.07,980.458 2120.72,980.458 2121.37,980.458 2122.01,980.458 2122.66,980.458 2123.3,980.458 2123.95,980.458 2124.6,980.458 2125.24,980.458 2125.89,980.458 \n",
       "  2126.54,980.458 2127.18,980.458 2127.83,980.458 2128.48,980.458 2129.12,980.458 2129.77,980.458 2130.41,980.458 2131.06,980.458 2131.71,980.458 2132.35,980.458 \n",
       "  2133,980.458 2133.65,980.458 2134.29,980.458 2134.94,980.458 2135.59,980.458 2136.23,980.458 2136.88,980.458 2137.52,980.458 2138.17,980.458 2138.82,980.458 \n",
       "  2139.46,980.458 2140.11,980.458 2140.76,980.458 2141.4,980.458 2142.05,980.458 2142.69,980.458 2143.34,980.458 2143.99,980.458 2144.63,980.458 2145.28,980.458 \n",
       "  2145.93,980.458 2146.57,980.458 2147.22,980.458 2147.87,980.458 2148.51,980.458 2149.16,980.458 2149.8,980.458 2150.45,980.458 2151.1,980.458 2151.74,980.458 \n",
       "  2152.39,980.458 2153.04,980.458 2153.68,980.458 2154.33,980.458 2154.98,980.458 2155.62,980.458 2156.27,980.458 2156.91,980.458 2157.56,980.458 2158.21,980.458 \n",
       "  2158.85,980.458 2159.5,980.458 2160.15,980.458 2160.79,980.458 2161.44,980.458 2162.09,980.458 2162.73,980.458 2163.38,980.458 2164.02,980.458 2164.67,980.458 \n",
       "  2165.32,980.458 2165.96,980.458 2166.61,980.458 2167.26,980.458 2167.9,980.458 2168.55,980.458 2169.19,980.458 2169.84,980.458 2170.49,980.458 2171.13,980.458 \n",
       "  2171.78,980.458 2172.43,980.458 2173.07,980.458 2173.72,980.458 2174.37,980.458 2175.01,980.458 2175.66,980.458 2176.3,980.458 2176.95,980.458 2177.6,980.458 \n",
       "  2178.24,980.458 2178.89,980.458 2179.54,980.458 2180.18,980.458 2180.83,980.458 2181.48,980.458 2182.12,980.458 2182.77,980.458 2183.41,980.458 2184.06,980.458 \n",
       "  2184.71,980.458 2185.35,980.458 2186,980.458 2186.65,980.458 2187.29,980.458 2187.94,980.458 2188.59,980.458 2189.23,980.458 2189.88,980.458 2190.52,980.458 \n",
       "  2191.17,980.458 2191.82,980.458 2192.46,980.458 2193.11,980.458 2193.76,980.458 2194.4,980.458 2195.05,980.458 2195.69,980.458 2196.34,980.458 2196.99,980.458 \n",
       "  2197.63,980.458 2198.28,980.458 2198.93,980.458 2199.57,980.458 2200.22,980.458 2200.87,980.458 2201.51,980.458 2202.16,980.458 2202.8,980.458 2203.45,980.458 \n",
       "  2204.1,980.458 2204.74,980.458 2205.39,980.458 2206.04,980.458 2206.68,980.458 2207.33,980.458 2207.98,980.458 2208.62,980.458 2209.27,980.458 2209.91,980.458 \n",
       "  2210.56,980.458 2211.21,980.458 2211.85,980.458 2212.5,980.458 2213.15,980.458 2213.79,980.458 2214.44,980.458 2215.09,980.458 2215.73,980.458 2216.38,980.458 \n",
       "  2217.02,980.458 2217.67,980.458 2218.32,980.458 2218.96,980.458 2219.61,980.458 2220.26,980.458 2220.9,980.458 2221.55,980.458 2222.19,980.458 2222.84,980.458 \n",
       "  2223.49,980.458 2224.13,980.458 2224.78,980.458 2225.43,980.458 2226.07,980.458 2226.72,980.458 2227.37,980.458 2228.01,980.458 2228.66,980.458 2229.3,980.458 \n",
       "  2229.95,980.458 2230.6,980.458 2231.24,980.458 2231.89,980.458 2232.54,980.458 2233.18,980.458 2233.83,980.458 2234.48,980.458 2235.12,980.458 2235.77,980.458 \n",
       "  2236.41,980.458 2237.06,980.458 2237.71,980.458 2238.35,980.458 2239,980.458 2239.65,980.458 2240.29,980.458 2240.94,980.458 2241.59,980.458 2242.23,980.458 \n",
       "  2242.88,980.458 2243.52,980.458 2244.17,980.458 2244.82,980.458 2245.46,980.458 2246.11,980.458 2246.76,980.458 2247.4,980.458 2248.05,980.458 2248.69,980.458 \n",
       "  2249.34,980.458 2249.99,980.458 2250.63,980.458 2251.28,980.458 2251.93,980.458 2252.57,980.458 2253.22,980.458 2253.87,980.458 2254.51,980.458 2255.16,980.458 \n",
       "  2255.8,980.458 2256.45,980.458 2257.1,980.458 2257.74,980.458 2258.39,980.458 2259.04,980.458 2259.68,980.458 2260.33,980.458 2260.98,980.458 2261.62,980.458 \n",
       "  2262.27,980.458 2262.91,980.458 2263.56,980.458 2264.21,980.458 2264.85,980.458 2265.5,980.458 2266.15,980.458 2266.79,980.458 2267.44,980.458 2268.09,980.458 \n",
       "  2268.73,980.458 2269.38,980.458 2270.02,980.458 2270.67,980.458 2271.32,980.458 2271.96,980.458 2272.61,980.458 2273.26,980.458 2273.9,980.458 2274.55,980.458 \n",
       "  2275.19,980.458 2275.84,980.458 2276.49,980.458 2277.13,980.458 2277.78,980.458 2278.43,980.458 2279.07,980.458 2279.72,980.458 2280.37,980.458 2281.01,980.458 \n",
       "  2281.66,980.458 2282.3,980.458 2282.95,980.458 2283.6,980.458 2284.24,980.458 2284.89,980.458 2285.54,980.458 2286.18,980.458 2286.83,980.458 2287.48,980.458 \n",
       "  2288.12,980.458 2288.77,980.458 2289.41,980.458 2290.06,980.458 2290.71,980.458 2291.35,980.458 2292,980.458 2292.65,980.458 2293.29,980.458 2293.94,980.458 \n",
       "  2294.59,980.458 \n",
       "  \"/>\n",
       "<path clip-path=\"url(#clip6800)\" d=\"\n",
       "M1853.56 386.635 L2280.76 386.635 L2280.76 205.195 L1853.56 205.195  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<polyline clip-path=\"url(#clip6800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1853.56,386.635 2280.76,386.635 2280.76,205.195 1853.56,205.195 1853.56,386.635 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip6800)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1877.56,265.675 2021.56,265.675 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip6800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 2045.56, 283.175)\" x=\"2045.56\" y=\"283.175\">Dynamic</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip6800)\" style=\"stroke:#e26f46; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1877.56,326.155 2021.56,326.155 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip6800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 2045.56, 343.655)\" x=\"2045.56\" y=\"343.655\">FEA</text>\n",
       "</g>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 100,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "setup=getSetup()\n",
    "displacementFEA=fea(setup)\n",
    "\n",
    "latticeSize=4\n",
    "numTimeSteps=3000\n",
    "setup=getSetup()\n",
    "displacements=[]\n",
    "save=true\n",
    "returnEvery=1\n",
    "runMetavoxelGPU!(setup,numTimeSteps,latticeSize,displacements,returnEvery,false)\n",
    "\n",
    "numTimeStepsRecorded=length(displacements)\n",
    "d=[]\n",
    "dFEA=[]\n",
    "j=length(displacements[end])\n",
    "step=1\n",
    "for i in 1:step:numTimeStepsRecorded\n",
    "    append!(d,displacements[i][j].y)\n",
    "    append!(dFEA,displacementFEA[j].y)\n",
    "end\n",
    "println(\"FEA displacement= $(displacementFEA[j].y),converged displacement= $(displacements[numTimeSteps][j].y)\")\n",
    "plot(1:step:numTimeStepsRecorded,d,label=\"Dynamic\",xlabel=\"timestep\",ylabel=\"displacement\",title=\"Convergence Study\")\n",
    "plot!(1:step:numTimeStepsRecorded,dFEA,label=\"FEA\")\n",
    "# savefig(\"convergence\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip5000\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<path clip-path=\"url(#clip5000)\" d=\"\n",
       "M0 1600 L2400 1600 L2400 0 L0 0  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip5001\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<path clip-path=\"url(#clip5000)\" d=\"\n",
       "M235.533 1487.47 L2352.76 1487.47 L2352.76 47.2441 L235.533 47.2441  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip5002\">\n",
       "    <rect x=\"235\" y=\"47\" width=\"2118\" height=\"1441\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  295.454,1487.47 295.454,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  794.799,1487.47 794.799,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1294.14,1487.47 1294.14,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1793.49,1487.47 1793.49,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2292.83,1487.47 2292.83,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  235.533,1396.74 2352.76,1396.74 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  235.533,960.492 2352.76,960.492 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  235.533,524.249 2352.76,524.249 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  235.533,88.0053 2352.76,88.0053 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,1487.47 2352.76,1487.47 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,1487.47 235.533,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  295.454,1487.47 295.454,1470.19 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  794.799,1487.47 794.799,1470.19 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1294.14,1487.47 1294.14,1470.19 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1793.49,1487.47 1793.49,1470.19 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2292.83,1487.47 2292.83,1470.19 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,1396.74 260.939,1396.74 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,960.492 260.939,960.492 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,524.249 260.939,524.249 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  235.533,88.0053 260.939,88.0053 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 295.454, 1541.47)\" x=\"295.454\" y=\"1541.47\">4</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 794.799, 1541.47)\" x=\"794.799\" y=\"1541.47\">5</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1294.14, 1541.47)\" x=\"1294.14\" y=\"1541.47\">6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1793.49, 1541.47)\" x=\"1793.49\" y=\"1541.47\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2292.83, 1541.47)\" x=\"2292.83\" y=\"1541.47\">8</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 211.533, 1414.24)\" x=\"211.533\" y=\"1414.24\">-0.003</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 211.533, 977.992)\" x=\"211.533\" y=\"977.992\">-0.002</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 211.533, 541.749)\" x=\"211.533\" y=\"541.749\">-0.001</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 211.533, 105.505)\" x=\"211.533\" y=\"105.505\">0.000</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip5002)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  295.454,88.0053 794.799,88.0053 1294.14,90.0895 1793.49,88.0053 2292.83,1446.71 \n",
       "  \"/>\n",
       "<path clip-path=\"url(#clip5000)\" d=\"\n",
       "M1947.14 251.724 L2280.76 251.724 L2280.76 130.764 L1947.14 130.764  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1947.14,251.724 2280.76,251.724 2280.76,130.764 1947.14,130.764 1947.14,251.724 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip5000)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1971.14,191.244 2115.14,191.244 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip5000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 2139.14, 208.744)\" x=\"2139.14\" y=\"208.744\">FEA</text>\n",
       "</g>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 94,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "n=[]\n",
    "nFEA=[]\n",
    "j=length(displacements[end])\n",
    "for i in 4:j\n",
    "    append!(n,displacements[end][i].y)\n",
    "    append!(nFEA,displacementFEA[i].y)\n",
    "end\n",
    "plot(4:j,nFEA,label=\"FEA\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "@webio": {
   "lastCommId": null,
   "lastKernelId": null
  },
  "kernelspec": {
   "display_name": "Julia 1.2.0",
   "language": "julia",
   "name": "julia-1.2"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "1.2.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}