01-15-2010, 11:36 AM
It is really interesting. I remember taking numerical methods classes and just assuming that only PhD students cared about such things. You are correct that looking into the equations is hard. I spend a lot of time reading papers on applications of X equation to Y problem. Trying to understand them sometimes forces me to go back to simpler things and work my way forward (Looking at poissions equation for our application the first time was a little intimidating:
and as a side note the book "Div, Grad, Curl and all that" is a good reference)
I have never actually heard of FREEFEM but I just looked at it briefly just now, and I suspect it runs like crap. Actually, it looks like something I might get, and have to work on to make run better as part of my job. Stuff like these arithmetic if statements within loops runs absolutely terrible!
// we remove the boundary link
for (int k=0;k<nt;k++)
for (int j=0;j<3;j++)
if (TonBoundary[k] & MaskEdge[j])
{
int s0,s1;
VerticesNumberOfEdge(triangles[k],j,s0,s1);
// cout << s0 << " " << s1 << " ld " << linkd[s0] << " " << linkd[s1] << " lg " << linkg[s0] << " " << linkg[s1] << " apres " ;
linkg[s0] = linkg[s0] != -1 ? -2 : -1;
linkg[s1] = linkg[s1] != -1 ? -2 : -1;
linkd[s1] = linkd[s1] != -1 ? -2 : -1;
linkd[s0] = linkd[s0] != -1 ? -2 : -1;
// cout << " ld " << linkd[s0] << " " << linkd[s1] << " lg" << linkg[s0] << " " << linkg[s1] << endl;
}
What are you doing these days? grad school I suppose? Where at? I am going to try to find a place where I can still do this kind of work, even if I quit this job, I think. It is a lot of fun getting some physics code to run 60% faster, and also staying sharp on any math skills I was starting to forget as well.
and as a side note the book "Div, Grad, Curl and all that" is a good reference) I have never actually heard of FREEFEM but I just looked at it briefly just now, and I suspect it runs like crap. Actually, it looks like something I might get, and have to work on to make run better as part of my job. Stuff like these arithmetic if statements within loops runs absolutely terrible!
// we remove the boundary link
for (int k=0;k<nt;k++)
for (int j=0;j<3;j++)
if (TonBoundary[k] & MaskEdge[j])
{
int s0,s1;
VerticesNumberOfEdge(triangles[k],j,s0,s1);
// cout << s0 << " " << s1 << " ld " << linkd[s0] << " " << linkd[s1] << " lg " << linkg[s0] << " " << linkg[s1] << " apres " ;
linkg[s0] = linkg[s0] != -1 ? -2 : -1;
linkg[s1] = linkg[s1] != -1 ? -2 : -1;
linkd[s1] = linkd[s1] != -1 ? -2 : -1;
linkd[s0] = linkd[s0] != -1 ? -2 : -1;
// cout << " ld " << linkd[s0] << " " << linkd[s1] << " lg" << linkg[s0] << " " << linkg[s1] << endl;
}
What are you doing these days? grad school I suppose? Where at? I am going to try to find a place where I can still do this kind of work, even if I quit this job, I think. It is a lot of fun getting some physics code to run 60% faster, and also staying sharp on any math skills I was starting to forget as well.
2013 Cadillac ATS....¶▅c●▄███████||▅▅▅▅▅▅▅▅▅▅▅▅▅▅▅▅||█~ ::~ :~ :►
2008 Chevy Malibu LT....▄██ ▲ █ █ ██▅▄▃▂
1986 Monte Carlo SS. ...███▲▲ █ █ ███████
1999 F250 SuperDuty...███████████████████►
1971 Monte Carlo SC ...◥☼▲⊙▲⊙▲⊙▲⊙▲⊙▲⊙▲⊙☼◤
2008 Chevy Malibu LT....▄██ ▲ █ █ ██▅▄▃▂
1986 Monte Carlo SS. ...███▲▲ █ █ ███████
1999 F250 SuperDuty...███████████████████►
1971 Monte Carlo SC ...◥☼▲⊙▲⊙▲⊙▲⊙▲⊙▲⊙▲⊙☼◤
