vector.cpp 2.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122
  1. //#include "stdafx.h"
  2. #include <stdlib.h>
  3. #include <math.h>
  4. #include "vector.h"
  5. /*************************** constructors **************************/
  6. Vector::Vector(int n) {
  7. if (n < 0) {
  8. cerr << "Vector: negative dimension." << endl;
  9. exit(1);
  10. }
  11. d = n;
  12. if (d > 0) {
  13. v = new double[d];
  14. while (n--) v[n] = 0.0;
  15. }
  16. else v = NULL;
  17. }
  18. Vector::Vector(int n, double v0) {
  19. if (n <= 0) {
  20. cerr << "Vector: undefined dimension." << endl;
  21. exit(1);
  22. }
  23. d = n;
  24. v = new double[d];
  25. while (n--) v[n] = v0;
  26. }
  27. Vector::Vector(int n, double* v0) {
  28. if (n <= 0) {
  29. cerr << "Vector: undefined dimension." << endl;
  30. exit(1);
  31. }
  32. d = n;
  33. v = new double[d];
  34. while (n--) v[n] = *(v0+n);
  35. }
  36. Vector::Vector(const Vector& v0) {
  37. d = v0.d;
  38. if (d > 0) {
  39. v = new double[d];
  40. for(int i = 0; i < d; i++) v[i] = v0.v[i];
  41. }
  42. else v = NULL;
  43. }
  44. /************************* members ****************************/
  45. /************** friends *********************/
  46. ostream& operator<<(ostream& out, const Vector& v) {
  47. for (int i = 0; i < v.d; i++)
  48. //out << v.v[i] << ' ' << endl;
  49. out << " [" << i << "] = " << v.v[i] << ' ' << endl;
  50. return (out);
  51. }
  52. istream& operator>>(istream& in, Vector& x) {
  53. int i = 0;
  54. while ((i < x.d) && (in >> x.v[i++]));
  55. return (in);
  56. }
  57. Vector operator*(double f, const Vector& v) {
  58. return (v*f);
  59. }
  60. double max(const Vector& vec) {
  61. int n = vec.dim();
  62. double c = vec[n-1];
  63. while (n--) c = (c > vec[n] ? c : vec[n]);
  64. return c;
  65. }
  66. double min(const Vector& vec) {
  67. int n = vec.dim();
  68. double c = vec[n-1];
  69. while (n--) c = (c < vec[n] ? c : vec[n]);
  70. return c;
  71. }
  72. /************** nonfriends *********************/
  73. double norm(const Vector& v) {
  74. double result = 0;
  75. int n = v.dim();
  76. while (n--) result += v[n];
  77. return (result);
  78. }
  79. double norm2(const Vector& v) {
  80. double result = 0;
  81. int n = v.dim();
  82. while (n--) result += v[n]*v[n];
  83. return (sqrt(result));
  84. }
  85. Vector operator^(const Vector& v,const Vector& w) {
  86. if (v.dim() != w.dim()) {
  87. cerr << "Vector=: vectors have different length" << endl;
  88. exit(1);
  89. }
  90. int n = v.dim();
  91. Vector result(n, 0.0);
  92. while (n--) result[n] = v[n] * w[n];
  93. return (result);
  94. }