30:T6c1,universe = { const renderer = new THREE.WebGLRenderer({antialias: true, alpha:true}); renderer.setSize(viewSize, viewSize); renderer.setPixelRatio(window.devicePixelRatio); renderer.domElement.style["background"] = `url(${bgImageURL}) no-repeat center center`; renderer.domElement.style["background-size"] = "cover"; const PI = Math.PI; const deg2rad = deg => PI*deg/180; const controls = new THREE.OrbitControls(camera, renderer.domElement); controls.target.set(0,0,camera.position.z+1); controls.maxPolarAngle = deg2rad(90+maxOrbitAngleDeg); controls.minPolarAngle = deg2rad(90-maxOrbitAngleDeg); controls.maxAzimuthAngle = deg2rad(180+maxOrbitAngleDeg); controls.minAzimuthAngle = deg2rad(180-maxOrbitAngleDeg); controls.rotateSpeed = -1; controls.enableRotate = true; controls.enableZoom = false; controls.enablePan = false; controls.addEventListener("change", () => { renderer.render(scene, camera); }); invalidation.then(() => (controls.dispose(), renderer.dispose())); applyCosmology(initialModel); renderer.render(scene, camera); return Object.assign(renderer.domElement, { updateVisibility(viz) { for(const T of visibleTypes) scene.getObjectByName(T).visible = viz.includes(T); renderer.render(scene, camera); }, updateView(cz) { // Linear ramp from 0.5 to 1.0 Gpc. const alpha = Math.max(0, Math.min(1, (cz/1000 - 0.5) / 0.5)); renderer.setClearAlpha(alpha); camera.position.set(0,0,cz); controls.target.set(0, 0, camera.position.z+1); controls.update(); renderer.render(scene, camera); }, updateCosmology(model) { applyCosmology(model); renderer.render(scene, camera); } }); }31:T610,class SceneGalaxy extends SceneObj { constructor(name, coords, minRedshift, maxRedshift, size) { super(name, coords, minRedshift, maxRedshift); const nChunks = textures[name].length; this.chunks = arange(nChunks + 1, {dtype:Int32Array}).map(n => Math.round(n*this.count/nChunks)); this.shuffle = RNG.shuffle(arange(this.count, {dtype:Int32Array})); this.xyzs = []; for(let j=0; j\"\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":927,"value":"observerText = {\n const D = 1e-3 * cameraDistance;\n const f = cameraDistance / universeSize;\n const lbtime = cosmologyModel.lookback_time(0, cameraRedshift);\n return html`You are at redshift z=${Number(cameraRedshift).toFixed(2)}, which is ${Number(D).toFixed(2)} Gpc from Earth and ${Number(100*f).toFixed(0)}% of the way to the edge of the observable universe. Light reaching Earth from here today left ${Number(lbtime).toFixed(1)} Gyrs ago.`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":874,"value":"modelText = {\n const curv = (cosmologyModel.R0 > 400000) ? \"spatially flat\" :\n ((cosmologyModel.kappa < 0) ? \"negatively\":\"positively\") + \" curved\";\n return html`This universe is ${Number(cosmologyModel.age).toFixed(2)} billion years old and ${curv}.`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":273,"value":"viewof showTypes = Checkbox(visibleTypes, {\n value:visibleTypes,\n format:T=>html`${T}`\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":47,"value":"viewof Om0 = Range([0.2,0.4], {value:Om0Default, step:0.005, label:html`Matter ${tex`\\Omega_{\\text{m},0}`}`})","pinned":false,"hidden":false,"mode":"ojs"},{"id":63,"value":"viewof Ode0 = Range([0.6,0.8], {value:Ode0Default, step:0.005, label:html`Dark energy ${tex`\\Omega_{\\Lambda,0}`}`})","pinned":false,"hidden":false,"mode":"ojs"},{"id":67,"value":"reset = {\n // See https://observablehq.com/@mbostock/views-are-mutable-values\n const form = html`\n`;\n form.benchmark.onclick = () => {\n viewof Om0.value = Om0Default; viewof Om0.dispatchEvent(new Event('input'));\n viewof Ode0.value = Ode0Default; viewof Ode0.dispatchEvent(new Event('input'));\n }\n return form;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":802,"value":"md`### Legacy Survey Image`","pinned":false,"hidden":false,"mode":"ojs"},{"id":819,"value":"fieldOfViewDeg = 4","pinned":true,"hidden":false,"mode":"ojs"},{"id":822,"value":"imageSize = viewSize","pinned":true,"hidden":false,"mode":"ojs"},{"id":805,"value":"imageURL = {\n const RA = 180 * catalog[\"RA0\"] / Math.PI;\n const DEC = 180 * catalog[\"DEC0\"] / Math.PI;\n const pixscale = 3600 * fieldOfViewDeg / imageSize;\n const url = new URL(\"https://www.legacysurvey.org/viewer/jpeg-cutout\");\n url.searchParams.set(\"ra\", Number(RA).toFixed(3));\n url.searchParams.set(\"dec\", Number(DEC).toFixed(3));\n url.searchParams.set(\"layer\", \"unwise-neo6\");\n url.searchParams.set(\"pixscale\", Number(pixscale).toFixed(5));\n url.searchParams.set(\"size\", imageSize);\n return url;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":829,"value":"html`Use this link to view the Legacy survey image for this data.`","pinned":false,"hidden":false,"mode":"ojs"},{"id":868,"value":"// blocked by CORS policy\n//imageData = await fetch(imageURL)","pinned":true,"hidden":false,"mode":"ojs"},{"id":1255,"value":"md`Use the version hosted on github instead of a FileAttachment`","pinned":true,"hidden":false,"mode":"ojs"},{"id":1205,"value":"bgImageURL = \"https://dkirkby.github.io/desi3d/img/rosette15-unwise-neo6.jpg\"\n//bgImageURL = await FileAttachment(\"rosette15-unwise-neo6.jpg\").url()","pinned":true,"hidden":false,"mode":"ojs"},{"id":12,"value":"md`### Geometry`","pinned":false,"hidden":false,"mode":"ojs"},{"id":759,"value":"viewSize = Math.min(width, 800)","pinned":true,"hidden":false,"mode":"ojs"},{"id":22,"value":"md`### 3D Rendering`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1029,"value":"maxOrbitAngleDeg = 15;","pinned":true,"hidden":false,"mode":"ojs"},{"id":24,"value":"camera = {\n const fov = 25;\n const aspect = 1;\n const near = 25;\n const far = 8000;\n const camera = new THREE.PerspectiveCamera(fov, aspect, near, far);\n camera.position.set(0,0,0);\n camera.lookAt(new THREE.Vector3(0,0,10000));\n return camera;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":319,"value":"sourceColor = ({\n BGS: \"0xffff00\",\n LRG: \"0xff0000\",\n ELG: \"0x00ff00\",\n QSO: \"0x00ffff\",\n LYA: \"0xddffff\"})","pinned":true,"hidden":false,"mode":"ojs"},{"id":467,"value":"visibleTypes = Object.keys(sourceColor)","pinned":true,"hidden":false,"mode":"ojs"},{"id":1108,"value":"md`File attachments must be named statically so we list them all here`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1269,"value":"textureURLs = {\n const numTextures = 4;\n const base = \"https://dkirkby.github.io/desi3d/img/\";\n return Object.fromEntries(Types.map(ttype => [\n ttype,\n Array.from({length:numTextures}, (_,i) => base + ttype + i + \".png\")\n ]));\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1268,"value":"textures = {\n const loader = new THREE.TextureLoader();\n const load = url => new Promise(async resolve => loader.load(url, resolve));\n const entries = Object.entries(textureURLs).map(async ([key, urls]) => {\n const textures = await Promise.all(urls.map(load));\n return [key, textures];\n });\n return Object.fromEntries(await Promise.all(entries));\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":838,"value":"class SceneObj {\n constructor(name, coords, minRedshift, maxRedshift) {\n // Find the indices passing the redshift cuts, assuming that the input coords are already sorted.\n const last = coords[\"z\"].filter(z => z <= maxRedshift).length;\n const first = coords[\"z\"].slice(0,last).filter(z => z <= minRedshift).length;\n this.unitX = Float32Array.from(coords[\"phi\"].slice(first,last), phi => Math.cos(phi));\n this.unitY = Float32Array.from(coords[\"phi\"].slice(first,last), phi => Math.sin(phi));\n this.theta = Float32Array.from(coords[\"theta\"].slice(first,last), mrad => 1e-3*mrad);\n this.redshift = Float32Array.from(coords[\"z\"].slice(first,last));\n this.count = this.redshift.length;\n this.obj = new THREE.Group();\n this.obj.name = name;\n }\n needsUpdate() {\n for(const child of this.obj.children) {\n child.geometry.attributes.position.needsUpdate = true;\n }\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1178,"value":"RNG = new RandomState('desi3d');","pinned":true,"hidden":false,"mode":"ojs"},{"id":488,"value":"$31","pinned":true,"hidden":false,"mode":"ojs"},{"id":593,"value":"$32","pinned":true,"hidden":false,"mode":"ojs"},{"id":562,"value":"SceneObjects = ({\n BGS: new SceneGalaxy(\"BGS\", catalog[\"BGS\"], 0.1, 0.5, 2),\n LRG: new SceneGalaxy(\"LRG\", catalog[\"LRG\"], 0.3, 1.4, 6),\n ELG: new SceneGalaxy(\"ELG\", catalog[\"ELG\"], 0.4, 1.9, 6),\n QSO: new SceneGalaxy(\"QSO\", catalog[\"QSO\"], 0.9, 5.0, 6),\n LYA: new SceneLineOfSight(\"LYA\", catalog[\"QSO\"], 2.1, 5.0)\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":1292,"value":"SceneObjects[\"BGS\"].count + SceneObjects[\"LRG\"].count + SceneObjects[\"ELG\"].count + SceneObjects[\"QSO\"].count","pinned":true,"hidden":false,"mode":"ojs"},{"id":1042,"value":"createRings = function(distances, radius, nSegments=48, color=0xffffff) {\n const nRings = distances.length;\n const xyz = zeros(6 * nSegments * nRings, {dtype:Float32Array});\n const dPhi = 2 * Math.PI / nSegments;\n let [x1, y1, x2, y2] = [0, 0, radius, 0];\n for(let i=0; i cosmo.Sk(d));\n for(const T of visibleTypes) SceneObjects[T].updateComoving(distanceGrid, transverseGrid);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":19,"value":"md`### Catalog`","pinned":false,"hidden":false,"mode":"ojs"},{"id":89,"value":"md`Load a catalog of observed angles and redshifts.`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1305,"value":"catalog = JSON.parse(\n await blob2text(\n await png2blob(\"https://dkirkby.github.io/desi3d/img/aleatoire.png\")))","pinned":true,"hidden":false,"mode":"ojs"},{"id":280,"value":"Types = catalog[\"ttypes\"]","pinned":true,"hidden":false,"mode":"ojs"},{"id":1281,"value":"counts = Types.map(T => catalog[T][\"z\"].length).reduce((a, b) => a + b, 0)","pinned":true,"hidden":false,"mode":"ojs"},{"id":128,"value":"md`Find the maximum redshift in the catalog`","pinned":false,"hidden":false,"mode":"ojs"},{"id":131,"value":"redshiftMax = Math.max(...Types.map(T => Math.max(...catalog[T][\"z\"])))","pinned":true,"hidden":false,"mode":"ojs"},{"id":1296,"value":"// Copied from @mootari/data-images to avoid a separate js fetch at page load\nasync function png2blob(src, {channels = 3} = {}) {\n const img = DOM.element('img', {crossorigin: 'anonymous', src});\n await img.decode();\n const {naturalWidth: width, naturalHeight: height} = img;\n const ctx = DOM.element('canvas', {width, height}).getContext('2d');\n ctx.drawImage(img, 0, 0);\n return decode(ctx.getImageData(0, 0, width, height), {channels}); \n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1298,"value":"// Copied from @mootari/data-images to avoid a separate js fetch at page load\nfunction decode(imgData, {channels = 4, eofMarker = EOF} = {}) {\n const d = imgData.data;\n const arr = new Uint8Array(d.length / 4 * channels);\n for(let i = 0, j = 0; i < d.length; (i += 4), (j += 4 - channels)) {\n arr.set(d.slice(i, i + channels), i - j);\n }\n\n let len = arr.length;\n while(len > 0 && arr[len-1] !== eofMarker) len--;\n return new Blob([arr.slice(0, len-1)]);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1302,"value":"// Copied from @mootari/data-images to avoid a separate js fetch at page load\nasync function blob2text(blob) {\n return new Response(blob).text();\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1300,"value":"EOF = 123","pinned":true,"hidden":false,"mode":"ojs"},{"id":5,"value":"md`---`","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"THREE = {\n // Using the recipe in https://observablehq.com/@observablehq/how-to-require-stubborn-modules\n const THREE = window.THREE = await require(\"three@0.127.0/build/three.min.js\");\n await require(\"three@0.127.0/examples/js/controls/OrbitControls.js\").catch(() => {});\n await require(\"three@0.127.0/examples/js/lines/LineMaterial.js\").catch(() => {});\n //const THREE = window.THREE = await require('three');\n //await require('three/examples/js/controls/OrbitControls.js').catch(() => {});\n return THREE;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":250,"value":"import {Scrubber} from \"@mbostock/scrubber\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":44,"value":"import {Range,Checkbox} from '@observablehq/inputs'","pinned":false,"hidden":false,"mode":"ojs"},{"id":39,"value":"import {Planck15,FluidsModel} from '@dkirkby/cosmocalc'","pinned":false,"hidden":false,"mode":"ojs"},{"id":147,"value":"import {zeros,linspace, interp, arange} from \"@dkirkby/jumpy\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1175,"value":"import {RandomState} from \"@dkirkby/random\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":189,"value":"import {interleave} from \"@dkirkby/array-utilities\"","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}