/// A Physical Vector Class . //#include #include #include #include #include #include #include using namespace std; using namespace std; typedef long double Real; const Real Pi = 3.14159265359; const Real Empty = 181; inline Real degtorad(Real deg) { return (deg*Pi)/180.0; } inline Real radtodeg(Real rad) { return (rad*180.0)/Pi; } struct Vec { Real i,j,k,mag,dx,dy,dz; //Complex ic,jc,kc,magc,dxc,dyc,dzc; }; //HINSTANCE hInst; //For use in Plot Function - Windows int nCmdShow; Vec vc; /// Vector Output Manipulators const ios_base::fmtflags detailed = (ios_base::fmtflags)0x1; const ios_base::fmtflags simple = (ios_base::fmtflags)0x0; bool __indetail__ = false; bool __read_using_ratio__ = true; ostream& detailed(ostream& os) { __indetail__ = true; return os; } ostream& simple(ostream& sio) { __indetail__ = false; return sio; } ostream& ratiowise(ostream& sio) { __read_using_ratio__ = true; return sio; } ostream& noratiowise(ostream& sio) { __read_using_ratio__ = false; return sio; } /*LRESULT CALLBACK WindowProcedure(HWND hwnd,UINT message,WPARAM wParam,LPARAM lParam) /// Window Procedure for Plot Function. Currently under work, but can be used after uncommenting. /// Currently Prints a vector in 2D with i^ or j^ <= 120. { Vec o = vc; int ox = 200,oy = 185; HFONT hF; HPEN hP; int stepx = static_cast(80/abs(o.i)); if(abs(o.i)>30) stepx=7; int stepy = static_cast(80/abs(o.j)); if(abs(o.j)>30) stepy=7; if(o.j==0) stepy=stepx; if(o.i==0) stepx=stepy; int step = static_cast((stepx + stepy)/2.0); int sx = (o.i>0)?1:-1; int sy = (o.j>0)?1:-1; switch(message) { case WM_PAINT: { HDC hDC; PAINTSTRUCT ps; hDC = BeginPaint(hwnd, &ps); hP = CreatePen(PS_SOLID,1,RGB(220,220,220)); SelectObject(hDC,hP); for(int i=oy;i<=400;i+=step) { MoveToEx(hDC,0,i,NULL); LineTo(hDC,400,i); } for(int i=oy;i>=0;i-=step) { MoveToEx(hDC,0,i,NULL); LineTo(hDC,400,i); } for(int i=ox;i<=400;i+=step) { MoveToEx(hDC,i,0,NULL); LineTo(hDC,i,400); } for(int i=ox;i>=0;i-=step) { MoveToEx(hDC,i,0,NULL); LineTo(hDC,i,400); } hP = CreatePen(PS_SOLID,1,RGB(0,0,0)); SelectObject(hDC,hP); MoveToEx(hDC,0,oy,NULL); LineTo(hDC,400,oy); MoveToEx(hDC,ox,0,NULL); LineTo(hDC,ox,400); hP = CreatePen(PS_SOLID,2,RGB(0,0,255)); SelectObject(hDC,hP); if(abs(o.i)<30&&abs(o.j)<30&&abs(o.i)>0&&abs(o.j)>0) { MoveToEx(hDC,ox,oy,NULL); LineTo(hDC,ox+(step*(o.i)),oy-(step*(o.j))); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,ox+step*(o.i),(oy-step*(o.j))+(sy*5)); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,(ox+step*(o.i))-(sx*5),oy-step*(o.j)); hP = CreatePen(PS_DOT,1,RGB(0,0,0)); SelectObject(hDC,hP); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,ox,oy-step*(o.j)); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,ox+step*(o.i),oy); } else if(o.i==0||o.j==0) { if(o.i==0&&(abs(o.j)<30)) { MoveToEx(hDC,ox,oy,NULL); LineTo(hDC,ox,oy-step*(o.j)); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,(ox+step*(o.i))+5,oy-step*(o.j)); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,(ox+step*(o.i))-5,oy-step*(o.j)); } else if(o.j==0&&abs(o.i)<30) { MoveToEx(hDC,ox,oy,NULL); LineTo(hDC,ox+step*(o.i),oy); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,(ox+step*(o.i)),oy-step*(o.j)+5); MoveToEx(hDC,ox+step*(o.i),oy-step*(o.j),NULL); LineTo(hDC,(ox+step*(o.i)),oy-step*(o.j)-5); } else if(abs(o.i)>30||abs(o.j)>30) { int xctr = ((o.i==0)?5:0); int yctr = ((o.j==0)?5:0); MoveToEx(hDC,ox,oy,NULL); LineTo(hDC,ox+(step*(o.i/5)),oy-(step*(o.j/5))); MoveToEx(hDC,ox+(step*(o.i/5)),oy-(step*(o.j/5)),NULL); LineTo(hDC,ox+(step*(o.i/5))+xctr,oy-(step*(o.j/5))+yctr); MoveToEx(hDC,ox+(step*(o.i/5)),oy-(step*(o.j/5)),NULL); LineTo(hDC,ox+(step*(o.i/5))-xctr,oy-(step*(o.j/5))-yctr); } } else { MoveToEx(hDC,ox,oy,NULL); LineTo(hDC,ox+step*(o.i/5),oy-step*(o.j/5)); MoveToEx(hDC,ox+step*(o.i/5),oy-step*(o.j/5),NULL); LineTo(hDC,ox+step*(o.i/5),(oy-step*(o.j/5))+(sy*5)); MoveToEx(hDC,ox+step*(o.i/5),oy-step*(o.j/5),NULL); LineTo(hDC,(ox+step*(o.i/5))-(sx*5),oy-step*(o.j/5)); hP = CreatePen(PS_DOT,1,RGB(0,0,0)); SelectObject(hDC,hP); MoveToEx(hDC,ox+step*(o.i/5),oy-step*(o.j/5),NULL); LineTo(hDC,ox,oy-step*(o.j/5)); MoveToEx(hDC,ox+step*(o.i/5),oy-step*(o.j/5),NULL); LineTo(hDC,ox+step*(o.i/5),oy); } DeleteObject(hP); DeleteObject(hF); EndPaint(hwnd, &ps); break; } case WM_DESTROY: { PostQuitMessage(0); break; } default: { return DefWindowProc(hwnd,message,wParam,lParam); break; } return 0; } } */ template class Vector { V i,j,k, mag,dx,dy,dz; static int VCtr; int ctr; public: V AX() {return dx;} /// Returns Angle With X-Axis V AY() {return dy;} /// Returns Angle With Y-Axis V AZ() {return dz;} /// Returns Angle With Z-Axis V l() {return cos(dx);} /// Returns Direction Cosine With X-Axis V m() {return cos(dy);} /// Returns Direction Cosine With Y-Axis V n() {return cos(dz);} /// Returns Direction Cosine With Z-Axis V Magnitude(){return mag;} /// Returns Magnitude Vector(): mag(0),dx(0),dy(0),dz(0),i(0),j(0),k(0) /// Default Constructor { VCtr++; ctr = VCtr; } Vector(V vi,V vj,V vk):i(vi),j(vj),k(vk) /// Initialize Vector using Direction Ratios { mag = sqrt(pow(i,2)+pow(j,2)+pow(k,2)); dx = acos((i/mag)); dy = acos((j/mag)); dz = acos((k/mag)); VCtr++; ctr = VCtr; } Vector(V mg,V x,V y,V z) /// Initialize Vector using Magnitude and Angles with Coordinate Axes { mag = mg; dx = degtorad(x); dy = degtorad(y); dz = degtorad(z); if(dx == Empty && dy!=Empty && dz!=Empty) { dx = acos(sqrt(1-pow(cos(dy),2)-pow(cos(dz),2))); } else if(dy == Empty && dx!=Empty && dz!=Empty) { dy = acos(sqrt(1-pow(cos(dx),2)-pow(cos(dz),2))); } else if(dz == Empty && dy!=Empty && dx!=Empty) { dz = acos(sqrt(1-pow(cos(dy),2)-pow(cos(dx),2))); } i = mag*l(); j = mag*m(); k = mag*n(); } Vector operator+(Vector b) /// Vector Addition { Vector res((i+b.i),(j+b.j),(k+b.k)); return res; } Vector operator-(Vector b) /// Vector Subtraction { Vector res((i-b.i),(j-b.j),(k-b.k)); return res; } V operator|(Vector b) /// Dot/Scalar Product { return ((i*b.i)+(j*b.j)*(k+b.k)); } Vector operator*(V var) /// Multiplication with a Scalar { Vector res((i*var),(j*var),(k*var)); return res; } Vector operator*(Vector b) /// Cross/Vector Product { Vector res((j*b.k-k*b.j),(k*b.i-i*b.k),(i*b.j-j*b.i)); return res; } Vector& operator+=(Vector b) /// Addition Assignment { *this = *this+b; return *this; } Vector& operator-=(Vector b) /// Subtraction Assignment { *this = *this-b; return *this; } Vector& operator*=(int k) /// Scalar Multiplication Assignment { *this = *this * k; return *this; } Vector& operator*=(Vector b) /// Vector Multiplication Assignment { *this = *this*b; return *this; } /* friend int Plot(Vector v) /// Function to Plot a Vector on a Small Scale Graph { HWND hwnd; MSG msgs; WNDCLASSEX Window; //if(typeid(V)==typeid(Complex)) //{ // vc.dxc=v.dx; vc.dyc=v.dy; vc.dzc=v.dz; // vc.magc=v.mag; vc.ic=v.i; // vc.jc=v.j; vc.kc=v.k; //} //else //{ vc.dx=v.dx; vc.dy=v.dy; vc.dz=v.dz; vc.mag=v.mag; vc.i=v.i; vc.j=v.j; vc.k=v.k; //} Window.hInstance = hInst; Window.lpszClassName = ("Vector "+to_string(v.ctr)).c_str(); Window.lpfnWndProc = WindowProcedure; Window.style = CS_VREDRAW | CS_HREDRAW; Window.cbSize = sizeof(WNDCLASSEX); Window.hIcon = LoadIcon (NULL, IDI_APPLICATION); Window.hIconSm = LoadIcon (NULL, IDI_APPLICATION); Window.hCursor = LoadCursor(NULL, IDC_ARROW); Window.lpszMenuName = NULL; Window.cbClsExtra = 0; Window.cbWndExtra = 0; Window.hbrBackground = CreateSolidBrush(RGB(255,255,255)); if(!RegisterClassEx(&Window)) { MessageBox(NULL,"Unable to Register the Class!","ERROR",MB_OK|MB_ICONERROR); return 0; } hwnd = CreateWindowEx( 0, ("Vector "+to_string(v.ctr)).c_str(), ("Vector "+to_string(v.ctr)).c_str(), WS_OVERLAPPED | WS_CAPTION | WS_SYSMENU, CW_USEDEFAULT, CW_USEDEFAULT, 400, 400, NULL, NULL, hInst, NULL); ShowWindow(hwnd,nCmdShow); while(GetMessage(&msgs,NULL,0,0)) { TranslateMessage(&msgs); DispatchMessage (&msgs); } return msgs.wParam; } */ friend ostream& operator<< (ostream& os, Vector v) { if(__indetail__||(os.flags()&detailed!=0)) { os.setf(ios::fixed); os.precision(2); os<<" V"<0?"(":" - (")<0?" + ":" - "); if(v.j!=0.0) os<<"("<0?" + ":" - "); if(v.k!=0.0) os<<"("<0?" + ":" - "); if(v.j!=0.0) os<<"("<0?" + ":" - "); if(v.k!=0.0) os<<"("<> (istream& is, Vector& v) { if(__read_using_ratio__) { is>>v.i>>v.j>>v.k; v.mag = sqrt(pow(v.i,2)+pow(v.j,2)+pow(v.k,2)); v.dx = acos((v.i/v.mag)); v.dy = acos((v.j/v.mag)); v.dz = acos((v.k/v.mag)); return is; } else { is>>v.mag>>v.dx>>v.dy>>v.dz; v.dx = degtorad(v.dx); v.dy = degtorad(v.dy); v.dz = degtorad(v.dz); v.i = v.mag*v.l(); v.j = v.mag*v.m(); v.k = v.mag*v.n(); return is; } } friend Vector operator*(int k, Vector v) /// Scalar Multiplication - friend { Vector res = v*k; return res; } }; template int Vector::VCtr = 0; /// Keeps a count of declared vectors. //int WINAPI WinMain(HINSTANCE hInstance,HINSTANCE hPrev,LPSTR CmdLine, int nCmdS) /// WinMain for use in Windows PC only... int main(int argc, char** argv) { try { //hInst = hInstance; nCmdShow = nCmdS; cout<(cout.flags())< v1(1,2,3), v2(2,1,1); /// Declares the Initial Vectors. /// There may be some problem on sololearn due to use of static members, though it works fine on a PC. cout.setf(detailed); /// Set to detail formatting. cout<(cout.flags())<>v1;cout<>v2;cout<(cout.flags())<(cout.flags())<