/* -------------------------------------------------------------------------- * * OpenMM * * -------------------------------------------------------------------------- * * This is part of the OpenMM molecular simulation toolkit originating from * * Simbios, the NIH National Center for Physics-Based Simulation of * * Biological Structures at Stanford, funded under the NIH Roadmap for * * Medical Research, grant U54 GM072970. See https://simtk.org. * * * * Portions copyright (c) 2008 Stanford University and the Authors. * * Authors: Peter Eastman * * Contributors: * * * * This program is free software: you can redistribute it and/or modify * * it under the terms of the GNU Lesser General Public License as published * * by the Free Software Foundation, either version 3 of the License, or * * (at your option) any later version. * * * * This program is distributed in the hope that it will be useful, * * but WITHOUT ANY WARRANTY; without even the implied warranty of * * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * * GNU Lesser General Public License for more details. * * * * You should have received a copy of the GNU Lesser General Public License * * along with this program. If not, see . * * -------------------------------------------------------------------------- */ #include "CudaPlatform.h" #include "CudaStreamFactory.h" #include "CudaStreamImpl.h" #include "openmm/OpenMMException.h" #include "openmm/internal/OpenMMContextImpl.h" #include "kernels/gputypes.h" using namespace OpenMM; StreamImpl* CudaStreamFactory::createStreamImpl(std::string name, int size, Stream::DataType type, const Platform& platform, OpenMMContextImpl& context) const { CudaPlatform::PlatformData& data = *static_cast(context.getPlatformData()); if (name == "particlePositions") { float padding[] = {100000.0f, 100000.0f, 100000.0f, 0.2f}; return new CudaStreamImpl(name, size, type, platform, data.gpu->psPosq4, 4, padding, data.gpu); } if (name == "particleVelocities") { float padding[] = {0.0f, 0.0f, 0.0f, 0.0f}; return new CudaStreamImpl(name, size, type, platform, data.gpu->psVelm4, 4, padding, data.gpu); } if (name == "particleForces") { float padding[] = {0.0f, 0.0f, 0.0f, 0.0f}; return new CudaStreamImpl(name, size, type, platform, data.gpu->psForce4, 4, padding, data.gpu); } switch (type) { case Stream::Float: case Stream::Double: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Float2: case Stream::Double2: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Float3: case Stream::Double3: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Float4: case Stream::Double4: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Integer: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Integer2: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Integer3: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); case Stream::Integer4: return new CudaStreamImpl(name, size, type, platform, 1, data.gpu); } throw OpenMMException("Tried to create a Stream with an illegal DataType."); }