« first day (159 days earlier)      last day (3070 days later) » 

8:00 PM
@evinda Generally, you pick a constant Courant number (something less than 1) and a particular spatial mesh. From that, you calculate the time step. If the mesh spacing is non-equidistant, you choose your dt from the maximal dx such that the Courant number is (hopefully) minimal.
 
I think Anne's problem was freewheeling her first (non-physics) year, then getting 50 out of 60 credits freewheeling her second year, and then she suddenly had to work for her credits
 
@evinda However, since you have a local advection velocity that is dependent on both space and time (2*t*x^2), the minimal Courant number is not only dependent on the spacing both also the local values, but that shouldn't matter if you choose a low enough Courant limit, which requires choosing, checking, and reducing (if necessary).
 
I have a former classmate from high school, who started physics with me, than got attrited. Now she's trying to restart, and she's failing again.
 
@AndrasDeak I can imagine. Anne's 23, in a class full of 17 year olds...
 
8:01 PM
@Adriaan Freewheeling? Is that another term for just not trying
 
I see... so isn't it like that?

I have thought that we could use the upper bounds for x and t. Then 2*νt_nx_j^2<= 8*ν<=1 => ν<=1/8.
 
@Ballbreaker I think he means passing without an effort. Isn't that right?
 
@Ballbreaker no, I meant it in the sense that she got all credits without effort
@AndrasDeak Yes.
 
Ahhhhhh alright
I like that term. I freewheeled 99/100ths of university lol
Which still to this day surprises me
 
@AndrasDeak I found that the exact solution is the following and I have checked if it satisfies the pde:
 
8:04 PM
@Ballbreaker that's what I did at Earth Sciences, got 82% average, then in physics got 61% average (a)
 
@evinda I'm glad it hasn't changed since yesterday.
 
@Adriaan Yeah I ended with around an 80% in engineering physics
 
@Adriaan rock lover
 
Doing that, but then again.. I think I'm kind of an outlier in that sense
 
I have implemented it like that pastebin.com/Z70BGvcx but I don't find an error. Is it wrong? @AndrasDeak
 
8:05 PM
@AndrasDeak actually 6 out of 6 poems for me yesterday were about me bringing my geological hammer everywhere to smash rocks :D
can't bring it to Budapest though, since I only have hand-luggage
 
@Adriaan Just put the hammer in your other hand!
I'm sure they'll let you bring that on the plane
(That was such an Andras joke)
 
@Ballbreaker please tell me it wasn't
 
OH it was.
It was.
 
I even almost got arrested on the ferry from Newcastle-upon-Tyne to IJmuiden, though I'm not sure it was the 1.5L whisky where 1L was allowed or the 4kg rocks, some out of national parks, or the hammer :P
 
8:08 PM
@Adriaan FWIW we don't have many rocks here in the city
maybe on Gellert hill, but that should be sedimentary stuff, nothing of interest
 
@TroyHaskin For N_x=100, N_t=1000 I get the following:
 
@Adriaan should've been more careful with those RFID national park rocks
 
@AndrasDeak You're delusional if you think that that joke wasn't an Andras joke.
 
@Ballbreaker of course I am:D
 
8:09 PM
@AndrasDeak that's more than we have in this piece of shit we call a country :P
 
@Adriaan Lots of rocks in Canada!
Rocks as far as the eye can see
 
the Netherlands consists of the crap Germany and France secrete into the Rhine/Meuse rivers
@Ballbreaker I like rocks.
 
@evinda I still don't like how your exact function is so steep
if it is, that won't behave well with the numerical solution
 
@AndrasDeak I like the potential well. Beautiful solution to the wave equation in simple terms :)
 
Did you see the function where I defined the exact solution? Isn't it right? What else could I have done wrong? @AndrasDeak
 
8:12 PM
@Adriaan potential well? solution to the wave eq?
 
@Adriaan I quantum tunneled from Canada to the Nederland to smash up all your rocks
 
@evinda I can't really imagine how that function looks like. I'm just saying that if that function really is the exact one, then your PDE is a weird one
 
@AndrasDeak pff, did you ever learn Quantum Mechanics? It's like the first or second example you get :P
 
@Ballbreaker well, it's possible...
 
@Ballbreaker they're 3km underneath the surface... :P
 
8:13 PM
@Adriaan OK:P
Oooh, I think I know what you mean
particles in a box
yeah, that's cute
 
The finite potential well (also known as the finite square well) is a concept from quantum mechanics. It is an extension of the infinite potential well, in which a particle is confined to a box, but one which has finite potential walls. Unlike the infinite potential well, there is a probability associated with the particle being found outside the box. The quantum mechanical interpretation is unlike the classical interpretation, where if the total energy of the particle is less than the potential energy barrier of the walls it cannot be found outside the box. In the quantum interpretation, there...
you learn the infinite one first though
SO PRETTY
 
it's just sin(k_x x)sin(k_y y), innit?
 
@AndrasDeak Damn right it's possible
 
I always figured it would be insane to think that there is a probabilistic chance that you might just disappear one day
 
8:16 PM
@AndrasDeak "The wavefunction of a 2D well with nx=4 and ny=4" is what the caption says. You're right, though you need to include the 1/sqrt(L1*L2) as well
 
That it may have already happened to someone/something on our planet, and we would just never know
 
@evinda the functions seems fairly symmetrical, not like your one
how do you call the exact2 function? You might be screwing up something
@Ballbreaker disappear, no
relocate, yes
possible, yes
plausible, hell no
even a single particle finds it almost impossible to move further than a few nanometers
and you have 60*1000*6e23 particles, which would have to move at once
we don't have a notation system to express the improbability of anything quantum happening on the macroscale:)
@Adriaan well normalization is trivial
 
I loved Feinman's quote: "There might just be one electron, piercing the present plane almost infinitely many times, therefore we perceive it as lots of electrons"
 
@Adriaan yeah that's a bit far out
although Hawking radiation makes me wonder
 
@AndrasDeak not if it concerns the Yang-Mills equations. You'd get the Nobel prize for that, if you'd been earlier than Gerard 't Hooft ;)
He actually still teachers particle physics in Utrecht
 
8:20 PM
@Adriaan cool
 
Im highly tempted to follow his course
 
@Adriaan you probably should
rare chance
 
I will do elementary particle physics in Delft, since Utrecht is more than an hour by public transport...
not really fun for a single class
 
:(
yeah, I can imagine
 
I did hear a talk by him once though, together with his PhD supervisor, Martinus Veltman. Veltman also got that nobel prize, and was the one who invented Mathematica
Serge Haroche spoke the year after that on "our" congress of the Dutch physical society
of which I am a proud member of course :P
 
8:22 PM
@AndrasDeak Well that's what I meant by disappear lol, thanks captain obvious
 
@Adriaan sure:P
@Ballbreaker not my fault if you don't say what you mean to say:P
 
As far as an earthly observer would know, the person disappeared and didn't relocate..
 
Robertus Henricus "Robbert" Dijkgraaf (Dutch pronunciation: [ roːˈbɛrtʏs ɦɛnˈrikʏs ˈrɔbərt ˈdɛikɣraːf]; born 24 January 1960) is a Dutch mathematical physicist and string theorist. He is tenured professor at the University of Amsterdam, and director and Leon Levy professor at the Institute for Advanced Study in Princeton. == Biography == Robertus Henricus Dijkgraaf was born on 24 January 1960 in Ridderkerk, Netherlands. He lives in Princeton, New Jersey. Dijkgraaf is married to the author Pia de Jong and has three children. Dijkgraaf went to Erasmiaans Gymnasium in Rotterdam, Netherlands. He...
 
@Ballbreaker so relocate where, smartypants?
> Dutch pronunciation: [ roːˈbɛrtʏs ɦɛnˈrikʏs ˈrɔbərt ˈdɛikɣraːf]
 
@AndrasDeak I'm an earthly observer, how should I know ;)
 
8:23 PM
*shiver*
 
@AndrasDeak he's the most awesome Dutch physicist I know of. I actually spoke to him a few times when he had talks here
 
@AndrasDeak lmao I was going to comment on that too
 
@AndrasDeak I'll pronounce it for you Saturday, you can do your mathematicians :P
 
@Adriaan Awwww, so cute :) <3
 
@Adriaan you mean Erdős and Rényi?
those are easy
you have impossible sounds
 
8:26 PM
I wonder if Shrodinger believed in Condoms
or if he was just hoping the whole pregnancy thing was like a black-box
 
@Ballbreaker "Schrödinger", barbarian
 
I don't got dem fancy letter thingamabobs on my computer
I like that google didn't say that "thingamabob" was incorrect spelling.. I guess that's a real word? O.o
 
@Ballbreaker you don't know if you came until you pull it out from the box?
 
thingamajig and thingamabob are words.. wow
 
and doohickey?
 
8:28 PM
lmao yes
 
@AndrasDeak I have found a mistake and changed it... Now I get the following graph:
 
@AndrasDeak The baby doesn't exist until I observe it
(best pro-abortion argument, ever)
 
@Ballbreaker As they say, the baby is in the eye of the beholder.
@evinda there you go
 
@AndrasDeak So is it right?
 
8:30 PM
and as I always say "A baby doesn't have a soul until it has a shadow"
akin to "Are all fetus' vampires, because they don't have a reflection" ?
 
@evinda well it's much better:P
 
@Ballbreaker and live off humans ;)
 
Is your final T the same in all graphs?
 
EXACTLY
I guess technically all undead don't have shadows right? So I guess a baby is really just unborn
the unborn and undead don't have shadows..
 
@Ballbreaker so you can still kill it if it comes out, as long as the room is well lit from all sides
 
8:31 PM
lmao
 
Why if we increase N_x and N_t, is the difference between the exact solution and the approximation growing? @AndrasDeak
 
@Ballbreaker put some LEDs on this baby and you're ready:
@evinda my previous question is related to your question
I asked first, so you answer first
@Ballbreaker that last one might be a bit too morbid to star:P
 
@AndrasDeak CATERPILLAR
I like Caterpillar trucks :D
 
@Adriaan is it?
I like them too; this is just the first pretty hit for wood chipper
 
Not sure, but it's yellow with a big bucket :p
 
8:34 PM
@Adriaan it says TRX CARTON babyremover 1260
 
@AndrasDeak hahahahahahahahahahahahahahahahahaha
 
I meant these beauties:
 
@Adriaan are you familiar with CAT phones?
 
such a pretty truck, such a pretty rock :D
 
Oh I didn't see it.. Yes, it's 1but I also plot the exact solution and the approximation in intermediate time steps.

Sorry, we don't increase N_x and N_t. We just increase the value of the time...
 
8:35 PM
@AndrasDeak nope, only CAT dozers, scrapes, diggers and trucks :p
 
@evinda that's why I asked
if the time is increasing, it's understandable that the solution diffuses away
 
I like that the output of that woodchipper is pointed back on you, just incase you want to cover yourself in your mistakes
 
because of the error introduced by the scheme, as Troy told you yesterday
 
@AndrasDeak Why does this happen?
 
errors accumulate
the longer the time, the larger the diffusion and the larger the error
 
8:36 PM
@AndrasDeak and then you get crap babies, because the weaker sperm wins
 
@Adriaan lol:D Not too related, but sure
don't work in a radar tower if you want sons
 
@AndrasDeak radar towers are for pussies. I do my work in a nuclear reactor
 
@Adriaan those are harmless
 
though pilot on intercontinental flights are 30 times worse in terms of radiation
 
it's EM radiation that seems to produce more daughters, at least anecdotally
and yeah, airplanes
and miners, probably
 
8:38 PM
@AndrasDeak MINING ENGINEERING
ROCK HAMMERING
HNNNGGG
 
friend of mine actually has loose asbestos and uranium ore (amongst 250 other samples) scatter across his living room, all well documented :D
@AndrasDeak CATERPILLAR
I like caterpillar
 
asbestos is not fun
 
but where's the dozer background in that
 
@Adriaan it runs android, you can customize it:P
 
8:40 PM
@AndrasDeak No,no,no, you don't understand. It should be pre-built-in :P
Mining equipment!
such good parodies
it's actually the largest man-made machine in existence
mines the brown-coal fields of Western Germany
 
@Adriaan not bad
Don't want to start a war, but I like those big-ass metro tunnel diggers and the rail replacing trains more
 
@AndrasDeak I see... :)
 
they're far more ingenious than the Bagger 288, but the Bagger 288 is just plain big, ugly and a behemoth :D
the best machine in terms of engineering quality is imho the Philips VHS player
 
@Adriaan I realized I've seen that bastard on TV
 
Do you maybe know if in general it holds that if the Courant number is closer to the upper bound of the cfl condition, then the approximation is better ? @TroyHaskin
 
8:45 PM
@AndrasDeak Lol my buddies phone sets off microphones and TVs when he receives text messages
I'm concerned for his testicles
 
@evinda isn't it the other way around? for a given courant number, the upper bound means "as fine a mesh as possible", which then causes smaller error?
I didn't understand this yesterday, and I don't understand it still
@Ballbreaker only be concerned for his sperm
and worst case is probably having a higher chance for daughters
girly sperm is tougher, but slower
but if the faster dude sperms die, then they can proliferate
maybe not slower, but tougher
 
if they were slower then there'd be a higher birth rate for boys
 
@AndrasDeak But at poth of the problems at which I implemented the upwind method, the approximation was better for ν=upper bound of cfl condition.
 
@AndrasDeak I stand by what I said..
@AndrasDeak That sounds like a pretty terrible worst case to me
 
8:49 PM
@AndrasDeak For N_x=50 , N_t=800 for example I get:
@AndrasDeak For N_x=100, N_t=800:
 
@evinda there's no contradiction
this is exactly what I'm saying
you have a finer mesh in the second case, and the finer the mesh, the better the approximation (as long as stability is OK)
ideally you should have as fine a mesh as possible
but stability and memory need and computation time restrict you in practice
besides, look at your youngest dataset: you have like ony a dozen points at your peak
so you have problems resolving the function
you're bound to have a lot of discretization error, even without the diffusive effect of upwind
ideally all your bumpy features should be as well resolved as the blue one
 
These graphs remind me of my physics laboratories with particle stuff and the such forth
 
@Ballbreaker those are the best, since they really are Gaussian peaks
they can be fit pretty well
 
@AndrasDeak Does this hold because if the mesh is finer, there are more points that approximate the exact solution?
 
@Adriaan I never actually saw the original. I just started watching: dude is incredible:D
@evinda I don't understand the question
 
9:02 PM
@AndrasDeak Yes. I want a burger now.
 
"I only thought of bacon!"
lol:D
 
The cheese is oozing!
 
@AndrasDeak In the second case is the approximation better since we have more points for the discretization of x and t ?
 
DAYUM
 
What's really insane is how sophisticated the dude actually is, it must be hard for him to play the primitive ganxta:D
@evinda I believe so, yes. But only for x.
And I can be wrong, I'm trying to make a not-so-educated guess
@TroyHaskin is your man if you want to make sure:P
 
9:07 PM
@AndrasDeak It's not very hard to put up a front of being less intelligent than you are :) I would know..
I HAVE YOU GUYS FOOLED!
 
I could watch John Cleese talk about anything
like sports:
 
@AndrasDeak that one's fantastic!
 
@AndrasDeak Hahaha screw you
 
9:22 PM
I like this guy
lmao this is awesome
Oh god hahahahah this comment
"Islam does give women a lot of rights. And a few lefts, a couple kicks, maybe a lashing or stoning to top it off."
l.m.f.a.o
 
@Ballbreaker "Why do the French have so many civil wars? So they can win one every once in a while" :D
 
9:38 PM
Hahahah ayeah
 
0
A: Reshaping Matlab Matrix

flawrYou can easily use the colon :, e.g. A = [1,2; 3,4; 5,6]; You get the output for A(:) ans = 1 3 5 2 4 6

@flawr this is a subpar newb question, and they need the row-wise reshape:P
I don't even know why this is upvoted
Please vote to close as dupe of stackoverflow.com/questions/15102216/…
or whatever
 
@AndrasDeak What question?
Oh that one.
I added the A(:).', you did not see the final verison I suppose=)
 
whoah I've never seen this gold badger who dupehammered it
@flawr that's not right
they need A.'(:) which is not valid matlab
 
No I think A(:).' is perfectly fine???
 
@flawr no:)
 
9:48 PM
Oh i see.
 
@AndrasDeak what about this one?
27
Q: How do you concatenate the rows of a matrix into a vector in MATLAB?

jimboFor an m-by-m (square) array, how do you concatenate all the rows into a column vector with size m^2 ?

 
@AndrasDeak Now something about the first problem...
The cfl condition was v<=1/2.
I want to apply the method for different v s.
I took also v=0.6 and v=0.8 and I got the following errors:

v=0.6: 828365242873271640000000000000000000000000000000000000000000000000000000000000000000000000000.000000

v= 0.8 : 3060735298198143200000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000.000000
 
@beaker sorry, gold badger already hammered it
 
@AndrasDeak The graph is the following:
 
@AndrasDeak I always was a bit slow ;)
 
9:49 PM
@AndrasDeak Can it be right?
 
@AndrasDeak Why the hell aren't most noobs not reading the docs???
 
@evinda those ridiculously huge integers up there? those can't be right, whatever those are:P
@flawr because they're lazy, and don't give a shit:) Which is why I advise against giving extensive answers: otherwise they'll never learn
@evinda I suggest that you observe the magnitude of your y axis
see anything fishy?
 
I only see something jellyfishy
 
@AndrasDeak What do you mean?
 
@Ballbreaker that's one black box I wouldn't put it into...:P
@evinda what do you think I mean?
Box jellyfish (class Cubozoa) are cnidarian invertebrates distinguished by their cube-shaped medusae. Some species of box jellyfish produce extremely potent venom: Chironex fleckeri, Carukia barnesi and Malo kingi. Stings from these and a few other species in the class are extremely painful and can be fatal to humans. == NomenclatureEdit == "Box jellyfish" and "sea wasp" are common names for the highly venomous Chironex fleckeri. However, these terms are ambiguous, as "sea wasp" and "marine stinger" are sometimes used to refer to other jellyfish. == AnatomyEdit == The medusa form of a ...
 
9:54 PM
@AndrasDeak Kekeke :D
"sea wasp" lol jesus
 
y is defined on a huge interval @AndrasDeak
 
@evinda that's sort of an understatement
you have y ~ 1e96
 
@AndrasDeak What could we wrong knowing that the results are right for v<=0.5...
 
@evinda what do you think can be wrong?
You establish that v<=0.5, then you are surprised that for v>0.5 things go wrong
Are you paying any attention in class?:)
 
@AndrasDeak For v>0/5 the method gets unstable.
So can't it be that the results are right?
 
9:58 PM
@evinda it could be, but they are clearly not
unstable = exactly what you're seeing
 
Can we talk about dead babies again, this conversation got boring
lol <(^_^)>
 
@AndrasDeak So seeing that for v>0.5 the errors are huge, the implementation of the method should be correct. Right?
 
@evinda not necessarily, I'm not sure how exact the Courant condition is
but your implementation seems fine already by comparison with the exact solution
 
10:13 PM
@evinda The CFL limit of 1 for stability being only a function of the temporal and spatial mesh spacing is true for constant velocity advection. You have a non-constant velocity (in fact, one that increases with time and quadratically with space), so the upper Courant limit for stability is decreasing with time.
@evinda Also, you're still not doing what I said which would make things easier: set a spatial mesh size and an upper Courant limit and then compute the temporal spacing from the Courant ratio.
@evinda Also, while these graphs are good, typically results are plotted at the same time for difference mesh densities, not at different times for the same mesh density. The former shows how the upwind diffusion decreases with increasing spatial density.
 
padarray!
 
10:30 PM
@beaker don't tell me there's not a dupe out there
 
there must be
@AndrasDeak how about that?
4
Q: How to zero pad a matlab array?

MorganWhat is the easiest way to (zero) pad a matlab array? i.e. given [1,2,3,4] and length 6 return [1,2,3,4,0,0] Background I have a data array which I would like to apply a windowing function to before running fft on the data. I use to pass data directly to fft which would zero pad to the next po...

 
I saw that, but this one's more specific
2
A: How to fill the matrix with 0 in matlab

rayryengIf you have the image processing toolbox, use padarray: >> A = [1 2; 3 4]; >> B = padarray(A, [1 1]) B = 0 0 0 0 0 1 2 0 0 3 4 0 0 0 0 0 The first input is the matrix you want to pad, and the second input is how many zeroe...

which is a dupe of
2
Q: How can I put margins in an image?

TrakerI have a binary image of 18x18 pixels and I want to put margins around this image with the purpose of obtaining an image 20x20 pixels. The image is binary and it can be represented by a matrix of 1s and 0s. The 0 pixels are in black colour and the 1 pixels are in white colour. I need to put ma...

 
@AndrasDeak the second one is better... it gives options in case you don't have the image processing toolbox
 
Indeed, and it's more according to protocol to link to the proto-dupe
voted to close as dupe
 
@AndrasDeak ditto
 
10:40 PM
upvoted you while I was there:P
 
i upvoted the newbie since their answer was more general
 
@beaker and because of Sportsmanship:P
(the badge, that is)
 
aaaah
 
11:09 PM
Doesn't it hold that the mesh is finer if and only if v is closer to the upper bound of the cfl condition? @AndrasDeak Or am I wrong?
@AndrasDeak @TroyHaskin @David Why don't we have also more points for the discretization of t?
 
I'm staying out of this!
 
@David LOL ;)
 
@evinda The CFL condition guarantees stability of the upwind method, but with a fixed CFL number, the diffusive nature of the upwind method can be mitigated by a denser spatial mesh (per that PDF I posted yesterday).
 

« first day (159 days earlier)      last day (3070 days later) »