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 it works as follows:
- 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
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 it works as follows:
- 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"
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 it works as follows:
- 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"
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 it works as follows:
- 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"
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: 1) 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)
2) 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"
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"
ah ta filozofia ;)
a jakie masz parametry? pewnie nie musisz nic przerabiac. moj program jest bardzo modulowy poza tym.
ja uzywam tak: Omega_b, Omega_CDM, Omega_Lambda, H_0, tau, n i warunki poczatkowe (tylko 2 wart. - adiabatyczne albo isocurv.) liczeniaa jest od groma i jeszcze troche to potrwa. ok. w czym napisany jest ten Twoj program do interpolacji i czy ten sam program robi ta parametryzacje widm (i co to jest za parametryzacja)? Na Twojej stronie w kazdym pliku (w tych o nazwach z wartosciami parametrow) jest chyba z kilkanascie roznych wartosci - czy to jest wlasnie ta parametryzacja ?
teraz mam czas wreszcie zeby sie zabrac do roboty. zastanawiam sie ostatnimi czasy jak to zrobic zeby przypasowac docelowo ok 280 tys widm teoretycznych do danych z mapa. (choc zaczymam myslec ze chyba lepiej bedzie zrobic to kilkuetapowo - metoda kolejnych przyblizen). Jest taki programik co liczy prawdopodobienstwo pasownosci danego widma do danych mapa. na stronach lambda. obczajales go moze ?
I jeszcze jedna kwestia. Jak wiadomo wielka bolaczka sa dlugotrwale obliczenia widm. ale ten proces mozna znacznie przyspieszyc mniej wiecej o rzad stosujac nowa metode ksplit - pewnie jak pisales magistere to jej pewnie jeszcze nie bylo (a moze byla). jesli ja uzywalas to co o niej sadzisz ? (daje dosc mala mozliwosc ingerencji w rozne parametry)
pozdrawiam Bartek
ps. imponujace podroze na Twojej www....
moja www to juz historia, from mysterious reasons. moze ja przeniose na astronomie.
fortran compiler we have been talking about is at: http://www.psrv.com/linux.html
good luck ;-)
Michal Frackowiak wrote:
fortran compiler we have been talking about is at: http://www.psrv.com/linux.html
good luck ;-)
@#$(%&#!!!! "This product is no longer being offered." a nice thing about non-opensource soft.
huehue, thats right ;) I've been there. dont know what happened is it commerical from some time on or what ? but there are some trial verisons over the net but need to get some licence key from different server.
******************************************************************************* Bartosz Lew e-mail: szajtan@astri.uni.torun.pl szajtan@phys.uni.torun.pl Centrum Astronomii UMK szajtan@poczta.onet.pl TORUN, POLAND www: http://www.phys.uni.torun.pl/~szajtan - no longer *******************************************************************************
On Fri, 14 Mar 2003, Michal Frackowiak wrote:
Michal Frackowiak wrote:
fortran compiler we have been talking about is at: http://www.psrv.com/linux.html
good luck ;-)
@#$(%&#!!!! "This product is no longer being offered." a nice thing about non-opensource soft.
There is some draft description of the code. For anyone interested.... - perhaps only Bartek ;-)
On Mon, 17 Mar 2003, Michal Frackowiak wrote:
There is some draft description of the code. For anyone interested.... - perhaps only Bartek ;-)
Hey, I'm interested too! I just don't have time (yet).
I'd rather do a circles analysis. Look up
-
Bartosz Lew -
Boud Roukema -
Michal Frackowiak -
szajtan -
Szajtan Odwieczny