dvect6.cpp

Go to the documentation of this file.

/*
 * $Id$
 *
 * Author: David Fournier
 * Copyright (c) 2008-2012 Regents of the University of California
 */
#include "fvar.hpp"

#ifdef __TURBOC__
  #pragma hdrstop
#endif

#include <math.h>

void shape_check(const dvector& v1, const dvector& v2,
  const char *function_name);

dvector sin(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=sin(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector cos(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=cos(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector tan(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=tan(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector sinh(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=sinh(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector cosh(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=cosh(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector tanh(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=tanh(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector pow(const dvector& v1, const double x)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=pow(v1.elem(i),x);             // ***
   }
   return(tmp);
 }

dvector pow(const dvector& v1, int x)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=pow(v1.elem(i),x);             // ***
   }
   return(tmp);
 }

dvector pow(const dvector& v1, const dvector& v2)// ***
 {
   shape_check(v1,v2,"dvector pow(const dvector& v1,constdvector& v2)");

   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=pow(v1.elem(i),v2.elem(i));             // ***
   }
   return(tmp);
 }

dvector pow(const double x, const dvector& v2)// ***
 {
   dvector tmp(v2.indexmin(),v2.indexmax());
   for (int i=v2.indexmin();i<=v2.indexmax();i++)
   {
     tmp.elem(i)=pow(x,v2.elem(i));             // ***
   }
   return(tmp);
 }

dvector asin(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=asin(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector acos(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=acos(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector atan(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=atan(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector log10(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=log10(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector sqrt(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=sqrt(v1.elem(i));             // ***
   }
   return(tmp);
 }

dvector sqr(const dvector& v1)// ***
 {
   dvector tmp(v1.indexmin(),v1.indexmax());
   for (int i=v1.indexmin();i<=v1.indexmax();i++)
   {
     tmp.elem(i)=sqrt(v1.elem(i));             // ***
   }
   return(tmp);
 }