a jakie masz parametry? pewnie nie musisz nic przerabiac. moj program jest bardzo modulowy poza tym. w_x1 i w_x to do parametryzacji zachowania "kwintesencji" czy raczej pola skalarnego, parametryzuja rownanie stanu dark energy. looknij do mojej magisterki, jesli Cie to interesuje. jest na mojej stronie http://www.ncac.torun.pl/~michalf dobra, ja jestem na kazdym seminarium w camku (czwartek), wiec mozemy sobie uciac pogawedke. fajnie.
nara michal
ps. imponujace podroze na Twojej www....
Bartosz Lew wrote:
to racja, ze generalnie bez sensu jest powtarzac to co zrobili inni juz wczesniej to tak jakby pisac wielki program i nie korzystac z gotowych bibliotek, ale... dwie kewstie:
- to jest moje pierwsze podejscie do problemu i tak czy siak musze sie
troche z tym obyc. Dlatego ... ymm :) tego, wlasnie no, nie ludze sie ze to bedzie jakas wielka odkrywcza praca. (nawaisem mowiac to mi az troche szkoda tak swietnych danych (WMAP'a) do pracy ktora musze zakonczyc w przeciagu gora 3 miesiecy)
- ja zapuszaczjac cmbfast robie siatke troche innych parametrow niz Ty.
ja mam ich w sumie 7 (no powiedzny 6). Ty zapodales jeszcze w_x1 i w_x ( co to w ogole jest ?) tak wiec i tak musialbym przerobic Twoj kod.
za trzeciej strony pisanie programow (choc w wciaga) jest jadnak strata czasu (ktorego mam nie za wiele). Chetnie wpadne pogadac, najlepiej w jakis czwartek na seminarum. (tylko jak bedzie cos ciekaweego ;) - dam znac.
pozdrawiam Bartek
ps O kwintesencji wiem nie duzo (i generalnie to nie jest zwiazane z moim tematem - choc zadanienie jest b. ciekawe), ale jesli chodzi o degeneracje w przestrzeni parametrow to przeciez normalka.
powiem Ci jako czlowiek moze bez duzego doswiadczenia,ale juz troche obyty - bez sensu, ze chcesz to sam oprogramowywac. masa ludzi juz to zrobila. tyle, ze nie udostepniaja kodow - choc jakbys poszukal, to bys ze 2 znalazl (i to wiem jakie). nie ma sensu wywazac otwartych drzwi - owszem, zobacz jak dziala cmbfast, poczytaj papiery, ale nie odwalaj czarnej roboty jako kolejny z rzedu facet. i nie zebym zachwalal moj kod, ale jest on chwilowo najciekawszy. zobacz DASH, to go nawet nie skompilujesz, nie mowiac o szybkosci i ze nie uwzglednia pol skalarnych. jest jeszcze jeden taki "szybki interpolator", ale ma te same "zalety". nie ma sensu odwalac tej samej roboty po raz kolejny. i, moze rozczarowywujac Cie troche, nie ma sensu dopasowywac parametrow kwientesencji do cmb. poza tym, ze juz kopa ludzi to zrobila, z czego 3/4 zle, po prostu sie nie da. dostaniesz tylko zdegenerowane wyniki.
jak chcesz pogadac, to mozemy sie jakos spotkac.
pozdr - michal
Szajtan Odwieczny wrote:
mmm, no tak. juz wiem jak to dziala, mialem do tego zupelnie inne podejscie. Te ok 4000 plikow na Twojej stronie to jest wlasnie ta siatka obliczonych widm ? Ja prawdopodobnie tez bede musial cos takiego napisac - ale poki co startuje powoli z cmbfast. Dlatego dzieki, ale wole sam wszystko oprogramowc :)
pozdrawiam bartek
odpowiedz jest tak prosta, ze az brutalna: interpolacja n-wymiarowa, np. 2-go rzedu, bo 1. troche za kiepsko, patrz wyklad z metod numerycznych i przenies do n-wymiarowej przestrzeni.
boud mowil, ze chcesz cos robic z cmb.... moze przyda Ci sie ten kod.
szajtan wrote:
ok. ciakawa koncepcja ale nie do konca rozumiem jak to dziala. masz obliczona i zparametryzowana (parametrami) w jakis sposob siatke modeli (o roznych parametrach kosmologicznych). tak ze na podstawie zadanych wartosci parametrow (nie kosmolog.) mozna szybko odtworzyc widmo. poszukiwane widmo w n wymiarowej przestrzeni parametrow jest otoczone przez 2^n modeli bazowych (o znanych zestawach wartosci parametrow kosmologicznych, i tych z parametryzacji). Mozna tez obliczyc jak daleko w przestrzeni parametrow kosmologicznych lezy poszukiwany model od modelow budujacych wezly siatki (tylko jak z tego skorzystac). Wiemy rowniez o ile roznia sie miedzy soba zestawy wartosci parametrow przypisanych do wezlow sieci. Ale nie chwytam kewstii obliczenia poszukiwanych wartosci parametrow na podstawie ktorych mozna odtworzyc widmo. Czyli innymi slowy jak to interpolowac skoro wiemy tylko ze poszukiwany zestaw parametrow lezy gdzies w oczku siatki parametrow.?
I have the code - not publicly available yet, but soon
- grid is produced and in nodes exact (well... cmbfast) cl spectrum is
calculated 2. for every node some parametrisation of the spectrum is done (scaling, chebyshev coeficients etc.) 3. for given cosmo parms those parameters are interpolated and the spectrum is reconstructed.
quite simple. but faster than other interpolation methods (for gridded models). that is the main idea...
polarisation also possible, but not yet implemented. I can send you the code with the spectrum database if you wish.
Bartosz Lew wrote:
>ok, what's the catch ? >and where is the code ? > >and btw. how this works ? > >On Mon, 17 Feb 2003, Michal Frackowiak wrote: > > > > > > > > > >>well, not completed, not even documented, but if anyone of you would >>like to use it to reconstruct cmb spectrum REALLY FAST (within >>miliseconds) - please try. >>http://coyote.ncac.torun.pl/~michalf/cmbrecon/cmbrecon.php >> >>well - you can treat it as a preview of what is really about to happen. >>a fast cmb-spectrum reconstruction code. with quintessence (or whatever) >>included (soon also other scalar-tensor things). for generating grid: >>server-client architecture, quite smart reconstruction algorithms.... >>in case of interest/doubts/will to use - feel free to contact me. code >>is quite complete and usable. >> >>michal frackowiak >> >> >> >> >> >> >> >> >> >> >> >------------------------------------------------------------------------ > >LISTNAME: shape-univ >HELP: send an email to sympa@astro.uni.torun.pl with "help" >WEB ARCHIVE: http://www.astro.uni.torun.pl/sympa/shape-univ/ >UNSUBSCRIBE: email to sympa@astro.uni.torun.pl with "unsubscribe > > > > > >
shape-univ"
> > > > > >
LISTNAME: shape-univ HELP: send an email to sympa@astro.uni.torun.pl with "help" WEB ARCHIVE: http://www.astro.uni.torun.pl/sympa/shape-univ/ UNSUBSCRIBE: email to sympa@astro.uni.torun.pl with "unsubscribe shape-univ"
--------------r-e-k-l-a-m-a-----------------
Tanie bilety lotnicze! http://samoloty.onet.pl
LISTNAME: shape-univ HELP: send an email to sympa@astro.uni.torun.pl with "help" WEB ARCHIVE: http://www.astro.uni.torun.pl/sympa/shape-univ/ UNSUBSCRIBE: email to sympa@astro.uni.torun.pl with "unsubscribe shape-univ"
LISTNAME: shape-univ HELP: send an email to sympa@astro.uni.torun.pl with "help" WEB ARCHIVE: http://www.astro.uni.torun.pl/sympa/shape-univ/ UNSUBSCRIBE: email to sympa@astro.uni.torun.pl with "unsubscribe shape-univ"
LISTNAME: shape-univ HELP: send an email to sympa@astro.uni.torun.pl with "help" WEB ARCHIVE: http://www.astro.uni.torun.pl/sympa/shape-univ/ UNSUBSCRIBE: email to sympa@astro.uni.torun.pl with "unsubscribe shape-univ"