Category

Development Platform
Home > SourceCode/Document > Linux-Unix program > View Code
matrix_point.cpp
Package: nurbs++-3.0.11.zip　[view]
Upload User: ysz425
Upload Date: 2008-03-14
Package Size: 675k
Code Size: 6k
Category:

Linux-Unix program

Development Platform:

Unix_Linux

matrix_point.cpp：Code Content
  1. /*=============================================================================
  2.         File: matrix.cpp
  3.      Purpose:       
  4.     Revision: $Id: matrix_point.cpp,v 1.2 2002/05/13 21:07:45 philosophil Exp $
  5.   Created by: Philippe Lavoie          (3 Oct, 1996)
  6.  Modified by: 
  8.  Copyright notice:
  9.           Copyright (C) 1996-1998 Philippe Lavoie
  10.  
  11.   This library is free software; you can redistribute it and/or
  12.   modify it under the terms of the GNU Library General Public
  13.   License as published by the Free Software Foundation; either
  14.   version 2 of the License, or (at your option) any later version.
  15.  
  16.   This library is distributed in the hope that it will be useful,
  17.   but WITHOUT ANY WARRANTY; without even the implied warranty of
  18.   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19.   Library General Public License for more details.
  20.  
  21.   You should have received a copy of the GNU Library General Public
  22.   License along with this library; if not, write to the Free
  23.   Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  24. =============================================================================*/
  26. #include "matrix.cpp"
  28. namespace PLib {
  30.   double
  31.     Matrix<Point3Df>::norm(void) {
  32.     int i,j ;
  33.     double sumX, sumY, sumZ, maxsum;
  34.     int init=0 ;
  35.     Point3Df *ptr ;
  36.     ptr = m-1 ;
  37.     maxsum = -1 ; // shuts up the warning messages 
  38.     for(i=0;i<rows();++i){
  39.       sumX = 0.0 ;
  40.       sumY = 0.0 ;
  41.       sumZ = 0.0 ;
  42.       for ( j = 0; j < cols(); ++j) {
  43. sumX += (*ptr).x() * (*ptr).x() ;
  44. sumY += (*ptr).y() * (*ptr).y() ;
  45. sumZ += (*ptr).z() * (*ptr).z() ;
  46.       }
  47.       if(init)
  48. maxsum = (maxsum>(sumX+sumY+sumZ)) ? maxsum : (sumX+sumY+sumZ);
  49.       else{
  50. maxsum = (sumX+sumY+sumZ) ;
  51. init = 1;
  52.       }
  53.       ++ptr ;
  54.     }
  55.     return sqrt(maxsum);
  56.   }
  57.   
  58.   double
  59.     Matrix<Point3Dd>::norm(void) {
  60.     int i,j ;
  61.     double sumX, sumY, sumZ, maxsum;
  62.     int init=0 ;
  63.     Point3Dd *ptr ;
  64.     ptr = m-1 ;
  65.     maxsum = -1 ; // shuts up the warning messages 
  66.     for(i=0;i<rows();++i){
  67.       sumX = 0.0 ;
  68.       sumY = 0.0 ;
  69.       sumZ = 0.0 ;
  70.       for ( j = 0; j < cols(); ++j) {
  71. sumX += (*ptr).x() * (*ptr).x() ;
  72. sumY += (*ptr).y() * (*ptr).y() ;
  73. sumZ += (*ptr).z() * (*ptr).z() ;
  74.       }
  75.       if(init)
  76. maxsum = (maxsum>(sumX+sumY+sumZ)) ? maxsum : (sumX+sumY+sumZ);
  77.       else{
  78. maxsum = (sumX+sumY+sumZ) ;
  79. init = 1;
  80.       }
  81.       ++ptr ;
  82.     }
  83.     return sqrt(maxsum);
  84.   }
  85.   
  86.   double
  87.     Matrix<Point2Df>::norm(void) {
  88.     int i,j ;
  89.     double sumX, sumY, sumZ, maxsum;
  90.     int init=0 ;
  91.     Point2Df *ptr ;
  92.     ptr = m-1 ;
  93.     maxsum = -1 ; // shuts up the warning messages 
  94.     for(i=0;i<rows();++i){
  95.       sumX = 0.0 ;
  96.       sumY = 0.0 ;
  97.       sumZ = 0.0 ;
  98.       for ( j = 0; j < cols(); ++j) {
  99. sumX += (*ptr).x() * (*ptr).x() ;
  100. sumY += (*ptr).y() * (*ptr).y() ;
  101. sumZ += (*ptr).z() * (*ptr).z() ;
  102.       }
  103.       if(init)
  104. maxsum = (maxsum>(sumX+sumY+sumZ)) ? maxsum : (sumX+sumY+sumZ);
  105.       else{
  106. maxsum = (sumX+sumY+sumZ) ;
  107. init = 1;
  108.       }
  109.       ++ptr ;
  110.     }
  111.     return sqrt(maxsum);
  112.   }
  114.   double
  115.     Matrix<Point2Dd>::norm(void) {
  116.     int i,j ;
  117.     double sumX, sumY, sumZ, maxsum;
  118.     int init=0 ;
  119.     Point2Dd *ptr ;
  120.     ptr = m-1 ;
  121.     maxsum = -1 ; // shuts up the warning messages 
  122.     for(i=0;i<rows();++i){
  123.       sumX = 0.0 ;
  124.       sumY = 0.0 ;
  125.       sumZ = 0.0 ;
  126.       for ( j = 0; j < cols(); ++j) {
  127. sumX += (*ptr).x() * (*ptr).x() ;
  128. sumY += (*ptr).y() * (*ptr).y() ;
  129. sumZ += (*ptr).z() * (*ptr).z() ;
  130.       }
  131.       if(init)
  132. maxsum = (maxsum>(sumX+sumY+sumZ)) ? maxsum : (sumX+sumY+sumZ);
  133.       else{
  134. maxsum = (sumX+sumY+sumZ) ;
  135. init = 1;
  136.       }
  137.       ++ptr ;
  138.     }
  139.     return sqrt(maxsum);
  140.   }
  141.   
  143. #ifdef NO_IMPLICIT_TEMPLATES
  145.   // Point3D instantiation
  146.   
  147.   template class Matrix<Point3Df> ;
  148.   
  149.   template Matrix<Point3Df> operator+(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  150.   template Matrix<Point3Df> operator-(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  151.   template Matrix<Point3Df> operator*(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  152.   template Matrix<Point3Df> operator*(const double,const Matrix<Point3Df>&);
  153.   template int operator==(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  154.   // template int operator!=(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  155.   template Matrix<Point3Df> comm(const Matrix<Point3Df>&,const Matrix<Point3Df>&);
  156.   
  157.   
  158.   template class Matrix<Point3Dd> ;
  159.   
  160.   template Matrix<Point3Dd> operator+(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  161.   template Matrix<Point3Dd> operator-(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  162.   template Matrix<Point3Dd> operator*(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  163.   template Matrix<Point3Dd> operator*(const double,const Matrix<Point3Dd>&);
  164.   template int operator==(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  165.   //template int operator!=(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  166.   template Matrix<Point3Dd> comm(const Matrix<Point3Dd>&,const Matrix<Point3Dd>&);
  167.   
  168.   // Point2D instantiation
  169.   
  170.   template class Matrix<Point2Df> ;
  171.   
  172.   template Matrix<Point2Df> operator+(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  173.   template Matrix<Point2Df> operator-(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  174.   template Matrix<Point2Df> operator*(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  175.   template Matrix<Point2Df> operator*(const double,const Matrix<Point2Df>&);
  176.   template int operator==(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  177.   //template int operator!=(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  178.   template Matrix<Point2Df> comm(const Matrix<Point2Df>&,const Matrix<Point2Df>&);
  179.   
  181.   template class Matrix<Point2Dd> ;
  182.   
  183.   template Matrix<Point2Dd> operator+(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  184.   template Matrix<Point2Dd> operator-(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  185.   template Matrix<Point2Dd> operator*(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  186.   template Matrix<Point2Dd> operator*(const double,const Matrix<Point2Dd>&);
  187.   template int operator==(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  188.   //template int operator!=(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  189.   template Matrix<Point2Dd> comm(const Matrix<Point2Dd>&,const Matrix<Point2Dd>&);
  191. #endif
  193. }