Annotation of /trunk/linden/indra/newview/llwind.cpp
Parent Directory
| 
Revision Log
Revision 61 - (view) (download)
| 1 : | mjm | 57 | /** |
| 2 : | * @file llwind.cpp | ||
| 3 : | * @brief LLWind class implementation | ||
| 4 : | * | ||
| 5 : | * $LicenseInfo:firstyear=2000&license=viewergpl$ | ||
| 6 : | * | ||
| 7 : | * Copyright (c) 2000-2008, Linden Research, Inc. | ||
| 8 : | * | ||
| 9 : | * Second Life Viewer Source Code | ||
| 10 : | * The source code in this file ("Source Code") is provided by Linden Lab | ||
| 11 : | * to you under the terms of the GNU General Public License, version 2.0 | ||
| 12 : | * ("GPL"), unless you have obtained a separate licensing agreement | ||
| 13 : | * ("Other License"), formally executed by you and Linden Lab. Terms of | ||
| 14 : | * the GPL can be found in doc/GPL-license.txt in this distribution, or | ||
| 15 : | * online at http://secondlifegrid.net/programs/open_source/licensing/gplv2 | ||
| 16 : | * | ||
| 17 : | * There are special exceptions to the terms and conditions of the GPL as | ||
| 18 : | * it is applied to this Source Code. View the full text of the exception | ||
| 19 : | * in the file doc/FLOSS-exception.txt in this software distribution, or | ||
| 20 : | * online at http://secondlifegrid.net/programs/open_source/licensing/flossexception | ||
| 21 : | * | ||
| 22 : | * By copying, modifying or distributing this software, you acknowledge | ||
| 23 : | * that you have read and understood your obligations described above, | ||
| 24 : | * and agree to abide by those obligations. | ||
| 25 : | * | ||
| 26 : | * ALL LINDEN LAB SOURCE CODE IS PROVIDED "AS IS." LINDEN LAB MAKES NO | ||
| 27 : | * WARRANTIES, EXPRESS, IMPLIED OR OTHERWISE, REGARDING ITS ACCURACY, | ||
| 28 : | * COMPLETENESS OR PERFORMANCE. | ||
| 29 : | * $/LicenseInfo$ | ||
| 30 : | */ | ||
| 31 : | |||
| 32 : | // Wind is a lattice. It is computed on the simulator, and transmitted to the viewer. | ||
| 33 : | // It drives special effects like smoke blowing, trees bending, and grass wiggling. | ||
| 34 : | // | ||
| 35 : | // Currently wind lattice does not interpolate correctly to neighbors. This will need | ||
| 36 : | // work. | ||
| 37 : | |||
| 38 : | #include "llviewerprecompiledheaders.h" | ||
| 39 : | #include "indra_constants.h" | ||
| 40 : | |||
| 41 : | #include "llwind.h" | ||
| 42 : | |||
| 43 : | // linden libraries | ||
| 44 : | #include "llgl.h" | ||
| 45 : | #include "patch_dct.h" | ||
| 46 : | #include "patch_code.h" | ||
| 47 : | |||
| 48 : | // viewer | ||
| 49 : | #include "noise.h" | ||
| 50 : | #include "v4color.h" | ||
| 51 : | #include "llagent.h" | ||
| 52 : | #include "llworld.h" | ||
| 53 : | |||
| 54 : | |||
| 55 : | const F32 CLOUD_DIVERGENCE_COEF = 0.5f; | ||
| 56 : | |||
| 57 : | |||
| 58 : | ////////////////////////////////////////////////////////////////////// | ||
| 59 : | // Construction/Destruction | ||
| 60 : | ////////////////////////////////////////////////////////////////////// | ||
| 61 : | |||
| 62 : | LLWind::LLWind() | ||
| 63 : | : mSize(16), | ||
| 64 : | mCloudDensityp(NULL) | ||
| 65 : | { | ||
| 66 : | init(); | ||
| 67 : | } | ||
| 68 : | |||
| 69 : | |||
| 70 : | LLWind::~LLWind() | ||
| 71 : | { | ||
| 72 : | delete [] mVelX; | ||
| 73 : | delete [] mVelY; | ||
| 74 : | delete [] mCloudVelX; | ||
| 75 : | delete [] mCloudVelY; | ||
| 76 : | } | ||
| 77 : | |||
| 78 : | |||
| 79 : | ////////////////////////////////////////////////////////////////////// | ||
| 80 : | // Public Methods | ||
| 81 : | ////////////////////////////////////////////////////////////////////// | ||
| 82 : | |||
| 83 : | |||
| 84 : | void LLWind::init() | ||
| 85 : | { | ||
| 86 : | // Initialize vector data | ||
| 87 : | mVelX = new F32[mSize*mSize]; | ||
| 88 : | mVelY = new F32[mSize*mSize]; | ||
| 89 : | |||
| 90 : | mCloudVelX = new F32[mSize*mSize]; | ||
| 91 : | mCloudVelY = new F32[mSize*mSize]; | ||
| 92 : | |||
| 93 : | S32 i; | ||
| 94 : | for (i = 0; i < mSize*mSize; i++) | ||
| 95 : | { | ||
| 96 : | mVelX[i] = 0.5f; | ||
| 97 : | mVelY[i] = 0.5f; | ||
| 98 : | mCloudVelX[i] = 0.0f; | ||
| 99 : | mCloudVelY[i] = 0.0f; | ||
| 100 : | } | ||
| 101 : | } | ||
| 102 : | |||
| 103 : | |||
| 104 : | void LLWind::decompress(LLBitPack &bitpack, LLGroupHeader *group_headerp) | ||
| 105 : | { | ||
| 106 : | if (!mCloudDensityp) | ||
| 107 : | { | ||
| 108 : | return; | ||
| 109 : | } | ||
| 110 : | |||
| 111 : | LLPatchHeader patch_header; | ||
| 112 : | S32 buffer[16*16]; | ||
| 113 : | |||
| 114 : | init_patch_decompressor(group_headerp->patch_size); | ||
| 115 : | |||
| 116 : | // Don't use the packed group_header stride because the strides used on | ||
| 117 : | // simulator and viewer are not equal. | ||
| 118 : | group_headerp->stride = group_headerp->patch_size; | ||
| 119 : | set_group_of_patch_header(group_headerp); | ||
| 120 : | |||
| 121 : | // X component | ||
| 122 : | decode_patch_header(bitpack, &patch_header); | ||
| 123 : | decode_patch(bitpack, buffer); | ||
| 124 : | decompress_patch(mVelX, buffer, &patch_header); | ||
| 125 : | |||
| 126 : | // Y component | ||
| 127 : | decode_patch_header(bitpack, &patch_header); | ||
| 128 : | decode_patch(bitpack, buffer); | ||
| 129 : | decompress_patch(mVelY, buffer, &patch_header); | ||
| 130 : | |||
| 131 : | |||
| 132 : | |||
| 133 : | S32 i, j, k; | ||
| 134 : | // HACK -- mCloudVelXY is the same as mVelXY, except we add a divergence | ||
| 135 : | // that is proportional to the gradient of the cloud density | ||
| 136 : | // ==> this helps to clump clouds together | ||
| 137 : | // NOTE ASSUMPTION: cloud density has the same dimensions as the wind field | ||
| 138 : | // This needs to be fixed... causes discrepency at region boundaries | ||
| 139 : | |||
| 140 : | for (j=1; j<mSize-1; j++) | ||
| 141 : | { | ||
| 142 : | for (i=1; i<mSize-1; i++) | ||
| 143 : | { | ||
| 144 : | k = i + j * mSize; | ||
| 145 : | *(mCloudVelX + k) = *(mVelX + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + 1) - *(mCloudDensityp + k - 1)); | ||
| 146 : | *(mCloudVelY + k) = *(mVelY + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + mSize) - *(mCloudDensityp + k - mSize)); | ||
| 147 : | } | ||
| 148 : | } | ||
| 149 : | |||
| 150 : | i = mSize - 1; | ||
| 151 : | for (j=1; j<mSize-1; j++) | ||
| 152 : | { | ||
| 153 : | k = i + j * mSize; | ||
| 154 : | *(mCloudVelX + k) = *(mVelX + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k) - *(mCloudDensityp + k - 2)); | ||
| 155 : | *(mCloudVelY + k) = *(mVelY + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + mSize) - *(mCloudDensityp + k - mSize)); | ||
| 156 : | } | ||
| 157 : | i = 0; | ||
| 158 : | for (j=1; j<mSize-1; j++) | ||
| 159 : | { | ||
| 160 : | k = i + j * mSize; | ||
| 161 : | *(mCloudVelX + k) = *(mVelX + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + 2) - *(mCloudDensityp + k)); | ||
| 162 : | *(mCloudVelY + k) = *(mVelY + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + mSize) - *(mCloudDensityp + k + mSize)); | ||
| 163 : | } | ||
| 164 : | j = mSize - 1; | ||
| 165 : | for (i=1; i<mSize-1; i++) | ||
| 166 : | { | ||
| 167 : | k = i + j * mSize; | ||
| 168 : | *(mCloudVelX + k) = *(mVelX + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + 1) - *(mCloudDensityp + k - 1)); | ||
| 169 : | *(mCloudVelY + k) = *(mVelY + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k) - *(mCloudDensityp + k - 2*mSize)); | ||
| 170 : | } | ||
| 171 : | j = 0; | ||
| 172 : | for (i=1; i<mSize-1; i++) | ||
| 173 : | { | ||
| 174 : | k = i + j * mSize; | ||
| 175 : | *(mCloudVelX + k) = *(mVelX + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + 1) - *(mCloudDensityp + k -1)); | ||
| 176 : | *(mCloudVelY + k) = *(mVelY + k) + CLOUD_DIVERGENCE_COEF * (*(mCloudDensityp + k + 2*mSize) - *(mCloudDensityp + k)); | ||
| 177 : | } | ||
| 178 : | } | ||
| 179 : | |||
| 180 : | |||
| 181 : | LLVector3 LLWind::getAverage() | ||
| 182 : | { | ||
| 183 : | // Returns in average_wind the average wind velocity | ||
| 184 : | LLVector3 average(0.0f, 0.0f, 0.0f); | ||
| 185 : | S32 i, grid_count; | ||
| 186 : | grid_count = mSize * mSize; | ||
| 187 : | for (i = 0; i < grid_count; i++) | ||
| 188 : | { | ||
| 189 : | average.mV[VX] += mVelX[i]; | ||
| 190 : | average.mV[VY] += mVelY[i]; | ||
| 191 : | } | ||
| 192 : | |||
| 193 : | average *= 1.f/((F32)(grid_count)) * WIND_SCALE_HACK; | ||
| 194 : | return average; | ||
| 195 : | } | ||
| 196 : | |||
| 197 : | |||
| 198 : | LLVector3 LLWind::getVelocityNoisy(const LLVector3 &pos_region, const F32 dim) | ||
| 199 : | { | ||
| 200 : | // Resolve a value, using fractal summing to perturb the returned value | ||
| 201 : | LLVector3 r_val(0,0,0); | ||
| 202 : | F32 norm = 1.0f; | ||
| 203 : | if (dim == 8) | ||
| 204 : | { | ||
| 205 : | norm = 1.875; | ||
| 206 : | } | ||
| 207 : | else if (dim == 4) | ||
| 208 : | { | ||
| 209 : | norm = 1.75; | ||
| 210 : | } | ||
| 211 : | else if (dim == 2) | ||
| 212 : | { | ||
| 213 : | norm = 1.5; | ||
| 214 : | } | ||
| 215 : | |||
| 216 : | F32 temp_dim = dim; | ||
| 217 : | while (temp_dim >= 1.0) | ||
| 218 : | { | ||
| 219 : | LLVector3 pos_region_scaled(pos_region * temp_dim); | ||
| 220 : | r_val += getVelocity(pos_region_scaled) * (1.0f/temp_dim); | ||
| 221 : | temp_dim /= 2.0; | ||
| 222 : | } | ||
| 223 : | |||
| 224 : | return r_val * (1.0f/norm) * WIND_SCALE_HACK; | ||
| 225 : | } | ||
| 226 : | |||
| 227 : | |||
| 228 : | LLVector3 LLWind::getVelocity(const LLVector3 &pos_region) | ||
| 229 : | { | ||
| 230 : | llassert(mSize == 16); | ||
| 231 : | // Resolves value of wind at a location relative to SW corner of region | ||
| 232 : | // | ||
| 233 : | // Returns wind magnitude in X,Y components of vector3 | ||
| 234 : | LLVector3 r_val; | ||
| 235 : | F32 dx,dy; | ||
| 236 : | S32 k; | ||
| 237 : | |||
| 238 : | LLVector3 pos_clamped_region(pos_region); | ||
| 239 : | |||
| 240 : | F32 region_width_meters = LLWorld::getInstance()->getRegionWidthInMeters(); | ||
| 241 : | |||
| 242 : | if (pos_clamped_region.mV[VX] < 0.f) | ||
| 243 : | { | ||
| 244 : | pos_clamped_region.mV[VX] = 0.f; | ||
| 245 : | } | ||
| 246 : | else if (pos_clamped_region.mV[VX] >= region_width_meters) | ||
| 247 : | { | ||
| 248 : | pos_clamped_region.mV[VX] = (F32) fmod(pos_clamped_region.mV[VX], region_width_meters); | ||
| 249 : | } | ||
| 250 : | |||
| 251 : | if (pos_clamped_region.mV[VY] < 0.f) | ||
| 252 : | { | ||
| 253 : | pos_clamped_region.mV[VY] = 0.f; | ||
| 254 : | } | ||
| 255 : | else if (pos_clamped_region.mV[VY] >= region_width_meters) | ||
| 256 : | { | ||
| 257 : | pos_clamped_region.mV[VY] = (F32) fmod(pos_clamped_region.mV[VY], region_width_meters); | ||
| 258 : | } | ||
| 259 : | |||
| 260 : | |||
| 261 : | S32 i = llfloor(pos_clamped_region.mV[VX] * mSize / region_width_meters); | ||
| 262 : | S32 j = llfloor(pos_clamped_region.mV[VY] * mSize / region_width_meters); | ||
| 263 : | k = i + j*mSize; | ||
| 264 : | dx = ((pos_clamped_region.mV[VX] * mSize / region_width_meters) - (F32) i); | ||
| 265 : | dy = ((pos_clamped_region.mV[VY] * mSize / region_width_meters) - (F32) j); | ||
| 266 : | |||
| 267 : | if ((i < mSize-1) && (j < mSize-1)) | ||
| 268 : | { | ||
| 269 : | // Interior points, no edges | ||
| 270 : | r_val.mV[VX] = mVelX[k]*(1.0f - dx)*(1.0f - dy) + | ||
| 271 : | mVelX[k + 1]*dx*(1.0f - dy) + | ||
| 272 : | mVelX[k + mSize]*dy*(1.0f - dx) + | ||
| 273 : | mVelX[k + mSize + 1]*dx*dy; | ||
| 274 : | r_val.mV[VY] = mVelY[k]*(1.0f - dx)*(1.0f - dy) + | ||
| 275 : | mVelY[k + 1]*dx*(1.0f - dy) + | ||
| 276 : | mVelY[k + mSize]*dy*(1.0f - dx) + | ||
| 277 : | mVelY[k + mSize + 1]*dx*dy; | ||
| 278 : | } | ||
| 279 : | else | ||
| 280 : | { | ||
| 281 : | r_val.mV[VX] = mVelX[k]; | ||
| 282 : | r_val.mV[VY] = mVelY[k]; | ||
| 283 : | } | ||
| 284 : | |||
| 285 : | r_val.mV[VZ] = 0.f; | ||
| 286 : | return r_val * WIND_SCALE_HACK; | ||
| 287 : | } | ||
| 288 : | |||
| 289 : | |||
| 290 : | LLVector3 LLWind::getCloudVelocity(const LLVector3 &pos_region) | ||
| 291 : | { | ||
| 292 : | llassert(mSize == 16); | ||
| 293 : | // Resolves value of wind at a location relative to SW corner of region | ||
| 294 : | // | ||
| 295 : | // Returns wind magnitude in X,Y components of vector3 | ||
| 296 : | LLVector3 r_val; | ||
| 297 : | F32 dx,dy; | ||
| 298 : | S32 k; | ||
| 299 : | |||
| 300 : | LLVector3 pos_clamped_region(pos_region); | ||
| 301 : | |||
| 302 : | F32 region_width_meters = LLWorld::getInstance()->getRegionWidthInMeters(); | ||
| 303 : | |||
| 304 : | if (pos_clamped_region.mV[VX] < 0.f) | ||
| 305 : | { | ||
| 306 : | pos_clamped_region.mV[VX] = 0.f; | ||
| 307 : | } | ||
| 308 : | else if (pos_clamped_region.mV[VX] >= region_width_meters) | ||
| 309 : | { | ||
| 310 : | pos_clamped_region.mV[VX] = (F32) fmod(pos_clamped_region.mV[VX], region_width_meters); | ||
| 311 : | } | ||
| 312 : | |||
| 313 : | if (pos_clamped_region.mV[VY] < 0.f) | ||
| 314 : | { | ||
| 315 : | pos_clamped_region.mV[VY] = 0.f; | ||
| 316 : | } | ||
| 317 : | else if (pos_clamped_region.mV[VY] >= region_width_meters) | ||
| 318 : | { | ||
| 319 : | pos_clamped_region.mV[VY] = (F32) fmod(pos_clamped_region.mV[VY], region_width_meters); | ||
| 320 : | } | ||
| 321 : | |||
| 322 : | |||
| 323 : | S32 i = llfloor(pos_clamped_region.mV[VX] * mSize / region_width_meters); | ||
| 324 : | S32 j = llfloor(pos_clamped_region.mV[VY] * mSize / region_width_meters); | ||
| 325 : | k = i + j*mSize; | ||
| 326 : | dx = ((pos_clamped_region.mV[VX] * mSize / region_width_meters) - (F32) i); | ||
| 327 : | dy = ((pos_clamped_region.mV[VY] * mSize / region_width_meters) - (F32) j); | ||
| 328 : | |||
| 329 : | if ((i < mSize-1) && (j < mSize-1)) | ||
| 330 : | { | ||
| 331 : | // Interior points, no edges | ||
| 332 : | r_val.mV[VX] = mCloudVelX[k]*(1.0f - dx)*(1.0f - dy) + | ||
| 333 : | mCloudVelX[k + 1]*dx*(1.0f - dy) + | ||
| 334 : | mCloudVelX[k + mSize]*dy*(1.0f - dx) + | ||
| 335 : | mCloudVelX[k + mSize + 1]*dx*dy; | ||
| 336 : | r_val.mV[VY] = mCloudVelY[k]*(1.0f - dx)*(1.0f - dy) + | ||
| 337 : | mCloudVelY[k + 1]*dx*(1.0f - dy) + | ||
| 338 : | mCloudVelY[k + mSize]*dy*(1.0f - dx) + | ||
| 339 : | mCloudVelY[k + mSize + 1]*dx*dy; | ||
| 340 : | } | ||
| 341 : | else | ||
| 342 : | { | ||
| 343 : | r_val.mV[VX] = mCloudVelX[k]; | ||
| 344 : | r_val.mV[VY] = mCloudVelY[k]; | ||
| 345 : | } | ||
| 346 : | |||
| 347 : | r_val.mV[VZ] = 0.f; | ||
| 348 : | return r_val * WIND_SCALE_HACK; | ||
| 349 : | } | ||
| 350 : | |||
| 351 : | |||
| 352 : | void LLWind::setCloudDensityPointer(F32 *densityp) | ||
| 353 : | { | ||
| 354 : | mCloudDensityp = densityp; | ||
| 355 : | } | ||
| 356 : | |||
| 357 : | void LLWind::setOriginGlobal(const LLVector3d &origin_global) | ||
| 358 : | { | ||
| 359 : | mOriginGlobal = origin_global; | ||
| 360 : | } | ||
| 361 : | |||
| 362 : |
| ViewVC Help | |
| Powered by ViewVC 1.0.0 |
