oscura 0.7.0__py3-none-any.whl → 0.10.0__py3-none-any.whl

oscura/analyzers/power/__init__.py

@@ -1,12 +1,13 @@
-"""Power analysis module for Oscura.
+"""Power analysis module for Oscura - IEEE 1459-2010 compliant.
 
 Provides comprehensive power analysis capabilities including:
-- Basic power measurements (instantaneous, average, RMS, peak)
-- AC power analysis (reactive, apparent, power factor)
+- Basic power measurements (instantaneous, average, RMS, peak) per IEEE 1459
+- AC power analysis (reactive, apparent, power factor, harmonics) per IEEE 1459
 - Switching loss analysis for power electronics
+- Conduction loss analysis (MOSFET, IGBT, diode)
 - Safe Operating Area (SOA) analysis
-- Ripple measurement
-- Efficiency calculations
+- Ripple measurement and analysis
+- Efficiency calculations (single and multi-output)
 
 
 Example:
@@ -14,6 +15,9 @@ Example:
     >>> power_trace = instantaneous_power(voltage_trace, current_trace)
     >>> stats = power_statistics(power_trace)
     >>> print(f"Average power: {stats['average']:.2f} W")
+
+References:
+    IEEE 1459-2010: Standard for Power Quality Definitions
 """
 
 from oscura.analyzers.power.ac_power import (
@@ -23,6 +27,7 @@ from oscura.analyzers.power.ac_power import (
     phase_angle,
     power_factor,
     reactive_power,
+    three_phase_power,
     total_harmonic_distortion_power,
 )
 from oscura.analyzers.power.basic import (
@@ -30,39 +35,51 @@ from oscura.analyzers.power.basic import (
     energy,
     instantaneous_power,
     peak_power,
+    power_profile,
     power_statistics,
     rms_power,
 )
 from oscura.analyzers.power.conduction import (
     conduction_loss,
+    diode_conduction_loss,
     duty_cycle_weighted_loss,
     forward_voltage,
+    igbt_conduction_loss,
     mosfet_conduction_loss,
     on_resistance,
+    temperature_derating,
 )
 from oscura.analyzers.power.efficiency import (
     efficiency,
+    efficiency_vs_load,
+    loss_breakdown,
     multi_output_efficiency,
     power_conversion_efficiency,
+    thermal_efficiency,
 )
 from oscura.analyzers.power.ripple import (
     extract_ripple,
     ripple,
+    ripple_envelope,
     ripple_frequency,
+    ripple_harmonics,
     ripple_percentage,
     ripple_statistics,
 )
 from oscura.analyzers.power.soa import (
     SOALimit,
+    SOAViolation,
     check_soa_violations,
     create_mosfet_soa,
     plot_soa,
     soa_analysis,
 )
 from oscura.analyzers.power.switching import (
+    SwitchingEvent,
     switching_energy,
     switching_frequency,
     switching_loss,
+    switching_times,
     total_switching_loss,
     turn_off_loss,
     turn_on_loss,
@@ -70,22 +87,25 @@ from oscura.analyzers.power.switching import (
 
 __all__ = [
     "SOALimit",
+    "SOAViolation",
+    "SwitchingEvent",
     "apparent_power",
     "average_power",
     "check_soa_violations",
-    # Conduction
     "conduction_loss",
     "create_mosfet_soa",
+    "diode_conduction_loss",
     "displacement_power_factor",
     "distortion_power_factor",
     "duty_cycle_weighted_loss",
-    # Efficiency
     "efficiency",
+    "efficiency_vs_load",
     "energy",
     "extract_ripple",
     "forward_voltage",
-    # Basic power
+    "igbt_conduction_loss",
     "instantaneous_power",
+    "loss_breakdown",
     "mosfet_conduction_loss",
     "multi_output_efficiency",
     "on_resistance",
@@ -94,21 +114,24 @@ __all__ = [
     "plot_soa",
     "power_conversion_efficiency",
     "power_factor",
+    "power_profile",
     "power_statistics",
-    # AC power
     "reactive_power",
-    # Ripple
     "ripple",
+    "ripple_envelope",
     "ripple_frequency",
+    "ripple_harmonics",
     "ripple_percentage",
     "ripple_statistics",
     "rms_power",
-    # SOA
     "soa_analysis",
     "switching_energy",
     "switching_frequency",
-    # Switching
     "switching_loss",
+    "switching_times",
+    "temperature_derating",
+    "thermal_efficiency",
+    "three_phase_power",
     "total_harmonic_distortion_power",
     "total_switching_loss",
     "turn_off_loss",

oscura/analyzers/power/basic.py

@@ -232,7 +232,7 @@ def energy(
         power = instantaneous_power(voltage, current)
 
     data = power.data
-    sample_period = power.metadata.time_base
+    sample_period = 1.0 / power.metadata.sample_rate
 
     # Apply time limits
     if start_time is not None or end_time is not None:
@@ -290,7 +290,7 @@ def power_statistics(
         power = instantaneous_power(voltage, current)
 
     data = power.data
-    sample_period = power.metadata.time_base
+    sample_period = 1.0 / power.metadata.sample_rate
 
     # Use scipy trapezoid for stable API across NumPy versions
     from scipy.integrate import trapezoid
@@ -330,7 +330,7 @@ def power_profile(
     """
     power = instantaneous_power(voltage, current)
     data = power.data
-    sample_period = power.metadata.time_base
+    sample_period = 1.0 / power.metadata.sample_rate
 
     if window_size is None:
         # Auto-select: ~1% of total samples or 100, whichever is larger

oscura/analyzers/power/soa.py

@@ -105,7 +105,7 @@ def soa_analysis(
     min_len = min(len(v_data), len(i_data))
     v_data = v_data[:min_len]
     i_data = i_data[:min_len]
-    sample_period = voltage.metadata.time_base
+    sample_period = 1.0 / voltage.metadata.sample_rate
 
     # Select applicable limits based on pulse width
     if pulse_width is None:

oscura/analyzers/power/switching.py

@@ -65,7 +65,7 @@ def _prepare_switching_data(
     v_data = voltage.data[:min_len]
     i_data = current.data[:min_len]
     p_data = power.data[:min_len]
-    sample_period = power.metadata.time_base
+    sample_period = 1.0 / power.metadata.sample_rate
     return v_data, i_data, p_data, sample_period
 
 
@@ -386,7 +386,7 @@ def switching_energy(
         >>> print(f"Switching energy: {e*1e9:.2f} nJ")
     """
     power = instantaneous_power(voltage, current)
-    sample_period = power.metadata.time_base
+    sample_period = 1.0 / power.metadata.sample_rate
     time_vector = np.arange(len(power.data)) * sample_period
 
     # Select time window
@@ -540,7 +540,7 @@ def switching_times(
     min_len = min(len(voltage.data), len(current.data))
     v_data = voltage.data[:min_len]
     i_data = current.data[:min_len]
-    sample_period = voltage.metadata.time_base
+    sample_period = 1.0 / voltage.metadata.sample_rate
 
     v_min, v_max = float(np.min(v_data)), float(np.max(v_data))
     i_min, i_max = float(np.min(i_data)), float(np.max(i_data))