🏆本文收录于《CSDN问答解惑-专业版》专栏，主要记录项目实战过程中的Bug之前因后果及提供真实有效的解决方案，希望能够助你一臂之力，帮你早日登顶实现财富自由🚀；同时，欢迎大家关注&&收藏&&订阅！持续更新中，up！up！up！！
问题描述
c# 直接使用c++ 类库文件。需要完整的改写 以下 c++ 代码和 c#完整调用例子
每一个函数都要有，不要省略代码。
class NendToCSharp
{
private:
int N; //复系数多项式次数
Complex *c; //复系数多项式存储空间首地址
public:
NendToCSharp(void);
~NendToCSharp(void);
NendToCSharp(int nn=0, Complex *pp=NULL) //构造函数
{ N = nn; c = pp; }
template <class T>
int gauss(T *a, T *b, int n)
{
int *js,k,i,j,is,p,q;
double d,t;
T s;
js=new int[n];
for (k=0;k<=n-2;k++) //消元过程
{
d=0.0; //全选主元
for (i=k;i<=n-1;i++)
for (j=k;j<=n-1;j++)
{
t=ffabs(a[i*n+j]);
if (t>d) { d=t; js[k]=j; is=i;}
}
if (d+1.0==1.0) //系数矩阵奇异，求解失败！
{
for (i=0; i<n; i++) b[i] = init(s);
cout <<"系数矩阵奇异，求解失败！\n";
delete[]js;
return 0;
}
if (js[k]!=k) //列交换
{
for (i=0;i<=n-1;i++)
{
p=i*n+k; q=i*n+js[k];
s=a[p]; a[p]=a[q]; a[q]=s;
}
}
if (is!=k) //行交换
{
for (j=k;j<=n-1;j++)
{
p=k*n+j; q=is*n+j;
s=a[p]; a[p]=a[q]; a[q]=s;
}
s=b[k]; b[k]=b[is]; b[is]=s;
}
s=a[k*n+k];
for (j=k+1;j<=n-1;j++) //归一化
{
p=k*n+j; a[p]=a[p]/s;
}
b[k]=b[k]/s;
for (i=k+1;i<=n-1;i++) //消元
{
for (j=k+1;j<=n-1;j++)
{
p=i*n+j;
a[p]=a[p]-a[i*n+k]*a[k*n+j];
}
b[i]=b[i]-a[i*n+k]*b[k];
}
}
s=a[(n-1)*n+n-1];
if (ffabs(s)+1.0==1.0) //系数矩阵奇异，求解失败！
{
for (i=0; i<n; i++) b[i] = init(s);
cout <<"系数矩阵奇异，求解失败！\n";
delete[] js;
return 0;
}
b[n-1]=b[n-1]/s; //回代过程
for (i=n-2;i>=0;i--)
{
s=init(s);
for (j=i+1;j<=n-1;j++) s=s+a[i*n+j]*b[j];
b[i]=b[i]-s;
}
js[n-1]=n-1;
for (k=n-1;k>=0;k--) //恢复
if (js[k]!=k)
{
s=b[k]; b[k]=b[js[k]]; b[js[k]]=s;
}
delete[] js;
return 1;
}
template <class T> //模板声明T为类型参数
int gauss_jordan(T *a, T *b, int n)
{
int *js,k,i,j,is,p,q;
double d,t;
T s;
js=new int[n];
for (k=0;k<=n-1;k++) //消去过程
{
d=0.0; //全选主元
for (i=k;i<=n-1;i++)
for (j=k;j<=n-1;j++)
{
t=ffabs(a[i*n+j]);
if (t>d) { d=t; js[k]=j; is=i;}
}
if (d+1.0==1.0) //系数矩阵奇异，求解失败！
{
cout <<"系数矩阵奇异，求解失败！\n";
for (i=0; i<n; i++) b[i]=init(s);
delete[]js;
return 0;
}
if (js[k]!=k) //列交换
{
for (i=0;i<=n-1;i++)
{
p=i*n+k; q=i*n+js[k];
s=a[p]; a[p]=a[q]; a[q]=s;
}
}
if (is!=k) //行交换
{
for (j=k;j<=n-1;j++)
{
p=k*n+j; q=is*n+j;
s=a[p]; a[p]=a[q]; a[q]=s;
}
s=b[k]; b[k]=b[is]; b[is]=s;
}
s=a[k*n+k];
for (j=k+1;j<=n-1;j++) //归一化
{
p=k*n+j; a[p]=a[p]/s;
}
b[k]=b[k]/s;
for (i=0;i<=n-1;i++) //消元
{
if (i!=k)
{
for (j=k+1;j<=n-1;j++)
{
p=i*n+j;
a[p]=a[p]-a[i*n+k]*a[k*n+j];
}
b[i]=b[i]-a[i*n+k]*b[k];
}
}
}
for (k=n-1;k>=0;k--) //恢复
if (js[k]!=k)
{
s=b[k]; b[k]=b[js[k]]; b[js[k]]=s;
}
delete[] js;
return 1;
}
int trde(double b[], int n, int m, double d[])
{
int k,j;
double s;
if (m!=(3*n-2)) return -2;
for (k=0;k<=n-2;k++)
{
j=3*k; s=b[j];
if (fabs(s)+1.0==1.0) return 0;
b[j+1]=b[j+1]/s;
d[k]=d[k]/s;
b[j+3]=b[j+3]-b[j+2]*b[j+1];
d[k+1]=d[k+1]-b[j+2]*d[k];
}
s=b[3*n-3];
if (fabs(s)+1.0==1.0) return 0;
d[n-1]=d[n-1]/s;
for (k=n-2;k>=0;k--) d[k]=d[k]-b[3*k+1]*d[k+1];
return(2);
}
int band(double b[], double d[], int n, int l, int il, int m)
{
int ls,k,i,j,is,u,v;
double p,t;
if (il!=(2*l+1)) return(-2);
ls=l;
for (k=0;k<=n-2;k++)
{
p=0.0;
for (i=k;i<=ls;i++)
{
t=fabs(b[i*il]);
if (t>p) {p=t; is=i;}
}
if (p+1.0==1.0) return(0);
for (j=0;j<=m-1;j++)
{
u=k*m+j; v=is*m+j;
t=d[u]; d[u]=d[v]; d[v]=t;
}
for (j=0;j<=il-1;j++)
{
u=k*il+j; v=is*il+j;
t=b[u]; b[u]=b[v]; b[v]=t;
}
for (j=0;j<=m-1;j++)
{
u=k*m+j; d[u]=d[u]/b[k*il];
}
for (j=1;j<=il-1;j++)
{
u=k*il+j; b[u]=b[u]/b[k*il];
}
for (i=k+1;i<=ls;i++)
{
t=b[i*il];
for (j=0;j<=m-1;j++)
{
u=i*m+j; v=k*m+j; d[u]=d[u]-t*d[v];
}
for (j=1;j<=il-1;j++)
{
u=i*il+j; v=k*il+j; b[u-1]=b[u]-t*b[v];
}
u=i*il+il-1; b[u]=0.0;
}
if (ls!=(n-1)) ls=ls+1;
}
p=b[(n-1)*il];
if (fabs(p)+1.0==1.0) return(0);
for (j=0;j<=m-1;j++)
{
u=(n-1)*m+j; d[u]=d[u]/p;
}
ls=1;
for (i=n-2;i>=0;i--)
{
for (k=0;k<=m-1;k++)
{
u=i*m+k;
for (j=1;j<=ls;j++)
{
v=i*il+j; is=(i+j)*m+k;
d[u]=d[u]-b[v]*d[is];
}
}
if (ls!=(il-1)) ls=ls+1;
}
return(2);
}
int ldle(double a[], int n, int m, double c[])
{
int i,j,l,k,u,v,w,k1,k2,k3;
double p;
//省略
return(2);
}
int chlk(double a[], int n, int m, double d[])
{
int i,j,k,u,v;
if ((a[0]+1.0==1.0)||(a[0]<0.0)) return(0);
a[0]=sqrt(a[0]);
//省略
return(2);
}
int tlvs(double t[], int n, double b[], double x[])
{
int i,j,k;
double a,beta,q,c,h,*y,*s;
s=new double[n];
y=new double[n];
a=t[0];
//省略
return(1);
}
int seidel(double *a, double *b, int n, double *x, double eps)
{
int i,j,u,v;
double p,t,s,q;
//省略
return 1;
}
void grad(double a[],int n,double b[],double eps,double x[])
{
Matrix mx;
int i,k;
double *p,*r,*s,*q,alpha,beta,d,e;
p=new double[n];
r=new double[n];
s=new double[n];
q=new double[n];
for (i=0; i<=n-1; i++)
{ x[i]=0.0; p[i]=b[i]; r[i]=b[i]; }
i=0;
//省略
delete[] p; delete[] r; delete[] s; delete[] q;
return;
}
int maqr(double a[],int m,int n,double q[])
{
int i,j,k,l,nn,p,jj;
double u,alpha,w,t;
//省略
return 1;
}
int gmqr(double a[], int m, int n, double b[], double q[])
{
int i,j;
double d,*c;
c=new double[n];
i=maqr(a,m,n,q);
if (i==0)
{
delete[] c; return(0);
}
for (i=0; i<=n-1; i++)
{
d=0.0;
for (j=0; j<=m-1; j++) d=d+q[j*m+i]*b[j];
c[i]=d;
}
b[n-1]=c[n-1]/a[n*n-1];
for (i=n-2; i>=0; i--)
{
d=0.0;
for (j=i+1; j<=n-1; j++) d=d+a[i*n+j]*b[j];
b[i]=(c[i]-d)/a[i*n+i];
}
delete[] c; return(1);
}
int gmiv(double a[], int m, int n, double b[], double x[],
double aa[], double eps, double u[], double v[], int ka)
{
//省略
return(1);
}
int bingt(double a[], int n, double b[], double eps, double x[])
{
int i,j,k;
double q, qq;
double *p,*r,*e;
p=new double[n*n];
r=new double[n];
e=new double[n];
k=0;
//省略
delete[] p; delete[] r; delete[] e; return 1;
}
};
如上问题有来自我自身项目开发，有的收集网站，有的来自读者…如有侵权，立马删除。
解决方案
如下是上述问题的解决方案，仅供参考：
将C++类库文件直接用于C#项目涉及到互操作性问题。C++和C#是两种不同的编程语言，它们使用不同的运行时环境和内存管理方式。为了在C#中使用C++编写的代码，您需要使用特定的互操作技术。以下是一些步骤和示例，帮助您实现这一目标：
1. 创建C++/CLI包装器
C++/CLI（Common Language Infrastructure）是C++的一个扩展，允许您编写可以与.NET语言（包括C#）互操作的代码。您可以创建一个C++/CLI包装器，将您的C++代码暴露给C#。
C++/CLI示例
假设您有一个C++类NendToCSharp，您可以这样创建包装器：
// NendToCSharpWrapper.cpp
#include "NendToCSharp.h"
using namespace System;
public ref class NendToCSharpWrapper
{
private:
NendToCSharp* nativeObject;
public:
NendToCSharpWrapper() : nativeObject(new NendToCSharp()) {}
~NendToCSharpWrapper()
{
delete nativeObject;
}
// 将其他NendToCSharp的方法暴露给C#
// 例如：
int Gauss(array<double>^ a, array<double>^ b)
{
pin_ptr<double> aPin = &a[0];
pin_ptr<double> bPin = &b[0];
double* aPtr = aPin;
double* bPtr = bPin;
return nativeObject->gauss(aPtr, bPtr, a->Length);
}
// 其他方法...
};
2. 在C#中调用C++/CLI包装器
一旦您有了C++/CLI包装器，就可以在C#项目中添加对该包装器的引用，并像调用任何其他C#类一样调用它。
C#调用示例
using System;
class Program
{
static void Main(string[] args)
{
// 创建C++/CLI包装器的实例
NendToCSharpWrapper wrapper = new NendToCSharpWrapper();
// 创建数组并初始化（示例数据）
double[] a = new double[] { /* ... 数据 ... */ };
double[] b = new double[a.Length];
// 调用包装器方法
int result = wrapper.Gauss(a, b);
Console.WriteLine("Gauss method returned: " + result);
// 处理结果...
}
}
注意事项
- 您需要将C++代码中的原始指针和数组转换为C#中的
array或List<T>。 - 处理C++和C#之间的数据类型转换，特别是对于复杂的数据结构。
- 确保C++/CLI项目配置正确，以便它可以与C#项目一起编译和运行。
替代方案
- 使用P/Invoke：对于简单的函数调用，您可以使用P/Invoke从C#调用C++ DLL中的函数。但这通常适用于较小的、不涉及复杂数据结构的函数。
- 使用COM：如果您的C++代码可以通过COM暴露，您可以在C#中使用
System.Runtime.InteropServices来调用它。
由于您的C++代码相当复杂，并且涉及到模板和复杂的数据结构，C++/CLI可能是最合适的解决方案。请注意，这只是一个高层次的概述，具体实现可能需要根据您的具体需求进行调整。
希望如上措施及解决方案能够帮到有需要的你。
PS：如若遇到采纳如下方案还是未解决的同学，希望不要抱怨&&急躁，毕竟影响因素众多，我写出来也是希望能够尽最大努力帮助到同类似问题的小伙伴，即把你未解决或者产生新Bug黏贴在评论区，我们大家一起来努力，一起帮你看看，可以不咯。
若有对当前Bug有与如下提供的方法不一致，有个不情之请，希望你能把你的新思路或新方法分享到评论区，一起学习，目的就是帮助更多所需要的同学，正所谓「赠人玫瑰，手留余香」。
☀️写在最后
ok，以上就是我这期的Bug修复内容啦，如果还想查找更多解决方案，你可以看看我专门收集Bug及提供解决方案的专栏《CSDN问答解惑-专业版》，都是实战中碰到的Bug，希望对你有所帮助。到此，咱们下期拜拜。
码字不易，如果这篇文章对你有所帮助，帮忙给 bug菌 来个一键三连(关注、点赞、收藏) ，您的支持就是我坚持写作分享知识点传播技术的最大动力。
同时也推荐大家关注我的硬核公众号:「猿圈奇妙屋」 ；以第一手学习bug菌的首发干货，不仅能学习更多技术硬货，还可白嫖最新BAT大厂面试真题、4000G Pdf技术书籍、万份简历/PPT模板、技术文章Markdown文档等海量资料，你想要的我都有！
📣关于我
我是bug菌，CSDN | 掘金 | InfoQ | 51CTO | 华为云 | 阿里云 | 腾讯云 等社区博客专家，C站博客之星Top30，华为云2023年度十佳博主，掘金多年度人气作者Top40，掘金等各大社区平台签约作者，51CTO年度博主Top12，掘金/InfoQ/51CTO等社区优质创作者；全网粉丝合计 30w+；硬核微信公众号「猿圈奇妙屋」，欢迎你的加入！免费白嫖最新BAT互联网公司面试真题、4000G PDF电子书籍、简历模板等海量资料，你想要的我都有，关键是你不来拿哇。
