0

作为我正在编写的程序的一部分,我需要找到一个双精度向量和一个复数双精度向量的叉积。我写了一个我觉得应该这样做的函数,但是当我调用它时,我收到以下错误:

error: no matching function for call to ‘CrossProduct1D(std::vector< double, std::allocator<double> >&, std::vector<std::complex<double>, std::allocator<std::complex<double> > >&)’

这是我的代码:

#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <vector>
#include <math.h>
#include <complex>
using namespace std;

//1D cross product
template <typename T>
vector<T> CrossProduct1D(vector<T> const &a, vector<T> const &b)
{
  vector<T> r (a.size());  
  r[0] = a[1]*b[2]-a[2]*b[1];
  r[1] = a[2]*b[0]-a[0]*b[2];
  r[2] = a[0]*b[1]-a[1]*b[0];
  return r;
}

//::::::::importing data from text::::::::::
  vector<string> ezivec;
  ezivec.reserve(4000);

  string ezidat("ez.i.txt");

  ifstream ezifile;
  ezifile.open(ezidat.c_str());

  if(!ezifile.is_open())
    {
      cerr<<"Error opening file : "<<ezidat.c_str()<<endl;
      return -1;
    }

  string ezistr; //store lines in vector
  while(getline(ezifile, ezistr, ';'))
    {
      ezivec.push_back(ezistr);
    }

  ezifile.close();

 //Converting from vector of strings to vector of floats
  vector<double> ezi (ezivec.size());
  for(int i = 0; i < ezivec.size(); ++i)
    {
      ezi[i] = string_to_T<double>(ezivec[i]);
    }

 //

  vector<string> ezrvec;
  ezrvec.reserve(4000); 

  string ezrdat("ez.r.txt");

  ifstream ezrfile;
  ezrfile.open(ezrdat.c_str());

  if(!ezrfile.is_open())
    {
      cerr<<"Error opening file : "<<ezrdat.c_str()<<endl;
      return -1;
    }

  string ezrstr;
  while(getline(ezrfile, ezrstr, ';'))
    {
      ezrvec.push_back(ezrstr);
    }

  ezrfile.close();

  vector<double> ezr (ezrvec.size());
  for(int i = 0; i < ezrvec.size(); ++i)
    {
      ezr[i] = string_to_T<double>(ezrvec[i]);
    }

//:::::::defining vectors:::::::
vector<vector<complex<double> > > E0 (ezi.size(), vector<complex<double> > (3));
for(int i = 0; i < ezi.size(); i++)
{
  E0[i][0].real() = 0.0;
  E0[i][0].imag() = 0.0;
  E0[i][1].real() = 0.0;
  E0[i][1].imag() = 0.0;;
  E0[i][2].real() = ezr[i];
  E0[i][2].imag() = ezi[i];
} 

vector<double> n_a (3);
n_a[0] = 1.0;
n_a[1] = 0.0;
n_a[2] = 0.0;

//:::::::calling cross product:::::::
for(int j = 1; j < jmax; j++)
    {
      M[j] = CrossProduct1D(n_a, E0[j]);
    }

“ez.i.txt”和“ez.r.txt”是由分号分隔的 4000 个数字的文本文件。

4

1 回答 1

10

您的模板函数在单个类型上参数化,T, 并接受两个vector<T>,但您试图将两种不同类型的向量传递给它,因此没有一个T可以选择。

您可以有两个模板参数,例如

template<class T, class U> CrossProduct1D(std::vector<T> const& a,
                                          std::vector<U> const& b)
于 2012-07-19T17:28:59.973 回答