Save od ofstream stifness matrix asn file.txt

frederichecht · August 25, 2025, 9:30am

remark, i can save sparse matrix see

examples/plugin/SaveHB_sample.edp

develop

// -*- FreeFem++ -*-
// Sample of SaveHB
// Author: fujiwara@acs.i.kyoto-u.ac.jp
mesh Th = square(3,3);
macro Grad(u) [dx(u),dy(u)] // EOM

fespace Vh(Th,P1);
varf vP(u,v) = int2d(Th)(Grad(u)'*Grad(v)) + int2d(Th)(v) + on(1,u=1);

matrix A = vP(Vh,Vh);
real[int] rhs = vP(0,Vh);

load "SaveHB" 
int ret = SaveHB("sample.hb", A, rhs, "sample of SaveHB");
A.CSR; 
cout << A << endl; 
cout << " A = " << A << endl;
cout << " rhs  = " << rhs << endl;

assert(ret == 0);

This file has been truncated. show original

Topic		Replies	Views
About stiffness martix General Discussion	1	323	September 24, 2022
Export matrix to matlab General Discussion	2	1389	March 9, 2026
Can freefem + + find the largest element of the stiffness matrix? General Discussion	6	662	May 16, 2022
Get global stiffness matrix from freefem General Discussion	2	1096	August 22, 2022
Read varlues created from FreeFem ofstream General Discussion	1	451	July 23, 2022

Save od ofstream stifness matrix asn file.txt

Related topics