BoatAttack/Assets/Unity Physics Items/Systems/GertsnerSystem.cs

using System.Collections;
using System.Collections.Generic;
using Unity.Burst;
using Unity.Collections;
using Unity.Entities;
using Unity.Jobs;
using Unity.Mathematics;
using Unity.Physics.Authoring;
using Unity.Transforms;
using UnityEngine;
using WaterSystem;
using WaterSystem.Data;

public class GertsnerSystem : JobComponentSystem
{
	private int _waveCount;
	private NativeArray<Wave> waveData;

	protected override void OnCreate()
	{
        base.OnCreate();

        _waveCount = Water.Instance._waves.Length;
		waveData = new NativeArray<Wave>(_waveCount, Allocator.Persistent);
		for (var i = 0; i < waveData.Length; i++)
		{
			waveData[i] = Water.Instance._waves[i];
		}
	}

    protected override void OnDestroy()
    {
        if(waveData != null)
            waveData.Dispose();

        base.OnDestroy();
    }
    protected override JobHandle OnUpdate(JobHandle inputDeps)
	{
		var query = GetEntityQuery(typeof(Translation), typeof(Rotation), typeof(BuoyancyNormal));
		var entities = query.ToEntityArray(Allocator.TempJob);

		var job = new HeightQueryJob
		{
			waveData = waveData,
			time = Time.time,
			entities = entities,
			translations = GetComponentDataFromEntity<Translation>(true),
			rotations = GetComponentDataFromEntity<Rotation>(true),
			normals = GetComponentDataFromEntity<BuoyancyNormal>(false),
			offsetBuffer = GetBufferFromEntity<VoxelOffset>(false),
			heightBuffer = GetBufferFromEntity<VoxelHeight>(false)
		};

		return job.Schedule(entities.Length, 1);
	}

	[BurstCompile]
	public struct HeightQueryJob : IJobParallelFor
	{
		[ReadOnly] public NativeArray<Wave> waveData; // wave data stroed in vec4's like the shader version but packed into one

		[ReadOnly] public float time;

        [DeallocateOnJobCompletion]
        [ReadOnly] public NativeArray<Entity> entities;
		[ReadOnly] public ComponentDataFromEntity<Translation> translations;
		[ReadOnly] public ComponentDataFromEntity<Rotation> rotations;
		[ReadOnly] public BufferFromEntity<VoxelOffset> offsetBuffer;

		[NativeDisableParallelForRestriction]
		public ComponentDataFromEntity<BuoyancyNormal> normals;

		[NativeDisableParallelForRestriction]
		public BufferFromEntity<VoxelHeight> heightBuffer;

		
		public void Execute(int index)
		{
			var entity = entities[index];

			DynamicBuffer<VoxelOffset> offsets = offsetBuffer[entity];
			DynamicBuffer<VoxelHeight> heights = heightBuffer[entity];

			var entityTransform = new RigidTransform(rotations[entity].Value, translations[entity].Value);

			for (int vi = 0; vi < offsets.Length; vi++)
			{
				float3 voxelPos = math.transform(entityTransform, offsets[vi].Value);

				var waveCountMulti = 1f / waveData.Length;
				float3 wavePos = new float3(0f, 0f, 0f);
				float3 waveNorm = new float3(0f, 0f, 0f);

				for (var wave = 0; wave < waveData.Length; wave++) // for each wave
				{
					var thisWave = waveData[wave];

					// Wave data vars
					var pos = new float2(voxelPos.x, voxelPos.z);

					var amplitude = thisWave.amplitude;
					var direction = thisWave.direction;
					var wavelength = thisWave.wavelength;
					var omniPos = thisWave.origin;
					////////////////////////////////wave value calculations//////////////////////////
					var w = 6.28318f / wavelength; // 2pi over wavelength(hardcoded)
					var wSpeed = math.sqrt(9.8f * w); // frequency of the wave based off wavelength
					var peak = 0.8f; // peak value, 1 is the sharpest peaks
					var qi = peak / (amplitude * w * waveData.Length);

					var windDir = new float2(0f, 0f);
					var dir = 0f;

					direction = math.radians(direction); // convert the incoming degrees to radians
					var windDirInput = new float2(math.sin(direction), math.cos(direction)) * (1 - thisWave.onmiDir); // calculate wind direction - TODO - currently radians
					var windOmniInput = (pos - omniPos) * thisWave.onmiDir;

					windDir += windDirInput;
					windDir += windOmniInput;
					windDir = math.normalize(windDir);
					dir = math.dot(windDir, pos - (omniPos * thisWave.onmiDir)); // calculate a gradient along the wind direction

					////////////////////////////position output calculations/////////////////////////
					var calc = dir * w + -time * wSpeed; // the wave calculation
					var cosCalc = math.cos(calc); // cosine version(used for horizontal undulation)
					var sinCalc = math.sin(calc); // sin version(used for vertical undulation)

					// calculate the offsets for the current point
					wavePos.x += qi * amplitude * windDir.x * cosCalc;
					wavePos.z += qi * amplitude * windDir.y * cosCalc;
					wavePos.y += ((sinCalc * amplitude)) * waveCountMulti; // the height is divided by the number of waves 

					if (offsets.Length == 1)
					{
						////////////////////////////normal output calculations/////////////////////////
						float wa = w * amplitude;
						// normal vector
						float3 norm = new float3(-(windDir.xy * wa * cosCalc),
										1 - (qi * wa * sinCalc));
						waveNorm += (norm * waveCountMulti) * amplitude;
					}
				}

				heights[vi] = new VoxelHeight { Value = wavePos };


				if (offsets.Length == 1)
					normals[entity] = new BuoyancyNormal { Value = math.normalize(waveNorm.xzy) };
			}
		}
	}
}