~njansson/dolfin/hpc

Viewing changes to src/kernel/la/NewVector.cpp

Committer: fenics
Date: 2005-02-23 03:46:01 UTC
Revision ID: devnull@localhost-20050223034601-az3dhrq0nqas5xza

Tailorized "2005-02-22 21:46:01 by fenics"
Remaining changes, everything didn't go through before.

Also changed so that PETSc 2.2.1 is now the default version.

files added:
src/demo/solvers/ode/homotopy/simple/plotsolution.m

src/kernel/ode/Homotopy.cpp

src/kernel/ode/HomotopyJacobian.cpp

src/kernel/ode/HomotopyODE.cpp

src/kernel/ode/dolfin/Homotopy.h

src/kernel/ode/dolfin/HomotopyJacobian.h

src/kernel/ode/dolfin/HomotopyODE.h

files modified:
ChangeLog

configure

configure.ac

src/config/Makefile.am

src/config/Makefile.in

src/demo/la/Makefile

src/demo/la/newla/main.cpp

src/demo/solvers/ode/homotopy/economy/main.cpp

src/demo/solvers/ode/homotopy/simple/main.cpp

src/demo/solvers/ode/stiff/Makefile

src/demo/test/main.cpp

src/kernel/fem/Makefile

src/kernel/la/NewGMRES.cpp

src/kernel/la/NewVector.cpp

src/kernel/la/PETScManager.cpp

src/kernel/la/VirtualMatrix.cpp

src/kernel/la/dolfin/NewGMRES.h

src/kernel/la/dolfin/NewMatrix.h

src/kernel/la/dolfin/NewPreconditioner.h

src/kernel/la/dolfin/NewVector.h

src/kernel/la/dolfin/VirtualMatrix.h

src/kernel/la/dolfin/pcimpl.h

src/kernel/ode/ComplexODE.cpp

src/kernel/ode/Makefile.am

src/kernel/ode/Makefile.in

src/kernel/ode/MonoAdaptiveNewtonSolver.cpp

src/kernel/ode/MonoAdaptiveTimeSlab.cpp

src/kernel/ode/MultiAdaptiveTimeSlab.cpp

src/kernel/ode/NewTimeStepper.cpp

src/kernel/ode/ODE.cpp

src/kernel/ode/dolfin/ComplexODE.h

src/kernel/ode/dolfin/Makefile.am

src/kernel/ode/dolfin/Makefile.in

src/kernel/ode/dolfin/MonoAdaptiveNewtonSolver.h

src/kernel/ode/dolfin/MonoAdaptiveTimeSlab.h

src/kernel/ode/dolfin/MultiAdaptiveTimeSlab.h

src/kernel/ode/dolfin/NewTimeSlab.h

src/kernel/ode/dolfin/NewTimeStepper.h

src/kernel/ode/dolfin/ODE.h

src/kernel/ode/dolfin/dolfin_ode.h

Show diffs side-by-side

added added

removed removed

src/kernel/la/NewVector.cpp

190

return *this;

191

}

192

//-----------------------------------------------------------------------------

193

const NewVector& NewVector::operator+= (const NewVector& x)

194

{

195

const real a = 1.0;

196

VecAXPY(&a, x.x, this->x);

197

198

return *this;

199

}

200

//-----------------------------------------------------------------------------

201

const NewVector& NewVector::operator-= (const NewVector& x)

202

{

203

const real a = -1.0;

204

VecAXPY(&a, x.x, this->x);

205

206

return *this;

207

}

208

//-----------------------------------------------------------------------------

209

const NewVector& NewVector::operator*= (real a)

210

{

211

VecScale(&a, x);

212

213

return *this;

214

}

215

//-----------------------------------------------------------------------------

216

const NewVector& NewVector::operator/= (real a)

217

{

218

dolfin_assert(a != 0.0);

219

const real b = 1.0 / a;

220

VecScale(&b, x);

221

222

return *this;

223

}

224

//-----------------------------------------------------------------------------

225

real NewVector::norm(NormType type) const

226

{

227

real value = 0.0;

228

229

switch (type) {

230

case l1:

231

VecNorm(x, NORM_1, &value);

232

break;

233

case l2:

234

VecNorm(x, NORM_2, &value);

235

break;

236

default:

237

VecNorm(x, NORM_INFINITY, &value);

238

}

239

240

return value;

241

}

242

//-----------------------------------------------------------------------------

193

243

void NewVector::disp() const

194

244

{

195

245

// FIXME: Maybe this could be an option?

Older »