Howdy, Stranger!
It looks like you're new here. If you want to get involved, click one of these buttons!
This forum is in permanent archive mode. Our new active community can be found here.
CTRL+V
PASTE IT! whatever random thing you have copied while reading this, paste it for the world to see. A world of 20 people, that is.
CTRL+V: http://img96.imageshack.us/img96/1102/desktop0pv.jpg
CTRL+V: http://img96.imageshack.us/img96/1102/desktop0pv.jpg
Comments
I already own a projector, but a friend thought that was cool.
I just finished reading the bg post.
WHERE NOT EXISTS
(SELECT * FROM Source
WHERE NOT EXISTS
(SELECT * FROM Target-Source
WHERE Target-Source.Target#=Target.Target#
AND Target-Source.Source#=Source.Source#))
I like the idea, since it will lead to lots of random junk...
This is a model for SQL database division given in a lecture that my friend was asking about.
I basically hit this thread randomly and see if anything's actually in my clipboard ;^)
//get neighborhood
unsigned int neighborhood = hashPhyxel(*i);
//calculate moment matrix A and inverse A^-1
Mat3 moment, invMoment;
//if there are too few neighbors (or if they're coplanar/colinear), the
//moment matrix will be singular (non-invertable)
//if(hashtable[neighborhood].size() < 4) {
// cout << "This is not a good moment: " << hashtable[neighborhood].size()
// << endl;
//}
for(list<Phyxel *>::iterator j = hashtable[neighborhood].begin();
j != hashtable[neighborhood].end(); j++) {
// for(list<Phyxel *>::iterator j = allPhyxels.begin();
// j != allPhyxels.end(); j++) {
if(*(*j) != *i) {
Vec3 diff = (*j)->getCurrentPosition() - i->getCurrentPosition();
Mat3 momentJ = ops::squareVector(diff);
momentJ *= weight(norm(diff),i->getSupportRadius());
moment += momentJ;
}
}
double aDet = invert(invMoment, moment);
if(aDet != 0 ) {
if(i->id == 499) {
//cout << "Phyxel #" << i->id << endl;
//cout << "Moment matrix: " << endl << moment << endl;
//cout << "Determinant: " << aDet << endl;
//cout << "Inverted moment: " << endl << invMoment << endl;
}
Vec3 sumU, sumV, sumW, sumF;
Vec3 dU, dV, dW;
//Equation 15: calculate derivatives
for(list<Phyxel *>::iterator j = hashtable[neighborhood].begin();
j != hashtable[neighborhood].end(); j++) {
//for(list<Phyxel *>::iterator j = allPhyxels.begin();
// j != allPhyxels.end(); j++) {
if(*(*j) != *i) {
Vec3 diff = (*j)->getCurrentPosition() - i->getCurrentPosition();
Vec3 spaceFactor = diff * weight(norm(diff),i->getSupportRadius());
double u = (*j)->getDisplacement()[0] - i->getDisplacement()[0];
sumU += u * spaceFactor;
double v = (*j)->getDisplacement()[1] - i->getDisplacement()[1];
sumV += v * spaceFactor;
double w = (*j)->getDisplacement()[2] - i->getDisplacement()[2];
sumW += w * spaceFactor;
}
}
//get x displacement derivative
dU = invMoment * sumU;
//get y displacement derivative
dV = invMoment * sumV;
//get z displacement derivative
dW = invMoment * sumW;
//assemble the Jacobian matrix J
Mat3 jacobian(dU, dV, dW);
jacobian += ops::IDENTITY;
//if(i->id == 499) {
// cout << "Jacobian: " << endl <<jacobian << endl;
//}
//calculate strain & stress tensors
Mat3 strain = transpose(jacobian) * jacobian - ops::IDENTITY;
Mat3 stress = strain * i->getModulusOfElasticity();
//calculate forces
//Eq 24-25: fixed term
Mat3 externalForce = -2 * i->getVolume() * jacobian * stress;
Mat3 conservation(cross(dV, dW), cross(dW, dU), cross(dU, dV));
Mat3 volumeForce = -(i->getVolume()) * k * (det(jacobian)-1) * conservation;
Mat3 fixedTerm = (externalForce + volumeForce) * invMoment;
//cout << i->id << " volume " << i->getVolume() << endl;
//cout << "External F: " << externalForce << endl;
//cout << "Volume F: " << volumeForce << endl;
//cout << "Force term: " << fixedTerm << endl;
for(list<Phyxel *>::iterator j = hashtable[neighborhood].begin();
j != hashtable[neighborhood].end(); j++) {
//for(list<Phyxel *>::iterator j = allPhyxels.begin();
// j != allPhyxels.end(); j++) {
if((*(*j)) != *i) {
Vec3 diff = (*j)->getCurrentPosition() - i->getCurrentPosition();
Vec3 spaceFactor = diff * weight(norm(diff),i->getSupportRadius());
sumF -= spaceFactor;
(*j)->addExternalForce( fixedTerm * spaceFactor);
}
}
i->addExternalForce(fixedTerm * -sumF);
//reset displacement
i->updatePlasticState();
} else {
//moment matrix couldn't be inverted
//Usually, this means the phyxel was kicked out of range of any
//neighbors. For now, I'm just trusting gravity to bring it back.
i->clearExternalForces();
i->updatePlasticState();
}
*Donut[AFK]* INSULT
*Eurakarte* RETORT
*Donut[AFK]* COUNTER-RETORT
*Eurakarte*QUESTIONING OF SEXUAL PREFERENCE
*Donut[AFK]* SUGGESTION TO SHUT THE FUCK UP
*Eurakarte*NOTATION THAT YOU CREATE A VACUUM
*Donut[AFK]* RIPOSTE
*Donut[AFK]* ADDON RIPOSTE
*Eurakarte* COUNTER-RIPOSTE
*Donut[AFK]* COUNTER-COUNTER RIPOSTE
*Eurakarte*NONSENSICAL STATEMENT INVOLVING PLANKTON
*Miles_Prower* RESPONSE TO RANDOM STATEMENT AND THREAT TO BAN OPPOSING SIDES
*Eurakarte* WORDS OF PRAISE FOR FISHFOOD
*Miles_Prower* ACKNOWLEDGEMENT AND ACCEPTENCE OF TERMS
*this was a conversation that I decided to post ,originally located here on my blog because I found much humor. Never in my life have I giggled excessively because that sums up the life of a true "freak" of society ( you know, the guys with pants dropped to their knees and walk like penguins. Those guys).
13915 EL Moro Ave.
La Mirada, CA 90638
(Lol, just pasted that into Facebook, Myspace, etc.)
TAXONOMY: 2085R0202X - RADIOLOGY - DIAGNOSTIC RADIOLOGY
[copied from address bar, pasted into Facebook]